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Tlie  purpose  of  this  stvidy  was  to  improve  measuring  techniques  in 
depth  dose  scanning  of  a  radiation  treatment  field  with  regard  to  in- 
creasing accuracy,,  speed  of  measurement  and  the  design  of  an  easily 
transportable  system.   A  new  measurement  technique  and  a  new  detector 
transport  system  have  been  developed  to  fulfill  these  objectives. 

The  new  measurement  technique,  called  the  repetitive  integral  method, 
estimates  the  input  current  flowing  from  a  radiation  detector  by  repeti- 
tive cycles  of  integrating  and  resetting  over  short  time  intervals  and 
has  been  shown  theoretically  to  be  superior  to  the  standard  rate 
measurement  method  for  the  real-time  dynamic  current  measurements  re- 
quired in  depth  dose  scanning.   Electronic  instrumentation  which 
functions  in  the  repetitive  integral  mode  has  been  developed  and  experi- 
mental results  follow  the  theory.   The  repetitive  integral  measurement 
method  has  the  following  advantages : 

1.  The  effect  of  the  amplifier  time  constant  is  eliminated. 

2.  The  method  is  better  suited  for  measurement  of  pulsed  radiation 
sources. 

3.  Long-term  amplifier  drift  may  be  eliminated. 

ix 


4.  The  integrating  circuit  acts  to  attenuate  higher-frequency  noise, 

5.  Amplifier  thermal-noise  output  is  negligible. 

6.  The  method  is  independent  of  scanning  speed. 

The  new  detector  transport  system  which  uses  stepping  motors  and 
digital  logic  guarantees  the  accuracy  and  reproducibility  of  the  new 
depth  dose  scanning  system.   Digital  techniques  permit  precise  detector 
positioning  to  within  a  small  fraction  of  a  millimeter  and  ensure 
perfect  correlation  between  detector  position  and  amplifier  measurement 
cycles. 

Application  of  the  repetitive  integral  measurement  method  in  con- 
junction with  the  detector  transport  system  enables  accurate  and  repro- 
ducible measurements  to  be  made  at  greater  scanning  speeds  than  previ- 
ously possible.   The  depth  dose  scanning  system  developed  is  completely 
self-contained,  is  housed  in  two  units  which  are  easily  transported  in 
a  small  van  and  provides  a  means  of  rapidly  determining  complete  depth 
dose  distributions  for  a  variety  of  radiation  therapy  equipment. 


CHAPTER  1 


INTRODUCTION 


Shortly  after  the  earliest  therapeutic  applications  of  x-rays, 
it  became  evident  that  some  measure  of  absorbed  energy,  or  dose,  was 
required  to  prevent  over-treatment  and  to  enable  the  compilation  of 
treatment  results  which  would  be  transferable  to  other  patients.   As 
the  energies  of  the  radiation  used  in  x-rny  therapy  increased  the  need 
developed  for  dose  data  at  various  depths  below  the  skin  surface. 
These  doses  below  the  skin  have  come  to  be  called  depth  doses. 

Complete  depth  dose  determination  for  a  radiation  treatment 
field  requires  measurement  at  many  points  throughout  the  field. 
Scanning  devices  which  move  a  radiation  detector  through  a  radiation 
field  have  been  designed  for  obtaining  depth  dose  data  automatically. 
These  scanning  devices  have  been  slow  in  operation,  limited  in 
accuracy  and  so  large  that  transport  from  one  location  to  another 
has  been  difficult. 

The  development  of  computer  techniques  in  x-ray  treatment 
planning  during  the  1960s  has  made  possible  the  computation  of 
absorbed  doses  from  very  complicated  treatment  plans  and  has  provided 
a  means  of  computation  for  individual  patients.   It  is  generally 
agreed  among  radiation  therapists  that  computerized  treatment  planning 
has  improved  radiation  therapy.   Over-treatment  complications  may  be 


avoided  and  new  treatment  plans,  not  previously  used  because  of 
calculation  complexity,  may  be  computed.   The  computer  techniques 
now  available  make  use  of  all  the  known  physical  parameters  concerning 
radiation  absorption  in  tissue,  with  the  limiting  factor  in  dose 
determination  being  the  accuracy  of  the  basic  depth  dose  data  for 
the  radiation  fields  used  in  treatment. 

The  aim  of  this  work  has  been  the  improvement  of  measuring 
techniques  used  in  depth  dose  scanning.  Primary  consideration  has 
been  given  to  increasing  the  accuracy  and  speed  of  measurement  and  to 
designing  a  system  easily  transportable.  This  has  been  accomplished 
with  the  use  of  digital  design  techniques  and  a  new  measuring  method. 

All  testing  has  been  limited  to  cobalt-60  gamma  radiation,  but 
the  results  also  apply  to  the  x-rays  produced  by  other  types  of 
therapy  equipment.  The  radiation  detection  system  has  been  designed 
to  be  compatible  with  the  radiation  intensities  encountered  in 
radiation  therapy  equipment.   All  references  to  scanning  in  later 
chapters  pertain  to  depth  dose  scanning  and  the  term  radiation  always 
refers  to  x  or  gamma  radiation  or  to  both  collectively. 

Reviews  of  the  development  of  radiation  therapy  and  the  various 
methods  of  depth  dose  measurement  are  given  in  Chapter  2.   Chapter  3 
discusses  the  ionization  measureT.ent  method  and  its  limitations  in 
depth  dose  scanning.   A  new  measurement  technique  is  developed  and 
described  in  Chapters  4  and  5.   A  new  electrometer  amplifier  is  dis- 
cussed in  Chapter  6.   Chapter  7  discusses  the  mechanical  scanning 
design  and  the  integration  of  mechanical  and  electronic  systems  into 
the  complete  depth  dose  scanning  system.   Measurement  results,  system 
evaluation  and  conclusions  are  given  in  Chapter  8. 


.CHAPTER  2 


RADIATION  DETECTION  AND  DEPTR  DOSE 


Professor  Wilhelm  Conrad  Roentgen  made  his  historic  discovery  of 
x-rays  on  November  8,  1895  and  published  his  results  soon  thereafter.^ 
Roentgen  had  used  his  hand  in  sone  of  his  experiments  to  test  the  pene- 
trating power  of  x-rays;  a  fact  which  suggested  a  possible  application 
in  medicine.   Just  23  days  following  the  public  announcement  of 
Roentgen's  discovery,  Grubbe,^  a  businessman-inventor,  performed  the 
first  x-ray  therapy  on  a  breast  carcinoma  in  Chicago.   The  treatment 
was  suggested  by  the  physician  treating  Grubbe  for  x-ray  dermatitis 
which  he  had  contracted  in  duplicating  Roentgen's  experiments.   Grubbe 
had  the  dubious  distinction  of  being  the  first  to  use  x-rays  thera- 
peutically and  the  first  to  experience  the  harmful  effects  of  excessive 
x-ray  exposure. 

In  the  early  days  of  radiation  therapy,  skin  reactions  were  ob- 
served in  both  physician  and  patient,  making  it  evident  that  some  unit 
of  dose  for  specification  of  quantity  was  necessary.   The  dosage  was  at 
first  judged  by  the  effect  on  the  skin,  with  erythema  and  epilation  dose 
becoming  the  mainstay  of  clinical  radiotherapy.  No  information  concern- 
ing the  dose  received  at  any  point  under  the  skin  v/as  available.   The 
skin  dose  was  of  primary  importance  as  it  was  the  limiting  factor  in 


therapy.   The  early  x-ray  apparatus  used  no  beam  filtration  and 
operated  at  low  voltage;  so  the  amount  of  soft  radiation  which  was 
readily  absorbed  by  the  skin  was  quite  large. 

Radiation  therapy  is  generally  considered  to  have  come  of  age 
around  1920  when  x-ray  tubes  had  been  improved  so  that  they  could 
tolerate  200  kv  and  the  value  of  beam  filtration  had  been  recognized. 
Beyond  this  point,  progress  was  rapid  in  the  development  of  higher- 
voltage  machines  which  produced  radiation  with  greater  penetration, 
and  the  uses  of  x-radiation  in  therapy  continued  to  grow.  Already  in 
1936  Elward^  reported  that  there  were  over  AOO  different  diseases 
which  were  treated  with  x-ray  therapy.   As  the  uses  of  x-ray  therapy 
have  expanded  and  as  higher-energy  beams  have  been  developed,  there 
have  been  demands  for  more  and  better  dose  determination.   The  skin 
dose  is  no  longer  adequate  as  it  was  in  the  early  days  because  with 
higher- energy  beams  the  doses  below  the  skin  nay  be  the  limiting 
factor  in  therapy.   Doses  at  points  throughout  the  treatment  volume 
are  thus  required.   A  most  interesting  review  of  the  development  of 
radiation  dosimetry  may  be  found  in  papers  by  Glasser'*  and  Victoreen.  ^»^ 

Phantoms  and  Depth  Dose 

As  no  direct  measurement  of  the  radiation  received  at  the  tumor 
is  possible  it  must  be  determined  by  measurement  in  simulated  tissue. 
The  tissue-simulating  material,  called  a  phantom,  should  be  a  material 
which  absorbs  and  scatters  radiation  as  does  human  tissue.   Spiers^ 
pointed  out  that  a  good  simulator  of  tissue  should  have  the  same 


density,  same  atomic  number,  and  contain  the  same  number  of  electrons 
per  gram  as  tissue.   Dry  pressed  wood  of  Masonite  or  hardboard  have 
been  used  as  phantoms,  but  these  are  not  reproducible  from  one  therapy 
center  to  another  and  do  not  absorb  exactly  as  tissue.   Several  other 
reports  »^  have  described  solid  phantom  materials  which  are  tissue 
equivalent.   Pressed  wood  and  other  solid  materials  are  very  limited 
in  their  usefulness  since  the  measuring  probe  may  not  be  moved  from 
one  point  to  another  within  the  medium.   Water  and  soft  tissue  have 
the  same  effective  atomic  number'^  of  7.42,  the  same  density  and  3.34 
X  10^-^  and  3.36  x  10^^  electrons  per  gram  respectively.   The  radiation 
absorption  of  v/atcr  and  soft  tissue  is  thus  essentially  the  same  and 
water  is  usually  chc  phantom  choice  since  it  is  readily  available,  is 
of  constant  consistancy  and  allows  movement  of  the  measurement  probe 
within  the  mediur.i.   Water  phantoms,  however,  require  the  use  of  probes 
sealed  against  moisture. 

Water  phantoms  are  generally  used  as  simulators  of  tissue  when 
the  radiation  dose  is  derived  at  some  depth  below  the  skin.   Measure- 
ments made  in  water  are  usually  specified  as  a  ratio  of  tlie  recorded 
value  at  the  depth  of  interest  to  the  value  at  some  other  fixed  point. 
Johns  and  Cunningham^ ^  define  the  percentage  depth  dose  at  a  depth  d 
as  follows : 


dose  at  depth  d 
Percentage  depth  dose  at  depth  d  =  100  


dose  at  maximum  buildup 


At  energies  less  than  400  kv  the  depth  of  maximum  buildup  is  practically 
zero  and  the  surface  dose  is  used  in  the  definition. 


Although  many  different  instruments  have  been  used  in  radiation 
detection  and  in  depth  dose  determination,  the  operating  principles 
are  few.   The  basic  requirement  of  any  radiation  detector  is  that  it 
interact  with  radiation  so  that  its  response  is  proportional  to  the 
radiation  dose.   Table  1  classifies  the  primary  types  of  radiation 
detectors  according  to  the  detector  medium  and  the  effect  produced 
in  the  detector.   Each  of  the  basic  types  of  instrumentation  will  be 
discussed  briefly,  particularly  for  therapy  applications  in  depth  dose 
determination  with  x  and  gamma  radiation. 

Chemical  Change 

Film 

Photographic  emulsions  have  been  used  for  radiation  dosimetry 
from  the  very  discovery  of  x-rays.   Film  is  a  rather  unique  dosimeter 
as  it  allov^s  dosimetry  of  ionizing  radiation  over  dose  ranges  of 
millirads  to  megarads,  over  time  periods  of  microseconds  to  months  and 
over  areas  of  square  microns  to  square  meters.   The  small  thickness  of 
film  makes  it  particularly  suited  to  radiation  dosimetry-  as  it  can  be 
introduced  into  a  radiation  field  v.'ith  a  minimum  of  spatial  distortion. 
Film  has  a  high  resolution,  an  advantage  which  is  used  in  alpha  and 
beta  particle  autoradiography. 

In  the  dosimetry  of  x  and  gamma  radiation  by  photographic  film, 
the  film  is  exposed  by  secondary  electrons  produced  by  one  or  more  of 
the  primary  processes  (photoelectric  effect,  Compton  effect,  pair 
production)  whereby  x  and  gamma  radiation  interact  x>rith  matter.   Johns 
and  Cunningham^ ^  and  Hendee^^  describe  the  interactions  of  x  and  gamma 


TABLE  1 
PRIMARY  EFFECTS  IN  RADIATION  DETECTORS  AND  DETECTOR  MEDIA 


Instrumentation 


Detector 
Medium 


Radiation 
Effect 


Film 

Chemical  Dosimeter 


Photographic  Emulsion  Chemical  Change 
Solid  or  Liquid 


Calorimeter  Solid  or  Liquid 

Thermoluminescent  Dosimeter  Crystal 
Scintillation  Counter        Crystal  or  Liquid 


Solid  State 
Ionization  Chamber 
Geiger  Counter 
Proportional  Counter 


Semiconductor 

Gas 

Gas 

Gas 


Heat  Production 

Thermolumlnescence 

Light  Production 

Electrical 

(Pulse  or  Current) 


radiation  in  detail.   The  important  fact  concerning  the  interaction  of 
X  and  gannna  radiation  with  the  film  emulsion  is  that  the  mass  attenu- 
ation coefficient  for  the  photoelectric  effect  varies  with  the  cube  of 
the  atomic  number.   The  film  emulsion  is  made  up  largely  of  silver 
bromide  which  has  a  high  atomic  number.   Thus  a  large  amount  of  the 
energy  deposited  in  the  silver  emulsion  is  from  low-energy  photons. 
The  primary  defect  in  film  dosimetry  of  x  and  gamma  radiation  is  the 
dependence  of  sensitivity  on  photon  energy. 

Photographic  methods  with  plates  and  paper  were  used  in  the  early 
1900s  for  dose  determination.   Strips  of  photographic  paper  were  placed 
on  the  patient's  skin  in  the  treatment  area  during  the  time  of  treatment 
with  x-rays.   The  degree  of  blackening  on  the  paper  was  determined  by 
comparison  with  a  standard  chart.   Paper  strips  covered  with  different 
thicknesses  of  aluminum  were  used  to  estimate  depch  dose.   Tliese 
methods  were  unsatisfactory  since  standard  paper,  standard  developer, 
standard  temperature  and  a  standard  scale  were  required.   Also  the 
method  gave  results  only  after  the  radiation  dose  had  been  given.   The 
energy  dependence  of  film  was  especially  important  at  that  tine,  as  no 
filtration  was  used,  and  the  beam  was  composed  largely  of  low-energy 
photons . 

The  energy  dependence  of  film  has  been  overcome  to  some  extent  by 
selective  filtering  and  the  method  has  again  been  used  for  depth  dose 
determinations.   Granke  et  al.-^^  have  used  film  dosimetry  for  determi- 
nation of  tissue  doses  with  2  iMev  x-rays  bv  placing  sheets  of  film 
bet^-T-een  different  sections  of  an  Alderson^^  RAMDO  (RAdiation  ANalog 
Dosimetry)  phantom,  a  tissue  equivalent  phantom  containing  a  human 
skeleton.   Isodose  curves  of  depth  dose  values  for  different  radiation 


fields  were  also  produced  with  film  dosimetry.   Loevinger  and  Spira^^ 
used  film  dosimetry  in  radium  therapy  and  developed  a  method  combining 
multiple  films  to  determine  the  dose  for  a  total  treatment. 

Even  under  the  best  of  conditions,  film  dosimetry  presents 
significant  problems.   Irregular  emulsion  thickness  contributes  to 
varying  sensitivity  over  different  portions  of  the  film.   Variations  in 
film  processing  and  storage  conditions  from  one  laboratory  to  another 
make  it  difficult  to  exchange  results  with  other  researchers.   For 
these  reasons  and  the  problem  of  energy  dependence  film  dosimetry  Is 
not  the  method  of  choice  when  accurate  results  are  required.   Film 
dosimetry  has  been  most  useful  in  personnel  monitoring  where  highly 
accurate  results  are  not  required.   Film  dosimetry  as  a  method  of  depth 
dose  determination  for  x  and  gamma  radiation  has  been  replaced  with 
more  reliable  and  easier  techniques. 

Cliemical  Dosimeters 

With  chemical  dosimeters  the  radiation  dose  is  determined  by  the 
chemical  change  in  a  particular  medium.   Roentgen^  in  his  first  report 
of  the  discovery  of  x-rays  called  attention  to  chemical  reactions 
created  by  the  rays.   Very  early  dosimetry  used  crude  chemical  dosime- 
ters which  changed  color  on  exposure  to  the  soft  x-rays  then  used  in 
therapy  and  diagnosis.   In  1902  small  discs  of  fused  potassium  chloride 
and  sodium  carbonate  which  became  discolored  on  radiation  exposure  were 
used  and  various  colorations  were  related  to  erythema  doses.   Barium 
platinocyanide  formed  the  basis  of  a  dosimeter  developed  in  190A  and 
was  used  extensively  for  many  years.   Again  a  color  change  indicated 
the  degree  of  x-ray  exposure.   This  system  was  highly  energy  dependent 
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because  of  the  high  atomic  number  of  the  dosimeter  material.   These 
early  chemical  systems  were  rather  insensitive  and  required  exposures 
of  1,000  roentgens  or  more  to  produce  visible  changes.   Due  to  these 
problems  and  the  rapid  growth  of  the  ionization  method,  interest  in 
chemical  dosimetry  waned,  but  was  rekindled  with  the  advent  of  atomic 
weapons,  nuclear  reactors  and  the  use  of  high-intensity  radiation 
sources.  These  developments  created  a  need  for  dosimeters  for  measuring 
large  doses  of  x  and  gamma  radiation  and  have  provided  an  impetus  to 
further  research  in  the  area  of  chemical  dosimetry. 

The  Fricke  dosimeter,  described  by  Fricke  and  Morse, ^   is  a 
ferrous-ferric  sulfate  system  whose  principal  drawbacks  are  a  lack  of 
sensitivity  and  a  dependency  on  linear  energy  transfer  of  the  radiation. 
The  Fricke  dosimeter  is  generally  considered  the  best  for  4-40  kilorad 
doses  of  X  and  gamma  radiation.   It  has  been  accepted  as  a  standard  in 
radiation  chemistry  and  it  may  be  used  to  measure  dose  in  absolute  units 
to  within  1-2  percent.   Methods  of  preparing  a  Fricke  dosimeter  are 
given  in  the  National  Bureau  of  Standards  Handbook  85.^   The  Fricke 
dosimeter  has  been  recommended  for  calibration  of  high-energy  electron 
beams  used  in  radiation  therapy. ^^ 

Chemical  dosimetry  is  now  used  principally  for  measuring  higher 
doses  than  those  encountered  in  radiation  therapy  but  has  found  some 
clinical  application  in  determining  the  doses  received  by  interior  body 
structures,  since  the  dosimeter  may  be  inserted  in  body  cavities.   An 
aqueous  solution  enclosed  in  a  suitable  container  is  excellent  for 
biological  dosimetry  since  it  closely  approximates  the  density  and 
atomic  composition  of  biological  tissues.   Chemical  dosimeters  have 
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been  particularly  useful  for  x  and  gamma  ray  dosimetry  in  animals 
exposed  in  atomic  field  tests. 

Stenstrom  and  Lohmann^^»^°  used  a  chemical  dosimeter  derived  from 
methylene  blue  for  determination  of  depth  doses  but,  in  general, 
chemical  dosimeters  have  found  very  little  use  for  depth  dose  measure- 
ments in  radiation  therapy  primarily  due  to  their  insensitivity. 
However,  some  researchers  have  used  the  ferrous-ferric  system. ^^ 
Presently  chemical  dosimeters,  which  may  be  used  for  doses  as  low  as 
AOO  rads,22  ^^g  available  and  thus  are  sensitive  in  the  range  of 
therapy  doses. 

Heat  Production 

Calorimetry  has  been  eniployed  to  measure  the  heat  produced  from 
the  dissipation  of  energy  by  radiation.   A  calorimeter  contains  an 
absorber  specifically  designed  to  convert  all  or  part  of  the  incident 
radiation  into  heat.   The  absorber  contains  a  temperature-sensitive 
element  and  is  thenaally  insulated.   Calorimetry  permits  measurements 
of  radiation  energy  in  teois  of  fundamental  energy  units  thus  providing 
an  absolute  basis  for  evaluation  and  comparison  of  experiments  or 
treatments. 

Calorimeters  have  for  many  years  been  applied  to  many  different 
problems  but  their  use  with  ionizing  radiation  has  been  limited  due 
to  the  small  magnitude  of  energy  involved.   The  calorimetric  method 
has  a  primary  advantage  over  other  dosimetric  methods  as  it  measures 
the  energy  deposited  directly  whereas  other  methods  use  secondary 
processes.   Calorimeters  also  have  some  disadvantages.   Their  sensitivity 
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is  adequate  for  radiation  intensities  and  dose  rates  used  in  diagnostic 
and  therapy  but  is  insufficient  for  use  in  radiation  protection.  Also 
in  some  types  of  radiation  processes  energy,  which  should  be  measured, 
may  be  lost  and  in  soiae  processes  energy  is  gained  from  sources  other 
than  the  radiation  beam.   Scattering  of  radiation  out  of  the  calorimeter, 
photonuclear  reaction  or  the  emission  of  a  penetrating  neutron, 
bremsstrahlung  production  by  secondary  electrons  and  endotheraic  or 
exothermic  chemical  reactions  are  such  processes. 

Calorimetry  was  used  by  early  investigators  in  radiation  measure- 
ment.  Curie  and  Laborde^^  used  the  method  in  the  first  determination 
of  the  rate  of  energy  release  by  radium  disintegration.   Stahel^** 
employed  the  calorimetric  method  in  measuring  the  energy  absorbed  by 
water  on  exposure  to  x  and  gamma  radiation.   The  skin  erythema  dose  in 
humans  was  related  to  these  measurements,  given  in  ergs  per  cubic 
centimeter.   Genna  and  Laughlin^^  demonstrated  the  determination  of 
absorbed  doses  in  radiation  therapy  with  cobalt  by  the  calorimetric 
method.   Rollo  et  al.^°  have  used  calorimetry  for  measurement  of  total 
energy  deposition  of  Grenz  rays.   The  calorimetric  method  has  not  played 
a  large  role  in  depth  dose  measurement  primarily  because  the  ionization 
method  provides  a  simpler  way  of  measurement.   As  more  interest  is 
expressed  in  absolute  units  the  calorimeter  will  surely  play  a  role  as 
a  primary  standard  for  these  units. 

Thermo luminescence 

Thei-moluninescent  dosimeters  have  been  classified  by  Fowler  and 
Attix   as  solid  state  integrating  dosimeters.   This  broad  classification 
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includes  photographic  film,  solid  scintillators  and  plastic  films 
containing  dye.   Semiconductor  devices  could  also  be  included  in  this 
classification  but  they  shall  be  discussed  in  a  separate  section. 
Nonelectrical  solid  state  dosimetry  systems  may  be  further  divided  into 
two  categories:   destructive  readout  and  nondestructive  readout,   llie 
nondestructive  readout  systems  are  those  in  which  the  effect  produced 
by  the  radiation  is  not  destroyed  by  measurement  of  the  effect. 
Nondestructive  systems  may  have  their  radiation  effect  measured  any 
number  of  times.  Destructive  readout  systems  may  have  their  radiation 
effect  measured  only  once  as  the  method  required  in  reading  destroys 
the  effect. 

Nondestructive  systems  include  glasses,  plastics  and  dyes  which 
change  color  upon  radiation  exposure.   Other  nondestructive  effects 
are  radiophotoluminescence  of  phosphate  glass,  degradation  of 
luminescence  of  anthracene  and  other  crystals  and  a  change  in  electron 
spin  resonance  in  alanine.   Most  of  the  nondestructive  system  materials 
are  useful  only  for  doses  greater  than  10*^  rads  but  some  dyes,  phosphate 
glass  and  alanine  are  suitable  for  doses  in  the  range  given  in  radiation 
therapy  (100  to  AOO  rads).   A  dyed  polymethyl  methacrylate,  called 
Perspex  Red  AOO,*  which  becomes  black  on  irradiation  is  commercially 
available.   It  has  been  found  useful  for  x-ray  and  electron  beam  dose 
for  depth  dose  measurements  in  x-ray  contact  therapy  where  the  dose 
rate  is  very  high.   Phosphate  glass  has  a  high  effective  atomic  number 
and  is  thus  quite  energy  dependent.   Bradshaw  et  al.^°  used  powdered 


*May  be  obtained  from  Imperial  Chemical  Industries  Ltd.,  Plastics 
Division,  Welw>-n  Garden  City,  Herts,  England. 
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alanine,  an  amino  acid,  for  depth  dose  measurements  in  animals  and 
satellites.   This  method  is  difficult  due  to  the  very  small  electron 
spin  resonance  signals  produced.  Nondestructive  readout  systems  for 
radiation  dosimetry  are  in  general  not  suitable  for  routine  depth 
dose  measurements  in  the  range  of  therapy  doses. 

Destructive  readout  systems  include  the  thermoluminescent 
materials  and  infrared-stimulated  luminescence.  Many  of  these  systems 
are  suitable  for  measurement  from  below  10  rads  to  lO**  rads.   The 
thermoluminescent  materials  have  had  the  greatest  impact  of  the  two 
systems  on  radiation  dosimetry  methods. 

Many  materials  are  thermoluminescent  but  only  a  few  have  been 
found  useful  for  therapeutic  dosimetry  applications.   The  ideal  thermo- 
luminescent phosphor  must  have  a  strong  light  output  and  be  able  to 
hold  trapped  electrons  for  some  period  of  time  at  a  temperature  suit- 
able for  the  particular  dosimetry  application.   The  three  most  impor- 
tant phosphors  used  in  radiation  dosimetry  in  the  range  of  therapy 
doses  and  lower  arc  calcium  sulfate,  calcium  fluoride  and  lithium 
fluoride. 

Calcium  sulfate  was  used  in  depth  dose  measurement  in  body  cavi- 
ties of  experimental  animals  as  early  as  1954.   This  phosphor  has  the 
remarkable  ability  to  measure  exposures  in  the  microroentgen  range. 
Its  main  disadvantage  is  significant  "fading"  at  room  temperature,  a 
factor  which  apparently  is  dependent  upon  the  method  used  in  phosphor 
preparation. 

Calcium  fluoride  has  a  linear  response  to  gamma  radiation  from 
a  few  milliroentgens  to  500  roentgens.   This  phosphor  does  not  have 
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significant  problems  of  "fading"  but  its  response  is  energy  dependent 
at  energies  below  0.1  Mev. 

Lithium  fluoride  dosimeters  may  be  found  in  a  number  of  shapes 
and  sizes  and  are  widely  used  in  clinical  dosimetry.   Lithium  fluoride 
dosimeters  are  approximately  tissue  equivalent  and  are  of  sizes  that 
may  be  placed  in  body  cavities  or  in  phantoips  for  depth  dose  studies. 
Vacirca  c t  a 1 . ^ ^  have  reported  a  f ilm-thcrmoluminescence  system  for 
determining  body  doses  in  diagnostic  radiography.   Thermoluminescent 
dosimeters  have  been  used  by  Puite  et  al.^^  for  the  purpose  of 
intercalibration  of  x-ray  units  at  many  different  institutions  so  that 
experimental  results  may  be  transferrable.   The  dosimeters  were  ideal 
for  this  purpose  as  they  could  be  sent  from  one  institution  to  another 
through  the  postal  service. 

The  search  for  other  useful  thermoluminescent  materials  continues 
with  the  hope  that  dosimeter  "fading,"  which  is  currently  a  problem, 
may  be  overcome. 

Dixon  and  Watts^^  have  reported  on  barium  fluoride  as  a  thermo- 
luminescent dosimeter  which  exhibits  good  characteristics,  but  "fading" 
is  still  a  serious  problem.   Jayachandran^^  reports  that  lithium 
borate  is  more  tissue  equivalent  than  other  thermoluminescent  dosime- 
ters and  should  be  evaluated  for  clinical  use.   Scarpa^ ^  has  suggested 
that  beryllium  oxide  may  be  a  useful  thermoluminescent  material. 

Light  Production 

One  of  the  earliest  means  of  radiation  detection  was  by  scin- 
tillation counting.  Rutherford  in  his  experiments  with  alpha  particles 
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used  a  zinc  sulfide  crystal  as  a  primary  detector  and  visibly  observed 
the  light  flashes  produced  by  the  alpha  particles  striking  the  crystal. 
This  was  very  tedious  and  with  the  development  of  gas-filled  detectors 
the  scintillation  method  fell  into  disuse.   Interest  in  the  method  was 
revived  in  the  late  1940s  and  phenomenal  development  has  taken  place. 

Many  solid  and  liquid  substances  emit  light  when  exposed  to  x, 
alpha,  beta  or  gamma  radiation.   Photomultiplier  tubes  have  made  it 
possible  to  detect  light  flashes  that  correspond  to  the  absorption  of 
a  single  quantum  or  a  single  charged  particle  and  stimulated  the  search 
for  better  scintillators.   Scintillators  may  be  classified  as  organic 
crystals  of  the  hydrocarbon  type,  inorganic  crystals  and  powders  of 
the  alkali  halide.and  zinc  sulfide  types,  liquid  and  plastic  solutions 
In  hydrocarbon  solutes,  noble  gas  types  and  glass  types.  Murray 3**  has 
listed  the  various  scintillators  available  and  describes  the  important 
characteristics  of  each. 

Two  methods  of  measurement  are  used  with  scintillation  detectors. 
With  the  first  method  individual  particles  or  quanta  are  counted.   In 
the  second  method  the  average  DC  current  is  measured  as  the  total 
light  output  of  the  scintillator.   The  DC  output  of  the  photomultiplier 
tube  is  proportional  to  the  rate  at  which  radiation  energy  is  absorbed 
in  the  scintillator  if  the  light  output  of  the  scintillator  is  propor- 
tional to  the  absorbed  energy.   The  inorganic  crystals  are  of  partic- 
ular importance  for  measuring  x  and  gamma  radiation.   The  alkali 
halides  are  most  often  used  since  they  can  be  grown  into  large  crystals 
of  good  transparency,  a  requirement  for  efficient  detection  of  x  and 
gamma  radiation.   The  inorganic  crystals  are  activated  by  addition  of 
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certain  metals  to  make  them  highly  efficient  scintillators.   Sodium 
iodide  activated  vrith  thallium  is  the  most  useful  scintillator  for 
X  and  gamma  radiation  detection. 

For  purposes  of  x  and  gamma  ray  dosimetry  it  is  desirable  to  have 
a  scintillator  whose  output  per  roentgen  is  constant  over  the  entire 
range  of  energies  for  which  the  detector  is  to  be  used.   The  output 
per  roentgen  is  a  comparison  between  the  scintillator  measurement  and 
the  results  obtained  from  the  standard  free-air  ionization  chamber. 
All  inorganic  scintillators  contain  elements  with  atomic  numbers 
higher  than  that  of  air  and  thus  have  a  greater  output  per  roentgen 
at  low  energies  due  to  the  photoelectric  effect.   At  high  energies 
a  larger  output  is  expected  due  to  pair  production.  With  organic 
crystals  the  output  per  roentgen  at  low  energies  is  determined  by  the 
photoelectric  absorption  of  carbon  which  is  less  than  with  air.  The 
output  per  roentgen  of  a  scintillator  will  be  constant  only  if  it  has 
the  same  effective  atomic  number  as  air,  that  is,  the  scintillator  must 
be  air  equivalent.   The  most  efficient  scintillator  with  the  closest 
effective  atomic  number  to  air  is  anthracene  with  an  effective  atomic 
number  of  5.8.   The  effective  atomic  number  for  air  is  7.64.  Plastic 
and  liquid  scintillators  have  atomic  numbers  similar  to  air  and  thus 
should  be  suited  for  biological  dosimetry.   The  inorganic  scintillators 
such  as  sodium  iodide  may  be  used  to  determine  depth  dose  for  high-ener- 
gy radiation  but  the  effects  of  low-energy  scatter  must  be  considered. 
Simple  pulse-counting  methods  cannot  provide  a  useful  measure  of  x 
or  gamma  ray  exposure  unless  pulse-height  discrimination  is  used. 
Scintillation  detectors  may  be  used  for  direct  measurement  of  depth 
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dose  in  body  cavities  but  calibration  is  required  with  radiation  of 
the  same  energy. 

Electrical 


Solid  State 

Electrical  solid  state  detectors  or  semiconductor  detectors  are 
solid  state  analogs  of  the  ionization  chamber.   The  ionizing  particle 
or  X  or  gamma  ray  reacts  with  the  sensitive  volume  of  the  detector  to 
produce  ionization.   The  ionization  produces  a  conductivity  through 
the  solid  which  is  a  function  of  absorbed  dose  rate  in  the  material. 
In  some  materials ' the  induced  conductivity  change  is  permanent  due  to 
radiation  damage  and  is  thus  a  measure  of  total  absorbed  dose. 

A  semiconductor  consists  of  a  p-n  junction  where  the  p-rcgion  is 
deficient  in  electrons  and  the  n-region  has  an  excess  of  electrons. 
When  an  electric  field  is  connected  across  the  semiconductor  the 
region  in  the  vicinity  of  the  junction  becomes  depleted  of  ions.   The 
presence  of  ionizing  radiation  in  this  region  causes  a  current  to  flow 
which  is  a  measure  of  the  amount  of  radiation.   The  size  of  the  voltage 
pulse  produced  is  proportional  to  the  energy  lost  in  the  detector  by 
the  incident  radiation. 

The  energy  required  to  produce  an  ion  pair  in  most  gases  is  30-40 
ev  whereas  it  is  only  3.5  ev  for  a  silicon  semiconductor.   Thus  many 
more  ions  are  produced  in  a  semiconductor  detector  than  in  a  gas 
detector  of  the  same  mass  for  the  same  amount  of  energy  absorbed.   This 
means  that  the  pulse  from  a  semiconductor  detector  is  about  ten  times 
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larger  than  that  fron  a  gas  detector.   Semiconductor  detectors  have 
a  linear  response,  excellent  resolution,  a  fast  rise  time  and  are 
small  in  size.   They  are  practically  100  percent  efficient  for  partic- 
ulate radiation  but  are  much  less  sensitive  to  x  and  gamma  radiation 
due  to  the  small  depletion  region.   Thicker  depletion  layers  may  be 
obtained  by  applying  a  higher  bias  voltage  or  by  the  use  of  silicon 
with  higher  specific  resistance.   Increased  bias  voltage  increases 
noise  and  thus  methods  of  increasing  the  resistivity  of  silicon  have 
been  developed.   In  the  "lithium  drifted"  silicon  detector  lithium  is 
diffused  into  the  silicon;  a  process  which  increases  the  resistivity 
about  1,000  times.   Thick  junctions  may  be  obtained  with  this  technique 
which  improves  sensitivity  to  x  and  gamma  radiation. 

Semiconductor  detectors  have  been  used  in  clinical  dosimetry  for 
depth  dose  measurements.^^  Silicon  detectors  are  particularly  useful 
for  measuring  depth  dose  distribution  in  bone  since  the  atomic  number 
of  silicon  (14)  is  near  that  of  bone  (13.8). 

Gas  Detectors 

Roentgen  reported  in  his  second  report^^  that  positive  and  negative 
charged  bodies  in  air  were  discharged  when  placed  in  an  x-ray  beam  and 
that  the  rate  of  discharge  was  proportional  to  the  intensity  of  the 
beam.   This  remarkable  observation  has  provided  the  basis  for  gas-filled 
radiation  detectors  used  today. 

All  gas-filled  detectors  operate  on  the  principle  of  collection  of 
the  ions  produced  in  the  gas  by  particulate,  x  or  gamma  radiation.  The 
detectors  consist  of  a  gas-filled  tube  with  a  center  electrode.   The 
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chamber  wall  forms  a  second  electrode.  With  the  proper  potential 
difference  electrons  are  collected  on  the  center  electrode.   The 
operation  of  the  gas  detector  depends  on  the  voltage  difference.   As 
the  voltage  difference  increases,  the  current  collected  increases 
until  it  approaches  asymptotically  the  saturation  current  for  the 
radiation  intensity  being  measured.  All  ions  produced  in  the  chamber 
are  collected  at  the  saturation  voltage  and  thus  the  current  is  pro- 
portional to  the  total  ionization  produced,  which  is  proportional  to 
the  amount  of  radiation  striking  the  chamber.   If  the  voltage  is 
increased  much  beyond  the  saturation  voltage,  secondary  ionization, 
caused  by  collision  of  the  original  ions  formed  by  the  radiation  with 
other  gas  molecules,  begins.   As  soon  as  secondary  ionization  occurs 
a  rapid  multiplication  of  ions  in  the  chamber  takes  place  and  the  total 
current  depends  very  strongly  on  the  applied  voltage.   A  detector 
operating  in  the  region  of  the  saturation  voltage  in  which  secondary 
ionization  is  not  taking  place  is  called  an  ionization  chamber,  or 
simply  an  ion  chamber.   Such  a  chamber  is  said  to  be  working  in  the 
ionization  chamber  region,  indicating  operation  at  saturation  voltage 
with  no  secondary  electron  emission.   The  current  produced  in  an  ion 
chamber  is  very  small  and  requires  amplification  before  measurement. 

The  Geiger  counter,  a  gas -filled  counter  operating  in  the  ion 
multiplication  region,  is  not  a  precise  instrument  for  measurement 
of  absorbed  dose  because  its  response  is  not  directly  proportional 
to  the  energy  absorbed  in  its  sensitive  volume.   It  has  a  detection 
efficiency  of  nearly  100  percent  for  alpha  and  beta  particles  and  is 
sensitive  to  x  and  gamma  radiation,  but  to  a  much  smaller  degree. 
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Geiger  counters  have  been  used  for  the  measurement  of  low-dose  rate 
radiation  fields.   In  this  case  the  recorded  counts  are  related 
indirectly  to  dose  by  considering  the  efficiency  with  wliich  secondary 
electrons  are  ejected  from  the  cathode  walls  as  a  function  of  x  or 
gamma  ray  flux.   These  results  are  then  related  to  the  dose  rate 
measured  with  an  air-wall  ion  chamber.   Counter  response  in  terms  of 
dose  rate  is  a  function  of  the  energy  of  the  x  or  gamma  radiation. 
Geiger  counters  are  not  suitable  for  measurement  of  low- energy  radia- 
tion due  to  the  photoelectric  Interaction  in  the  chamber  wall.   Neither 
are  they  suitable  for  depth  dose  measurement  in  a  scattering  medium 
because  of  the  mixture  of  energies  present  at  different  depths. 
Sinclair^^  has  discussed  the  use  of  Geiger  counters  in  the  measurement 
of  radiation  fields  and  the  problems  involved. 

Proportional  counters  have  not  been  used  for  measurement  of 
radiation  fields  due  to  the  complex  instrumentation  required  to  ensure 
higher-stability  high-voltage  power  supplies  and  amplifiers  of  great 
sensitivity. 

Because  of  the  many  difficulties  in  making  both  direct  and  indirect 
measurement  of  radiation  as  indicated  in  previous  sections,  the  prin- 
cipal method  of  dosimetry  has  been  based  on  the  ionization  of  gases, 
particularly  ionization  of  air.   Although  air  ionization  is  only  an 
indirect  measure  of  energy  absorption  in  other  media  it  has  proven  to 
be  of  the  greatest  value  in  medical  radiology.   The  effective  atomic 
number  of  both  soft  tissue  and  water  is  7.A2  which  is  quite  close  to 
that  of  air  (7.64).   Thus  the  absorption  of  ionizing  radiation  per  gram 
of  air  is  almost  the  same  as  the  absorption  in  tissue  or  water,  and 
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ionization  in  air  offers  an  accurate  and  reproducible  method  of  depth 
dose  dosimetry  in  water  phantoms.   The  unit  of  radiation  exposure,  the 
roentgen,  was  defined  in  1928  at  the  Second  International  Congress  of 
Radiology  in  terms  of  ionization  in  air.   The  definition  of  the  unit 
of  exposure  was  modified  slightly  in  1937^^  but  is  still  essentially 
as  originally  defined.   The  acceptance  of  air  ionization  as  the  means 
of  defining  and  measuring  x  or  gamma  ray  exposure  was  due  largely  to 
the  simplicity  and  reproducibility  of  the  method.  Another  important 
consideration  was  the  fact  that  the  measured  energy  absorption,  al- 
though not  in  tissue,  is  in  an  approximate  tissue-equivalent  material. 
Air  is  a  particularly  suitable  gas  to  use  in  ionization  chambers 
because  of  its  effective  atomic  number,  its  availability  and  its 
constancy  of  composition. 

Air  ionization  chambers  operate  in  the  saturation  voltage  region, 
where  the  output  signal  is  proportional  to  the  energy  deposited  in  the 
chamber  by  the  impinging  radiation.   Air  ionization  chambers  are  used 
routinely  in  equipment  calibration  and  in  depth  dose  measurements  with 
water  phantoms . 

Summary 

Many  methods  of  radiation  detection  are  available  and  all  have 
been  used  at  one  time  or  another  for  determination  of  depth  dose.   The 
principal  method  of  depth  dose  determination  has  been  based  on  ion- 
ization of  air.   This  method  is  simple  to  implement  and  is  easily 
reproducible  from  one  therapy  center  to  another.   The  method  is  partic- 
ularly suited  to  measurement  of  data  to  be  used  in  treatment  planning 
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for  x-ray  therapy  patients  since  the  output  of  an  ionization  chamber 
is  proportional  to  the  ejnergy  absorbed  in  the  air  medium.   Air  has 
essentially  the  same  effective  atomic  number  as  soft  tissue  and  water 
and  thus  absorbs  radiation  as  these  two  media.   It  is  because  of 
these  many  advantages  that  air  ionization  is  the  method  applied  in 
depth  dose  measurement  in  this  research. 

Further  details  of  the  air  ionization  method  and  its  limitations 
as  applied  in  depth  dose  scanning  are  presented  in  the  following 
chapter. 


CHAPTER  3 


THE  IONIZATION  METHOD 


The  physical  phenomenon  of  air  ionization  has  formed  the  basis 
of  the  most  widely  used  and  reliable  method  of  measurement  of  x  and 
gamma  radiation.   The  ionization  of  air  by  x-rays  was  observed  by 
Roentgen  immediately  after  the  discovery  of  the  rays  and  was  discussed 
in  his  early  publications.   The  Curies  also  observed  air  ionization  by 
gamma  rays  from  radium  and  noted  that  the  time  required  to  discharge  a 
gold-leaf  electroscope  was  proportional  to  the  intensity  of  the  radia- 
tion incident  upon  it.   The  electroscope  may  be  considered  as  the  first 
ionization  chamber.  Many  refinements  have  been  made  to  ionization 
chambers  over  the  years,  but  the  principle  of  operation  remains  the 


same. 


Ionization  Chambers  and  Radiation  Mea 


surement 


By  1925  ionization  measurements  had  become  standardized  to  the 
extent  that  agreement  had  been  reached  as  to  what  quantity  should  be 
measured.   Serious  problems  of  wavelength  dependence  due  to  wall  effects 
were,  however,  still  present.   The  early  ionization  chambers  designed 
for  dose  measurements  in  radiation  therapy  were  made  of  various  materi- 
als, usually  metal,  and  were  quickly  discovered  to  be  unsatisfactory 
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because  of  photoelectrons  emitted  from  the  walls.   The  wall  effects 
were  a  particular  problem  because  ionization  chambers  fabricated  at 
different  therapy  centers  did  not  produce  the  same  results  due  to 
different  materials  and  wall  thicknesses.  Most  radiation  physicists 
were  already  at  this  time  using  the  unit  which  v;as  later  defined  as 
the  roentgen  as  a  means  of  comparing  results.   This  unit  called  for 
measurement  of  "the  associated  corpuscular  emission  per  cubic  centi- 
meter of  air  at  standard  temperature  and  pressure,"  which  meant  that 
the  measurement  must  depend  only  on  air  ionization  with  no  ionization 
resulting  from  wall  interactions  Included. 

Duane^'  realized  that  in  order  to  measure  the  ionization  in  1 
cubic  centimeter  of  air  correctly  it  should  be  surrounded  by  a  large 
volume  of  air.   Duane  constructed  a  large  free-air  chamber  but  Glasser'*^ 
pointed  out  that  the  chamber  was  still  too  small  for  correct  measure- 
ment.  Failla"*^  and  several  European  physicists  independently  developed 
several  different  types  of  free-air  chambers  which  were  large  enough. 
Taylor, '*2  after  studying  these  chambers,  developed  the  American  stand- 
ard free-air  chamber.   The  development  of  the  standard  free-air  chamber 
provided  a  means  of  standardized  air  ionization  measurement  and  set  the 
stage  for  the  definition  of  the  unit  of  radiation  exposure  in  1928. 

In  1932  Taylor'*^  reported  agreement  between  the  standard  free-air 
chambers  of  several  national  laboratories  to  within  ±  1  percent.   This 
meant  a  standard  of  radiation  exposure  was  defined  on  an  international 
basis  and  it  then  became  possible  to  exchange,  internationally,  results 
of  radiation  therapy.   A  vast  body  of  knowledge  of  dose-effect  re- 
lationships has  been  built  up  over  the  years  since  the  establishment 
of  the  standard  ionization  measurement. 


Standard  free-air  ionization  chaobers  are  large,  bulky  and  very 
sensitive  instruments  and  are  not  suited  for  routine  radiation  measure- 
ment.  The  purpose  of  a  free-air  ionization  chan±ier  is  to  provide  a 
primary  standard  which  can  be  used  to  calibrate  other  instruments  which 
may  be  used  routinely.   A  dosimeter  used  in  dose  determination  in 
radiation  therapy  must  respond  like  tissue  to  radiation  exposure  and 
must  be  small  enough  to  determine  the  dose  in  a  small  volume.  First 
attempts  to  produce  small  ionization  chambers,  called  "thimble  chambers' 
because  of  their  size  and  shape,  were  not  very  successful.   The 
chambers  were  highly  wavelength  dependent  and  thus  did  not  respond  as 
tissue  to  the  radiation.   Fricke  and  Glasser'*'*  concluded  that  a  cham- 
ber having  a  wall  with  effective  atomic  number  of  7.69  would  have  the 
same  effective  atomic  number  as  air  and  would  thus  produce  the  same 
ionization  per  cubic  centimeter  as  a  standard  free-air  chamber.   This 
conclusion  was  however  not  borne  out  in  practice.   The  first  successful 
thimble  chamber  was  developed  by  Victoreen  in  1927  and  is  discussed  in 
his  review^  of  the  development  of  thimble  chambers.   His  chamber  was 
fabricated  from  carbon  with  an  aluminum  electrode  and  was  designed  to 
give  the  same  output  per  cubic  centimeter  as  a  free-air  chamber.   Each 
chamber  was  balanced  to  the  proper  response  by  calibration  with  a 
free-air  chamber.   The  Victoreen  chamber  vas  attached  to  a  condenser 
which  was  charged  prior  to  the  radiation  measurement.   The  loss  of 
charge  during  irradiation  was  proportional  to  the  dose.   This  type  of 
cavity  or  thimble  chamber  has  been  made  available  commercially*  and 


"This  instrument  is  available  from  Victoreen  Instrument  Company, 
Cleveland,  Ohio. 
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has  been  the  primary  instrument  for  dose  determination  in  radiation 
therapy  for  many  years.   The  instrument  is  commonly  called  today  the 
"condenser  r-meter"  or  simply  the  "r -meter." 

The  design  requirements  of  ionization  chambers  used  in  depth 
dose  scanning  are  not  quite  as  rigid  as  those  discussed  above.   Thimble 
chambers  are  used  to  measure  absolute  valuer  of  radiation  exposure  and 
thus  must  provide  the  same  result  as  would  a  standard  free-air  chamber. 
The  definition  of  depth  dose  given  earlier  involves  the  ratio  of  the 
measured  result  at  the  depth  in  question  to  the  measured  result  at  the 
depth  of  maximum  buildup.   Since  a  ratio  is  involved  it  is  not  essential 
that  the  ionization  chamber  used  in  depth  dose  scanning  produce  the 
same  result  as  would  a  standard  free-air  ionization  chamber.   In  depth 
dose  scanning  both  primary  and  scattered  radiation  are  measured  and 
so  the  ionization  chamber  walls  should  be  made  of  a  material  which  will 
not  produce  significant  wall  effects. 

The  Victoreen  r-meter  used  the  voltage  drop  on  a  condenser  as  a 
measure  of  radiation  exposure.   Later,  ionization  chambers  were  devel- 
oped which  operated  at  constant  potential  supplied  by  a  power  supply 
and  the  ionization  current  was  measured  directly  by  an  electrometer 
during  the  time  of  irradiation.   The  constant  potential  ionization 
chamber  is  the  type  used  in  depth  dose  scanning. 

Measurement  of  Small  Currents 


Ions  produced  by  radiation  interactions  are  collected  in  ionization 
chambers  under  the  influence  of  a  polarizing  voltage.   The  collection 
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of  charge  over  a  period  of  tin,e  constitutes  a  current  flow.   The 
currents  of  interest  In  radiation  measurement  are  very  s^U,  ranging 
fro»  10-  to  10-ie  ,„,„,,  „,^,  ,„-,e  ,,  ,,.,,  ^^^^^^^   ^^^^^^^  ^^ 

those  encountered  in  depth  dose  scanning.  These  small  currents  place 
.cvere  restrictions  upon  circuit  design  and  ma.e  electrometers  manda- 
tory. An  electrometer  is  an  instrument  with  verv  h,>h  ,■ 

wicn  very  high  input  impedance 

and  is  used  to  measure  stall  currents. 

One  Of  the  earliest  electrometers  to  he  used  in  radiation  measure- 
ment was  the  ,uart.  fiher  electrometer,  whose  operation  depends  upon 

the  attractive  or  repulsive  force  existina  Ko^ 

lorce  existing  between  two  conductors,  one 

Of  Which  is  charged  hy  the  signal  to  be  measured  and  the  other  hy  a 
polarising  voltage.   The  original  Victoreen  r-meter  used  this  type,  of 
electrometer.   These  electrometers  are  sensitive  to  shoe,  and  vihration 
and  have  heen  largely  replaced  hy  vacuum  tuhe  electrometers  which  are 
-ore  rugged  and  less  sensitive  to  damage.   Vacuum  tuhe  electrometers 
are  however  prone  to  drift.   In  an  effort  to  overcome  the  drift  proh- 
lem.  the  vihrating-capacitor  electrometer  was  developed.   This  electrom- 
eter^is  extremely  stable  and  highly  sensitive  so  that  currents  do.,  to 

10   ^'  to  10  16  a^pergg  ^^^  ^^   measured.   These  excP^^..^  • 

inese  excellent  instruments 

are  however  rather  hnlVv  =„^  ^ 

duwei  DuiKy  and  expensive. 

solid  state  electrometer  devices  were  developed  slowly  because 
solid  state  components  with  high  input  resistances  were  not  available 
the  field-effect  transistor  has  now  solved  this  problem.,  Mauderli  and 
Bruno-  developed  the  first  solid  state  electrometer  amplifier  using 
^-Id-effect  transistors.   This  amplifier  has  good  sensitivity  and  is 

elose  proximity  to  the  ionization  chamber. 
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The  small  currents  produced  in  ionization  chambers  present  some 
measurement  problems  not  encountered  with  larger  currents.   The  sensi- 
tivity of  an  electrometer  makes  it  particularly  susceptible  to  insu- 
lation difficulties  not  present  in  less  sensitive  current  measurement 
devices.   The  input  to  the  electrometer  must  be  highly  insulated  from 
its  support.   The  insulation  must  be  stable  with  time  and  should  be 
impervious  to  water.   The  electrometer  is  placed  as  near  the  ionization 
chamber  as  possible  and  thus  must  be  ionizatlon-proof .   No  air  spaces 
should  surround  the  electrometer  or  a  second  source  of  ionization  and 
collection  may  occur.   The  electrometer  must  be  surrounded  by  some 
highly  insulating  material  to  exclude  all  air  spaces.   Ceresin  wax  is 
often  used  due  to.  its  high  insulating  properties  and  low  melting  point. 

Measurement  Methods 

Electrical  current  measurements  are  performed  by  one  of  two 
methods,  the  IR-drop  method  or  the  loss-of -charge  method,  also  called 
the  rate-of -drift  method.   Figure  1  illustrates  a  circuit  typical  of 
small  current  measurement.   The  ionization  chamber  operates  at  constant 
potential  supplied  by  a  power  supply,  a  battery  in  this  case.   The 
collection  of  ions  at  the  center  electrode  causes  a  current  flow  which 
is  measured  by  the  electrometer.   The  resistance  R  and  the  capacitance 
C  include  the  total  resistance  and  capacitance  present  in  the  circuit. 
If  all  circuit  capacitance  is  ignored,  that  is ,  C  =  0 ,  the  essential 
characteristics  of  the  IR-drop  method  are  easily  seen.   The  resistor  R, 
in  this  case,  is  the  load  resistor  and  is  an  added  resistance  used  only 
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in  the  IR-drop  circuit.   When  a  polarizing  voltage  sufficient  to  cause 
saturation  in  the  chamber  is  applied,  the  magnitude  of  the  ionization 
current  is  proportional  to  the  rate  of  ion  production  which  in  turn 
is  proportional  to  the  radiation  intensity.   If  I  is  the  constant 
ionization  current,  then  the  voltage  drop,  as  a  function  of  time, 
apparent  at  the  electrometer  is,  by  Ohm's  law, 

V(t)  =  I  R  (1) 

and  the  origin  of  the  name  of  the  method  becomes  quite  clear. 

The  loss-of-charge  method  of  current  measurement  may  be  demon- 
strated from  Figure  1  if  the  total  circuit  resistance  is  made  infinite. 
The  definition  of  capacity  C,  is 


c=-a  (2) 


where  C  is  the  capacity  in  farads,  q  is  the  charge,  in  coulombs,  stored 
in  the  capacitor  and  V  is  the  potential,  in  volts,  at  the  terminals  of 
the  capacitor.   The  current  I  is  defined  as 


1=^  (3) 

d  t 


that  is,  the  rate  of  change  of  charge  with  time.  Differentiating 
equation  (2)  with  respect  to  time  t,  and  substituting  in  equation  (3) 
yields 


iVitl.i  („ 


Integrating  equation  (4)  the  voltage  at  time  t  is 


/ 


V(t)  =  ^  c  dt  =  —  (5) 

0 


assuming  the  capacitor  has  no  charge  at  time  zero  and  the  ionization 
current  I  is  constant.   The  term  loss-of-charge  is  used  because  the 
current  flow  is  in  the  direction  to  neutralize  the  charge  placed  on 
the  chamber  by  the  polarizing  voltage. 

The  actual  current  measurement  circuit  is  neither  of  the  two 
idealized  circuits  just  discussed  but  is  a  combination  of  the  two 
as  shown  in  Figure  1.   In  the  real  case  C  represents  the  total 
capacitance  of  the  circuit  including  that  of  the  ionization  chamber 
and  R  is  the  total  resistance  including  any  load  resistance.   Applying 
Kirchhoff's  law  to  the  circuit  in  Figure  1 


c  ^VUI  .  voo  _  I  =  0  (6) 

dt      R  ^  '' 


and  a  differential  equation  relating  voltage  and  current  results. 
For  constant  current  I  the  solution  is  easily  shown  to  be 

V(t)  =  I  R  [1  -  exp  (-t/RC)]  (7) 

if  the  capacitor  has  no  initial  charge.   If  the  current  is  not  constant 


33 


but  a  function  of  time  the  solution  of  equation  (6)  is 


/; 


V(t)  =  [exp  (-t/RC)/C]    I   exp  (x/RC)  T(x)  dx        (8) 

Jo 

if  the  capacitor  has  no  initial  charge.   Equations  (7)  and  (8)  express 
the  voltage  change  at  the  electrometer. 

Both  the  IR-drop  method  and  loss-of-charge  method  are  used 
routinely  in  radiation  measurement,  depending  on  the  particular 
application.   With  the  idealized  IR-drop  circuit,  it  is  clear  from 
equation  (1)  that  the  voltage  output  of  the  electrometer  is  pro- 
portional to  the  ionization  current  at  any  instant  of  time.   The 
voltage  produced  at  a  particular  time  is  in  no  way  affected  by  any 
changes  in  the  ionization  current,  which  occur  before  that  point  in 
time.   If  the  method  is  used  to  sample  the  ionization  current  at 
prescribed  points  in  time,  the  measured  values  when  plotted  may  not 
come  even  close  to  approximating  the  true  current  function. 

The  loss-of-charge  circuit  in  its  ideal  realization  will  result 
in  a  measurement  value  which  is  the  area  under  a  curve  describing 
the  input  signal.   The  measurement  is  made  over  some  time  interval  as 
opposed  to  an  instantaneous  measurement  with  the  IR-drop  method.   In 
this  case,  the  signal  change  at  every  instant  of  time,  during  the 
measuring  interval,  is  considered  in  the  total  measurement  and  hence 
this  method  is  inherently  the  most  sensitive  and  contains  the  most 
information  about  the  signal.   It  will  later  be  shown  that  noise  is 
negligible  compared  to  the  signal  with  the  loss-of-charge  method. 
The  method  cannot  give  exact  information  concerning  the  signal  at  any 
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instant  of  time  but  the  measurement  value  is  proportional  to  the 
average  of  the  input  signal  over  the  measurement  interval.   The  IR-drop 
method  and  the  loss-of-charge  method  are  often  referred  to  as  the  rate 
method  and  integral  method,  respectively. 

The  two  measurement  methods  have  been  discussed  up  to  this  point 
in  terms  of  the  idealized  circuits.  However,  no  matter  which  of  the 
two  methods  is  chosen  for  a  radiation  measurement,  the  actual  electrom- 
eter circuit  is  that  sho-vn  in  Figure  1  where  both  a  capacitance  and 
resistance  are  present. 

The  capacity  of  the  ionization  chamber  and  other  components  is 
always  present  and  the  circuit  can  never  have  infinite  resistance. 
V.'hen  a  rate  measurement  (IR-drop)  is  desired  a  load  resistor  is  added 
and  the  capacity  is  that  inherent  to  the  circuit.  When  an  integral 
measurement  (loss-of-charge)  is  desired  additional  capacity  is  added 
and  the  load  resistance  is  removed,  although  some  circuit  resistance 
is  still  present.   In  either  case  both  the  R  and  the  C  components  of 
Figure  1  are  present  and  thus  the  circuit  has  an  associated  RC  time 
constant.   Equation  (7)  clearly  shows  the  effect  of  this  time  constant 
for  constant  current.   The  output  from  the  actual  current  measuring 
circuit  reaches  the  value  of  the  idealized  IR-drop  circuit  only  after 
a  certain  time  period  depending  on  the  value  of  the  time  constant.   The 
output  voltage  reaches  99  percent  of  its  final  value  after  a  time  of 

4.6  time  constants.   The  measured  value  differs  from  the  final  value  by 
less  than  IQ-S  ,,i,  ,,,^^  ^   ^.^^   ^^  ^^^^   ^^^^   ^^   ^^^   ^^^^^^^^^   ^^^ 

elapsed.   If  the  input  current  varies  with  time,  equation  (8)  expresses 
the  measurement  result  for  the  circuit  in  Figure  1. 
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It  would  be  desirable  to  have  C  in  Figure  1  equal  to  zero  when 
making  rate  measurements,  in  which  case  equation  (8)  reduces  to 
equation  (1) .  Likewise  it  would  be  desirable  to  have  R  in  Figure  1 
equal  to  infinity  when  making  integral  measurements,  in  which  case 
equation  (8)  reduces  to  equation  (5).   These  ideal  conditions  are 
impossible  to  achieve.   VHicn  a  rate  measurement  is  desired  the  value 
of  the  time  constant  is  made  as  small  as  possible  but  is  limited  by 
the  size  of  the  load  resistor  which  must  be  used  and  the  circuit  stray 
capacitance.  When  an  integral  measurement  is  desired  the  value  of 
the  time  constant  is  made  as  large  as  possible  but  this  is  again 
limited  by  the  size  of  the  capacitor  which  must  be  used  and  the  leakage 
resistance. 

The  rate  method  has  always  been  used  in  depth  dose  scanning 
because  the  value  of  the  ionization  current  must  be  approximated  at 
many  points  in  time.  The  integral  method  has  never  been  used  in  depth 
dose  scanning  due  to  the  nature  of  the  measurement.   The  presence  of 
a  finite  nonzero  time  constant  makes  it  impossible  to  make  measurements 
in  a  radiation  field  instantaneously.   If  the  radiation  field  is  constant 
during  the  time  of  measurement  the  proper  measured  value  may  easily  be 
obtained  by  waiting  an  appropriate  length  of  time  until  the  instruiuent 
has  reached  its  final  value.  However,  if  the  radiation  field  is 
changing,  as  is  the  case  when  depth  dose  scanning  is  performed,  the  use 
of  the  rate  method  may  lead  to  significant  errors  particularly  in  the 
region  where  the  radiation  field  is  changing  most  rapidly. 
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Depth  Dose  Scanning 

In  the  1920s  and  early  1930s  x-ray  equipment  had  reached  a  stage 
of  development  such  that  the  equipment  was  reliable  and  capable  of 
delivering  doses  which  were  reproducible.   The  value  of  beam  filtration 
had  been  recognized  and  x-ray  equipment  was  able  to  operate  at  higher 
voltages.  With  beam  filtration  and  higher-energy  operations,  skin  and 
air  doses  did  not  correlate  with  observed  tumor  reactions  and  thus  the 
need  for  dose  measurements  beneath  the  skin  became  apparent. 

The  first  depth  dose  measurements  were  made  using  Victoreen 
condenser  ionization  chambers  embedded  within  various  tissue-simulating 
phantom  materials.   This  method  was  very  cumbersome  and  extremely  time 
consuming.   With  solid  phantom  materials  it  was  difficult  to  position 
the  chamber  precisely  at  the  desired  depth  making  it  necessary  to  use 
water  phantoms.   The  use  of  water  presented  still  another  problem,  that 
of  waterproofing  the  condenser  ionization  chamber.   Depth  dose  measure- 
ments performed  in  this  manner  were  limited  to  the  central  beam  axis 
because  of  the  great  amount  of  time  required  in  charging  the  chamber, 
waterproofing  it,  placing  it  in  the  phantom  and  reading  the  chamber. 
Central-axis  depth  dose  measurements  did  not  adequately  describe  radia- 
tion fields  used  in  therapy,  particularly  for  isotope  machines  with 
large  penumbra.   The  lack  of  sharp  beam  edges  due  to  large  source  size 
in  isotope  machines  made  it  imperative  that  entire  radiation  distribu- 
tions be  determined  in  the  form  of  isodose  curves. 

Kernp^"  developed  an  automatic  method  of  scanning  a  radiation  field 
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v/ith  a  diract-reading  ionization  chamber.   The  scanner  plotted  isodose 
curves  during  the  measurement  through  a  selsjm  motor  arrangement  and 
a  balancing  circuit.   Kis  device  though  automatic  was  very  slow,  moving 
at  a  speed  of  only  3  centimeters  per  minute.   An  entire  set  of  10 
isodose  curves  for  1  radiation  field  required  about  2  hours.   Several 
other  researchers'*'"^''  later  used  depth  dose  scanners  which  were 
essentially  the  same  as  that  developed  by  Kemp.   In  1954  Mauchel  and 
Johns^  developed  an  automatic  scanner  capable  of  measuring  a  set  of 
isodose  curves  in  10  minutes.   Berman  et  al.^^  reported  a  scanner  using 
a  ratio  circuit  to  eliminate  variations  in  response  due  to  variations 
of  output  from  an  x-ray  source. 

In  1966  Fitzgerald  et  al.^^  reported  the  first  digital  scanner 
capable  of  recording  data  at  discrete  points  during  scanning.   This 
scanner,  rather  than  plot  isodose  curves,  measured  data  at  5  millimeter 
intervals  so  that  a  dose  matrix  was  formed  representing  the  dose  dis- 
tribution.  The  dose  values  were  recorded  automatically  on  punched 
cards,  a  form  suitable  for  use  in  radiation  treatment  planning  computer 
programs.   This  scanner  was  the  first  which  recorded  data  suitable  for 
immediate  computer  processing.   If  required,  standard  isodose  curves 
could  be  produced  from  the  measured  data  by  a  digital  plotter  attached 
to  a  computer.   A  complete  radiation  field  distribution  could  be  re- 
corded in  digital  matrix  form  in  about  15  minutes.   This  scanner  used 
a  solid  state  electrometer  amplifier"*^  and  a  small  constant-potential 
ionization  chamber.   In  1969  tCing^'*  reported  a  similar  scanning  device 
which  produced  isodose  curves  during  the  scan.   Glenn  et  al.^^  developed 
a  computer-controlled  scanner  which  plotted  isodose  curves  during  meas- 
urement and  collected  data  in  digital  form  for  treatment  planning. 
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A  number  of  depth  dose  scanners*  are  now  available  conmercially. 
There  are  many  different  features  available  on  these  scanners  but  they 
all  have  one  thing  in  common,  the  measureinent  method.   All  use  the  rate 
method  of  measurement,  a  distinct  disadvantage. 

Summary 

The  ionization  method  of  radiation  measurement  had  become  the 
method  of  choice  by  1925.   The  development  of  the  standard  free-air 
ionization  chamber  paved  the  way  for  the  definition  of  the  roentgen 
as  the  unit  of  radiation  exposure  and  established  a  means  of  standard 
measurement  so  that  transfer  of  radiation  therapy  results  was  possible 
on  an  international  basis. 

The  currents  produced  in  an  ionization  chamber  on  exposure  to 
radiation  are  small  and  an  electrometer  is  required  to  measure  them. 
Current  measurements  are  made  by  either  the  rate  (IR-drop)  method  or 
the  integral  (loss-of-charge)  method.   The  rate  method  has  always  been 
used  in  depth  dose  scanning  and  the  RC  time  constant  associated  with 
the  measuring  circuit  has  posed  a  problem  in  this  application. 

The  limitations  of  the  rate  method  and  the  development  of  a  new 
measuring  technique  based  on  an  integral  measurement  method  is  discussed 
in  the  following  chapter. 


*Depth  dose  scanners  are  available  from: 

Victoreen  Instrument  Company,  Cleveland,  Ohio 
Artronix  Medical  Systems,  St.  Louis,  Missouri 
SHM  Nuclear  Corporation,  Sunnyvale,  California 
Scanditronix,  Uppsala,  Sweden 


CHAPTER  4 


THE  RATE  MEASUREMENT  AND  THE  REPETITIVE 
INTEGRAL  MEASUREMENT  METHODS 


The  Lateral  Dose  Function 

A  radiation  detector  moving  across  a  typical  s>Tnmetric  teletherapy 
radiation  field  should  produce  a  response,  which  when  plotted  against 
detector  position,  results  in  a  curve  of  the  type  shown  in  Figure  2. 
This  curve  is  representative  of  all  symmetric  radiation  fields  and  is 
the  idealized  response  of  the  detector,  that  is,  no  time  constant  error 
is  present. 

Daniel  and  Wood,^^  in  a  review  of  functions  producing  certain 
basic  shapes,  include  the  function 


h(x)  =  L 


Cx 


A  +  B 


For  A  and  B  positive,  this  function  takes  values  between  zero  and  a 
maximum  of  1/A.   The  range  of  x  required  to  produce  values  from  zero 
to  the  maximum  depends  on  the  values  of  B  and  C,   This  function  may  be 
used  to  approximate  as  closely  as  desired  the  radiation  detector  re- 
sponse as  the  detector  moves  into  a  radiation  field.   By  properly 
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choosing  the  values  of  A,  B  and  C,  the  function  h(x)  can  closely 
approximate  that  portion  of  the  curve  shown  in  Figure  2  corresponding 
to  abscissa  values  ranging  from  0  to  10  centimeters.  The  detector 
response  may  be  approximated  only  as  it  moves  from  outside  the  radi- 
ation field  into  the  field  since  once  h(x)  reaches  the  maximum  value 
of  1/A  it  remains  at  this  value  for  all  larger  values  of  x.   The 
other  one-half  of  the  response  curve,  representing  the  change  as  the 
detector  moves  from  inside  to  outside  the  field,  may  be  approximated 
by  subtracting  a  function  of  the  same  type  as  h(x)  from  h(x). 

A  new  function  was  defined  by  combining  two  functions  of  the  type 
given  by  h{x).  The  new  function,  f (x) ,  may  be  used  to  approximate  the 
entire  detector  response  as  it  moves  across  the  radiation  field.  This 
function  is  given  by 


f  (x)  = ^ for  X  >  0  (9) 

A  +  B^<^-^>   A  -.-  bC(''-^-> 


The  function  f(x)  has  been  defined  as  the  lateral  dose  function  and 
was  developed  for  use  in  mathematical  analyses  of  detector  response. 
The  curve  in  Figure  2  may  be  reproduced  with  this  function  with 
A  =  0.01,  B  =  2.0,  C  =  -lA.O,  D  =  1.0,  and  E  =  6.0. 

Depth  Dose  Scanning  and  the  Rate  Measurement 

Depth  dose  scanners  have  greatly  improved  over  the  years.   The 
mechanical  scanning  apparatus  has  been  improved  in  order  to  increase 
the  scanning  speed,  automate  the  scanning,  and  ensure  the  reproduc- 
ibility of  positioning  of  the  ionization  chamber.   It  is  generally 
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agreed  that  the  present  limitation  on  scanning  speed  is  the  amplifier 
time  constant.   This  limitation  was  recognized  by  the  developers  of 
early  scanners  and  was  responsible  for  the  long  measuring  times . 

In  the  author's^^  experience  of  depth  dose  scanning,  the  amplifier 
time  constant  is  normally  a  problem  in  that  area  of  the  radiation  field 
where  a  large  intensity  gradient  exists.   Figure  2  shows  a  plot  repre- 
sentative of  a  single  lateral  scan  across  a  radiation  field  and  corre- 
sponds to  a  curve  which  would  be  determined  by  an  idealized  rate 
measurement.   It  is  clear  that  a  rapid  change  of  signal  occurs  as  the 
detector  Dioves  into  the  radiation  field.   The  change  from  almost  zero 
signal  to  the  maximum  signal  occurs  within  a  few  millimeters.   A 
reverse  change  occurs  at  the  other  field  edge  as  the  ionization  chamber 
moves  out  of  the  beam.   It  is  in  these  two  regions  that  the  amplifier 
time  constant  presents  a  problem. 

The  curve  in  Figure  2  is  expressed  in  normalized  units  as  a  func- 
tion of  detector  (ionization  chamber)  position.   Hie  curve  represents 
the  voltage  output  of  an  electrometer  capable  of  making  idealized  rate 
measurements  or  the  ionization  current  at  the  electrometer  input.   When 
the  curve  represents  the  ionization  current,  the  response  of  a  nonideal- 
ized  rate-measurement  circuit  to  this  input  current,  given  by  equation 
(9) ,  may  be  mathematically  computed  by  evaluation  of  equation  (8) . 

The  dotted  and  dashed  curves  in  Figure  3  are  the  computed  responses 
of  a  rate-measuring  circuit,  with  a  finite  time  constant,  to  the  input 
current  given  by  equation  (9)  and  correspond  to  scans  across  a  radiation 
field  in  opposite  directions  at  the  same  depth.   Equation  (9)  is  eval- 
uated with  A  =  0.01,  B  =  2.0,  C  =  -14.0,  D  =  1.0,  and  E  =  6.0.   The 
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dotted  and  dashed  curves  are  obtained  by  numeric  integration*  of 
equation  (8)  with  the  time  constant  equal  to  300  milliseconds  and  a 
scanning  speed  of  2.54  centimeters  per  second.   Because  of  the  nonzero 
time  constant,  the  curves  show  displacement  toward  the  direction  of 
the  scan  in  areas  of  rapid  signal  change.   The  curves  correspond  with 
each  other  and  with  the  solid  curve,  which  represents  the  ideal 
response,  in  the  plateau  areas  since  the  input  does  not  change  and 
the  output  reaches  its  final  value  in  these  regions.   If  no  time 
constant  error  were  present  the  scans  made  in  opposite  directions 
would  be  identical  and  would  lie  on  the  solid  curve  in  Figure  3. 

The  deviations  of  the  dotted  and  dashed  curves  from  each  other 
and  from  the  correct  response,  given  by  the  solid  curve,  are  greater 
in  Figure  3  than  would  be  encountered  in  depth  dose  scanning  due  to 
the  long  time  constant  and  fast  scanning  speed  chosen  to  more  clearly 
demonstrate  the  problem  of  the  time  constant.   Figure  4  shows  two 
scans,  in  opposite  directions,  of  a  10  x  10  centimeter  cobalt  treatment 
field  performed  with  the  digital  depth  dose  scanner^-^  previously 
mentioned.   The  time  constant  in  this  case  was  about  100  milliseconds 
and  the  scanning  speed  was  just  over  1  centimeter  per  second.   The 
shifting  effect,  indicative  of  time  constant  error,  is  still  clearly 
evident  which  shows  that  the  rate  measurement  method,  as  applied  in 
depth  dose  scanning,  has  severe  limitations. 


*The  method  of  numeric  integration  was  three-point  quadrature, 
generally  known  as  Simpson's  Rule  and  described  in  all  texts  of 
integral  calculus. 
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Limitations  of  the  Rate  Measurement  Method 

The  accuracy  of  the  rate  measurement  method,  as  applied  in  depth 
dose  scanning,  is  limited  by  the  effect  of  the  nonzero  amplifier 
time  constant.   The  synergistic  effect  of  several  factors  produces 
this  problem.   The  most  important  of  these  factors  are  the  scanning 
speed,  the  magnitude  of  the  measured  currents  and  the  volume  of  the 
ionization  chamber  used. 

The  Scanning  Speed 

The  inaccuracy  introduced  by  the  nonzero  time  constant  may  be 
completely  eliminated  by  sufficiently  reducing  the  scanning  speed. 
If  scanning  is  so  slow  that  the  amplifier  output  reaches  its  final 
value  in  a  small  enough  time  interval  such  that  the  ionization  chamber 
movement  is  insignificant,  the  nonzero  time  constant  ceases  to  effect 
the  measuring  accuracy.   Slow  scanning,  however,  defeats  the  very 
purpose  for  which  depth  dose  scanners  were  developed.   If  adequate 
computer  treatment  planning  is  to  be  done,  radiation  depth  dose  data 
of  several  hundred  different  treatment  fields  are  required.   A  slow 
scanner  would  be  very  im.practical  in  view  of  the  man-hours  required 
for  depth  dose  measurements.   If  a  scanning  ionization  chamber  is  to 
move  no  more  than  0.5  millimeter  during  the  period  of  amplifier  output 
buildup  to  99  percent  of  its  final  value,  the  maximum  scanner  speed 
is  approximately  0.11  centimeter  per  second  when  the  amplifier  time 
constant  is  100  milliseconds.   In  this  case  the  complete  scan  of  one 
radiation  field  would  require  almost  3  hours.   Such  a  long  scanning 
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time  is  impractical  particularily  with  linear  accelerators  since  the 
generating  tubes  have  limited  operational  life. 

The  Current  Magnitude 

The  magnitude. of  the  current  produced  by  the  radiation-induced 
ionization  in  an  ionization  chamber  is  directly  related  to  the 
intensity  of  the  radiation  and  the  volume  of  the  detecting  chamber. 
A  cobalt  teletherapy  unit  typically  has  a  radiation  output  of  about 
80  roentgens  per  minute  at  a  treatment  distance  of  100  centimeters. 
This  radiation  output  would  produce  a  charge  of  80  electrostatic 
units  per  minute  per  cubic  centimeter  in  an  ionization  chamber.   The 
collected  ions  result  in  a  current  I  expressed  by 

I  =  80  esu/min/cc ^   ^_^^  ^  Iq-IO  c/sec/cc 

[3  X  lOy  esu/C]  [60  sec/min) 

•=  4.44  X  10"^°  A/cc 

The  ionization  current  may  be  made  as  large  as  desired,  depending  on 
the  chamber  volume. 

The  Cham.ber  Volume 

The  volume  of  an  ionization  chamber  used  in  depth  dose  scanning 
must  be  kept  as  small  as  possible  to  ensure  good  resolution.   The 
chamber  response  is  a  measure  of  the  average  ionization  occurring  in 
the  chamber  in  that  portion  of  the  radiation  field  covering  the 
sensitive  volume  and  with  a  smaller  chamber  volume  a  smaller  radiation 
gradient  is  present  across  the  chamber.   Cylindrical  ionization 
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chambers  are  most  often  used  in  depth  dose  scanning.   The  cylinder 
axis  is  perpendicular  to  the  direction  of  scanning  and  the  radiation 
beam  axis  so  as  to  minimize  the  radiation  gradient  across  the  chamber. 
The  scanning  chamber  used  in  this  study  has  an  inside  diameter  of 
3.8  millimeters  and  is  26.5  millimeters  in  length.   The  sensitive 
volume  is  0.21  cubic  centimeter.   The  current  delivered  to  the 
electrometer  from  the  ionization  chamber  is  9  x  10~^^  amperes  with  a 
radiation  output  of  80  roentgens  per  minute. 

The  Time  Constant 

The  amplifier  time  constant  i  is  defined  as 


T  =  RC 


where  C  is  the  total  circuit  capacitance  and  R  is  the  feedback 
resistance  in  the  circuit.   The  value  of  C  may  be  1  picofarad  in  a 
typical  rate-measuring  circuit.   With  a  rate  circuit,  the  voltage 
developed  at  the  electrometer  amplifier  is  the  product  of  the  input 
current  in  amperes  and  the  feedback  resistance  in  ohms.   This  means, 
that  with  the  ionization  current  calculated  above,  a  resistor  of  about 
10^^  ohms  is  required  to  produce  a  10  volt  output  at  the  amplifier 
output.   The  required  selection  of  R  completely  determines  the  value 
of  T  since  C  is  fixed.   The  value  of  the  time  constant  under  these 
conditions  is 
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T  =  RC  =  [1  X  lO^^fi  ]  [1  X  10  12  f]  =  0.1  sec  =  100  msec 

This  value  of  t  is  fixed  by  the  desired  output  voltage  and  may  not  be 
reduced. 


Correction  Technique 

Since  the  time  constant  limitation  of  the  rate  measurement  method 
cannot  be  eliminated  in  a  satisfactory  way,  a  means  of  correcting  the 
measured  scan  data  was  developed.   Equation  (6)  is  a  differential 
equation  relating -current  and  voltage  in  an  electrometer  amplifier  and 
has  the  form 


dt      R 


This  may  be  rewritten  as 


l(t)  R  =  V(t)  +  RC  ^Ktl  =  v(t)  +  ^  dV(t) 

dt  dt 


This  equation  relates  the  final  output  voltage  (I(t)  R)  to  the 
electrometer  output  voltage  and  its  derivative.   The  derivative  of 
the  output  voltage  may  be  approximated  from  output  voltage  values. 
If  Vi  and  V'i  +  ^  are  two  values  of  the  electrometer  output  voltage 
at  two  sampling  points  separated  by  a  distance  ^x,    the  derivative  of 
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the  output  voltage  midway  between  the  two  measuring  points  may  be 
approximated  by 


Vi  +  1  -  Vi 


AX 

V 


where  v  is  the  scanner  speed. 

The  value  of  the  output  voltage  at  the  point  midway  between  the  two 

sampled  points  may  be  approximated  by 


Vi  +  1  +  V. 


The  approximation  to  the  correct  output  voltage  is  given  by 

f 

Vi  =  I(t)  R  ~  Vj  +  1  +  ^i  +  T  Vi  +  1  -  ^i  (10) 

2  AX 

V 

where  x  is  the  time  constant  of  the  amplifier. 

This  correction  technique  was  applied  to  each  sampled  point  in 
each  depth  dose  scan  performed  with  the  previously  described  digital 
scanner.   The  correction  greatly  reduced  the  effect  of  the  time  constant 
and  produced  scans  which  were  essentially  symmetric  in  both  directions. 
This  correction  technique  was  particularly  attractive  since  the 
measured  data  could  be  corrected  after  measurement  with  no  increase 
in  scanning  time  required.   Over  700  different  radiation  fields  were 
scanned  and  corrected  as  above.   Many  of  the  isodose  curves  plotted 
from  these  data  are  included  in  the  International  Atomic  Energy  Agency 
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Dosimetry  Catalog^^  and  have  had  worldwide  distribution. 

The  correction  technique  provides  only  an  approximation  to  the 
correct  chamber  response  since  the  differential  equation  is  evaluated 
using  discrete  data  points.   This  correction  technique  then  becomes 
less  satisfactory  as  the  scanning  speed  increases.   The  fact  that  the 
rate  measurement  method  is  unsuitable  for  the  rapid  sampling  required 
in  depth  dose  scanning  and  the  need  for  faster  depth  dose  scanners 
has  promoted  a  search  for  a  new  measurement  technique. 

The  Integral  Method 

An  integral  method  of  measurement  is  used  in  radiation  therapy 
when  the  total  dose  administered  over  some  time  period  is  desired. 
The  Victoreen  r-nieter,  developed  in  the  1920s  and  still  the  primary 
instrument  for  equipment  calibration,  is  such  an  integrating  device. 
The  loss  of  charge  on  the  ionization  chamber  capacitor  is  the  measure 
of  total  accumraulated  dose  during  the  time  of  exposure.   The  time  of 
exposure  in  this  case  is  determined  by  a  separate  timing  device. 

Integrating  dosimeters  have  been  used  primarily  in  dose  monitor 
devices  which  turn  off  the  x-ray  unit  after  a  preset  dose  is  delivered. 
These  monitors  have  used  both  electromechanical  and  electronic  methods 
of  integration.   Wheatly^^  and  Farr^^  reported  the  use  of  a  small  motor 
driven  by  an  amplifier  whose  output  voltage  was  proportional  to  the 
ionization  current.   The  motor  speed  varied  linearly  with  applied 
voltage  and  the  number  of  motor  revolutions-  made  in  a  period  of  time 
was  proportional  to  the  integral  of  the  ionization  current.   After 
a  preset  number  of  revolutions  the  x-ray  unit  was  turned  off.   This 
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electromechanical  means  of  digital  integration  has  not  been  widely 
used  because  the  motor  speed  was  not  quite  linear  with  voltage  at  low 
voltage  levels  and  primarily  because  better  techniques  became  available. 

When  integration  is  performed  electronically  the  integrating 
capacitor  is  charged  by  the  ionization  current  and  the  voltage  present 
after  some  time  is  proportional  to  the  integral  of  the  current.   Inte- 
gration may  also  be  accomplished  by  initially  charging  the  capacitor  to 
some  fixed  voltage  and  letting  the  capacitor  be  discharged  by  the 
ionization  current.  The  charge  lost  during  the  measuring  time  is  pro- 
portional to  the  integral  of  the  current.   Both  methods  have  been  used 
in  dose  monitors  to  turn  off  an  x-ray  unit  when  the  voltage  on  the 
capacitor  reached  a  preset  voltage  level.   The  dose  monitor  had  to 
work  over  a  wide  dose  range  from  just  a  few  roentgens  to  several  hundred. 
It  was  difficult  to  cover  the  entire  required  dose  range  wi  th  a  single 
integrating  capacitor  and  thus  a  method  of  digital  integration  was 
developed. 

A  single  integrating  capacitor  may  be  used  in  a  dose  monitor, 
operating  over  a  wide  dose  range,  if  the  capacitor  can  be  recharged 
when  the  voltage  on  the  capacitor  drops  to  a  preset  voltage  or,  with 
the  other  method  of  integrating,  if  the  capacitor  can  be  reset  to  zero 
when  it  reaches  a  preset  voltage.   The  number  of  times  the  capacitor 
goes  through  its  voltage  excursion  and  recharge  or  discharge  cycles  is 
a  digital  measure  of  the  integral  of  the  ionization  current.   Farmer^^ 
developed  a  measuring  circuit  of  this  type  which  used  an  electrostatic 
relay  to  actuate  a  counter  and  at  the  same  time  recharge  the  integrating 
capacitor.   This  device  was  limited  to  measuring  dose  rates  of  less  than 
40  roentgens  per  minute  due  to  the  recharge  time  of  the  capacitor. 
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Watson^ ^  used  a  neon  tube,  which  flashed  after  a  particular  voltage 
built  up  on  the  capacitor,  to  discharge  the  capacitor.   This  method 
was  used  in  a  Geiger-counter  type  personnel  monitor  worn  by  persons 
working  in  radiation  therapy.   This  device  could  be  adjusted  to  operate 
over  a  dose  range  of  0.1  milliroentgen  to  10  roentgens  per  minute. 
Watson   also  developed  a  dose  monitor  for  use  in  radiation  therapy. 
In  his  design,  a  moving-coil  relay  operated  to  reset  the  capacitor. 
The  device  produced  one  count  for  every  10  roentgens  and  could  operate 
at  up  to  30  cycles  per  minute.  Naylor,^^  in  a  similar  circuit,  used 
a  cold  cathode  trigger  tube  to  discharge  the  capacitor.   Kemp  et  al.^** 
and  Kemplay^^  have  used  solid  state  electronics  to  reset  the  capacitor 
with  cycle  repetition  rates  of  100  or  more  per  minute.   Lorenz  and 
Mauderli^^  reported  a  similar  solid  state  system  for  use  with  cobalt 
teletherapy  units. 

Neither  the  electromechanical  nor  the  electronic  methods  of 
digital  integration  discussed  above  are  applicable  to  depth  dose  scanning 
for  several  reasons.   Electromechanical  methods  are  not  sufficiently 
accurate  because  of  nonlinear  motor  response  to  voltage.   The  best  of 
the  electronic  methods  has  a  repetition  rate  of  just  over  a  hundred 
per  minute  which  is  too  slow  for  measuring  while  scanning.   The  magni- 
tude of  the  ionization  currents  measured  in  the  dose  monitor  devices 
is  of  the  order  of  10  °   amperes,  with  10  ^  amperes  minimum.   These 
currents  are  1,000  to  10,000  times  greater  than  those  encountered  in 
depth  dose  scanning.   The  methods  of  recharging  or  discharging  the 
capacitor  involve  switching  at  the  input  to-  the  electrometer  amplifier. 
Switching  transients  and  leakage  currents  become  quite  important  when 
currents  below  10  '  amperes  are  measured  and  these  problems  have 
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limited  the  presently  used  digital  integration  techniques  to  measure- 
ment of  large  currents.   Furthermore,  the  digital  integration  tech- 
niques discussed  above  are  not  sufficiently  accurate  for  depth  dose 
scanning  since  the  measurement  may  be  in  error  by  as  much  as  one  com- 
plete cycle. 

The  Repetitive  Integral  Measurement  Method 

A  new  method  of  measurement  based  on  repetitive  integrate-reset 
cycles  has  been  developed.   Significant  modification  of  the  techniques 
used  previously  in  dose  monitors  allows  the  new  method  to  be  applied 
to  depth  dose  scanning.   The  method  has  been  named  the  "Repetitive 
Integral  Measurement  Method." 

The  details  of  the  new  method  may  best  be  described  by  reconsider- 
ation of  some  of  the  results  developed  in  Chapter  3.   The  general 
solution  of  the  differential  equation  (6) ,  relating  ionization  current 
to  the  voltage  output  of  the  electrometer  amplifier,  is  given  by 
equation  (8)  as 


V(t)  =  [exp(-t/RC)/C]   I   exp(x/RC)  I(x)  dx 

'0 


•-'0 


when  RC  becomes  large  this  becomes 


X 


V(t)  =  ^   J  I(x)  dx  (11) 


C 

0 


showing  that  the  voltage  output  of  the  electrometer  is  proportional 
to  the  integral  of  the  input  current  from  time  zero  to  time  t.   The 
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average  current  1(0,  t)  in  the  interval  [0,  t]  is 


1(0 


,  t)  =  1   I  I(x)  dx 

t  Jo 


or 


Jo 


tl(0.  t)  =   I  I(x)  dx  (12) 


Substituting  equation  (12)  into  equation  (11) 


v(t)  =  -^iio-ll 


or 


1(0,  t)  =  I  V(t)  (13) 


and  it  is  seen  that  V(t),  the  voltage  present  at  the  electrometer 
output  at  time  t,  is  a  measure  of  the  average  current  over  the  time 
interval  [0,  t]. 

The  integral  measurement  is  a  sum  of  all  values  of  the  input 
current  over  the  interval  of  integration  and  thus  the  value  of  the 
integral  can  give  no  information  concerning  the  variation  of  the  input 
current  at  any  particular  tirae.   As  nothing  can  be  said  concerning  the 
actual  value  of  the  input  current  at  a  specific  time,  it  is  reasonable 
to  estimate  the  value  of  the  input  current  at  the  center  of  the  interval 
of  integration  by  the  average  current  1(0,  t) .   This  may  be  expressed, 
using  equation  (13) ,  as 
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I(t/2)  ~  1(0,  t)  =  ^  v(t) 


(14) 


Equation  (14)  shows  that  the  voltage  V(t)  at  the  electrometer  output 
at  time  t  is  a  measure  of  the  input  current  at  the  center  of  the 
measuring  interval. 

An  input  current  may  be  approximated  over  a  total  time  t  by 
estimating  the  current  value  at  the  center  of  each  of  m  measuring 
intervals  by  the  average  current  in  each  of  the  intervals.   The  input 
current  function  may  be  approximated  by  the  average  current  in  each  of 
m  equal  intervals  as  indicated  by 


/(n  +  1)  At 
I(x)dx._£v[(„+1,  .tl 
nA  t- 


(n  +  1)  At 
I(x 

nAt 

n  =  0,  1,  2,  .  .  .  m  -  1 


where  At  =  t/m  is  the  length  of  each  of  the  m  intervals.   In  the  above 
expression,  the  value  of  the  input  current  at  time  (n  +  1/2)  At  is 
approximated  by  the  average  current  in  the  interval  [nAt,  (n  +  1)  At] 
which  is,  as  shown  by  equation  (14).  proportional  to  the  output  voltage 
of  the  electrometer  at  time  (n  +  1)  At. 

The  voltage  V[ (n  +  1)  At]  in  the  above  equation  is  a  measure  of 
the  average  current  in  the  interval  [nAt,  (n  +  1)  At]  only  if  the 
integrating  capacitor  was  in  a  reset  state  prior  to  start  of  the  inte- 
gration at  time  nAt.   The  value  of  the  integral  over  each  of  the  m 
intervals  must  be  independent  of  the  value  in  any  other  interval. 

The  relationship  between  output  voltage  and  the  average  current 
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in  an  interval  may  be  shown  by  combining  equation  (11)  and  equation 
(12). 


(n  +   1)    At 


/(n  +   1) 
I(x)   dx 
nAt 


V[(n  +   1)    At]    =  i        j      I(x)   dx  =  Ai  I[nAt,    (n  +   1)    At] 


The  output  voltage  at  time  (n  +  1)  At  is  proportional  to  the  average 
current  in  the  interval  [nAt,  (n  +  1)  At]  when  the  voltage  at  time 
nAt  is  zero.   If  the  capacitor  is  not  reset  at  time  (n  +  1)  At  and 
Integration  continues  until  (n  +  2)  At  then 


/ 


(n  +  2)  At 
V[(n  +  2)  At]  =  ^       I   I(x)  dx  =  ^^   I[nAt,  (n  +  2)  At] 

nAt 


and  V[(n  +  2)  At]  is  proportional  to  the  average  current  over  the 
Interval  [nAt,  (n  +  2)  At].   The  necessity  of  resetting  the  capacitor 
prior  to  integration  in  each  interval  is  clear.   Without  resetting, 
the  output  voltage  reflects  the  average  current  from  the  time  inte- 
gration began  and  not  the  average  current  in  each  separate  time  inter- 
val. 

In  summary,  the  repetitive  integral  measurement  method  approxi- 
mates an  input  current  through  repetitive  cycles  of  integrating  and 
resetting.  The  output  voltage  generated  in  each  integration  cycle 
is  proportional  to  the  average  current  over  the  interval  and  is  used 
to  approximate  the  value  of  the  input  current  at  the  center  of  the 
measuring  interval. 

Constant  and  Linear  Current 

When  the  input  current  function  is  constant  (I(x)  =  Iq)  the 
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average  current  in  any  interval  [nAt,  (n  +  1)  At]  is  given  by 


(a  +  1)  At 
I[nAt,  (n  +  1)  At]  =  -^   I   Iq  dx  =  Iq 


1  fV 

•f    nAt 


and  in  this  case  the  current  value  in  the  middle  of  the  measuring 
interval  is  equal  to  the  average  current. 

Likewise  for  a  linear  function  (I(x)  =  kx,  for  k  constant)  the 
average  current  in  any  interval  [nAt,  (n  +  1)  At]  is  given  by 


/(n  +  1) 


(n  +  1)  At 
I[nAt,  (n  +  1)  At]  =  -^  |   kx  dx  =  k(n  +  1/2)  At 


and  again  the  value  of  the  average  current  over  the  measuring  inter- 
val is  equal  to  the  current  at  the  midpoint  of  the  interval. 

These  results  show  that  the  repetitive  integral  measurement 
method  will  give  exact  results  with  a  constant  or  linear  function  of 
current.   This  also  implies  that  if  it  is  possible  to  choose  the  inte- 
gration intervals  so  small  that  the  function  is  essentially  constant 
or  linear  within  the  interval  then  the  repetitive  integral  measurement 
method  will  result  in  very  good  approximations  to  the  actual  values  of 
the  input  function  at  the  center  of  the  measuring  interval. 

Summary 

The  rate  measurement  method  is  limited  in  depth  dose  scanning 
by  the  time  constant  of  the  amplifier.   The  time  constant  is  determined 
by  the  magnitude  of  the  current  to  be  measured  and  is  of  the  order  of 
100  milliseconds  with  the  currents  encountered  in  depth  dose  scanning. 
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With  a  time  constant  of  this  magnitude,  measurement  values  suffer 
significant  signal  lag  in  regions  where  there  is  rapid  change  in  the' 
input  function.   Correction  techniques  must  be  applied  to  measured 
data,  if  scanning  speeds  exceeding  a  few  millimeters  per  second  are 
used. 

A  new  measurement  technique  based  on  integral  measurement  has 
been  developed  for  use  in  depth  dose  scanning.   The  method  is  called 
the  repetitive  integral  measurement  method  and  approximates  an  input 
current  through  repetitive  cycles  of  Integrating  and  resetting  over 
short  time  intervals.   The  output  voltage  at  the  end  of  each  inte- 
gration cycle  is  proportional  to  the  average  current  over  the  interval 
and  is  used  as  an  approximation  of  the  input  current  at  the  center 
of  the  interval  of  integration. 

A  comparison  of  both  measurement  methods  Is  made  in  the  following 
chapter . 


CHAPTER  5 


COMPARISON  OF  THE  RATE  AMD  REPETITIVE  IUTEG8AL 
MEASUREMENT  >!ETHODS 


The  Step  Fimctior.  Regponga 

The  rate  measurement  ..ethod  has  been  shown  to  present  a  problem 
In  those  parts  of  the  radiation  field  where  rapid  changes  occur,  .^en 
the  input  signal  increases  or  decreases  rapidly  the  amplifier  output 
is  late  I„  increasing  or  decreasing  due  to  the  nonzero  time  constant. 

The  unit  step  function  has  the  value  .ero  up  to  a  prescribed 
point  at  Which  the  function  value  becomes  unity.   It  is  an  example  of 
a  function  with  the  most  rapid  change  possible  and  Is  useful  m 
evaluating  the  response  of  an  amplifying  system  to  a  changing  signal. 
The  unit  step  function  has  been  used  to  compare  the  responses  of  the 
-te  measurement  and  repetitive  integral  measurement  methods. 

The  step  function  may  be  defined  as  a  function  of  time  in  the 
following  way 

r 

ift^  =.•  i   t  >  T  ^ 

|0   t  <  T  "^^^""^   ^'  ^s  ^  fi^ed  time 

^e  responses  of  both  methods  of  measurenent  to  the  above  step  function 

are  determined  by  evalua^^•no  t-v,^ 

^y    t:vaj.uating  the  expressions 
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w  n 


R(t)  =  [exp(-t/T)/C]   I  exp(x/T)  j (x)  dx  (15) 

0 

for  the  rate  measurement,  where  t  =  RC,  and 

(n  +  1)  At 
S[(n  +  1/2)  At]  ={exp[-(n  +  1)  At/T]/C}   I   expCx/x)  j (x)  dx      (16) 


n  -  0,  1,  2, 


/(n  H 
exi 
nAt 


for  the  repetitive  integral  method.   The  time  unit  in  the  above 
equations  is  arbitrary.   Equation  (16)  produces  a  value  which  is  pro- 
portional to  the  amplitude  of  the  step.   For  a  meaningful  comparison 
with  the  rate  measurement  this  result  must  be  normalized  so  that 
equation  (16)  will  produce  a  value  equal  to  the  step  amplitude  for 
any  measuring  interval  fully  inside  the  step.  This  normalization  is 
performed  by  dividing  the  result  of  equation  (16)  obtained  in  each 
measuring  interval  [nAt,  (n  +  1)  At]  by  the  result  obtained  in  any 
one  interval  contained  entirely  within  the  step. 

Suppose  the  interval  (kAt,  (k  +  1)  At]  lies  entirely  within  the 
step  (kAt  >  T)  then  the  result  of  equation  (16)  in  this  interval  is 

(k  +  1)  At 
S[(k  +  1/2)  At]  =  {exp[-(k  +  1)  At/T]/C}    I   exp(x/T)  dx 


/(k  H 
exp 
kAt 


=  1[1  -  exp(-At/T)] 


The  normalized  values  are  given  by 
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„,,   ^,/^,  ,,1    S[(n  +  1/2)  At] 
N[(n  +  l/2)  At]  =  s>[(k  +  l/2)  At] 


/ 


(n  +  1)  At 


exp[-(n  +  1)  Lt/j] 

t[1  -  exp(-At/T)]    I   exp(x/T)  j(x)  dx    (17) 


nAt 


n  =  0,  1,  2, 


where  the  interval  [kAt,  (k  +  1)  At]  lies  entirely  within  the  step. 
Equation  (17)  gives  values  which  may  be  compared  directly  with  the 
results  produced  by  equation  (15) .   The  values  resulting  from  evalu- 
ation of  equation  (17)  in  each  of  the  intervals  is  interpreted  as  the 
value  of  the  input  function  j [ (n  +  1/2)  At]  in  the  middle  of  the 
measuring  interval.   Equation  (15)  is  evaluated  at  corresponding 
points  so  that  comparison  is  possible.   It  is  to  be  noted  that 
equations  (16)  and  (17)  contain  exponential  terms  not  present  in 
equation  (11),  the  equation  on  which  the  definition  of  the  repetitive 
Integral  measurement  method  is  based.   The  exponential  terms  are 
included  so  that  the  effect  of  the  time  constant  t  on  the  repetitive 
integral  method  may  be  observed.   When  x  becomes  large  the  exponential 
terms  reduce  to  unity  as  in  other  discussions  of  the  measurement  method. 

Point  Detector  Response 

Figure  5  shows  a  series  of  results  for  each  measurement  method 
for  several  values  of  the  time  constant  and  several  positions  of 
measurement  points  relative  to  the  step  location.   The  figure  repre- 
sents the  response  of  a  point  detector  to  a  time-dependent  step  input 
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Figure  5.   Point  Detector  Response  to  Step  Function  Input. 
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when  measurements  are  made  at  equally  spaced  time  intervals ,   All 
points  are  determined  at  time  intervals  separated  by  0.2  units  and 
straight  lines  connect  each  series  of  points.   For  example,  the 
calculated  values  at  time  0,  0.2,  0.4,  0.6,  0.8,  1.0  and  1.2  correspond 
to  the  same  series  of  measurements  and  are  connected  by  a  straight 
line.   Likewise  the  calculated  values  at  time  0.04,  0.24,  0.44,  0.64, 
0.84  and  1.04  correspond.   It  is  emphasized  that  each  value  determined 
with  the  repetitive  integral  measurement  method  is  the  normalized 
result  of  an  integration  over  a  time  interval.  The  value  is  plotted 
at  the  center  of  the  integration  interval.   In  each  interval  of 
measurement  At  =  0.2.  The  integral  expressions  in  both  equations  (15) 
and  (17)  were  evaluated  by  three-point  quadrature.   Equation  (15)  was 
evaluated  at  points  corresponding  to  the  center  of  the  interval  of 
integration  for  equation  (17) . 

Figures  5a  and  5b  show  the  results  of  the  two  measurement  methods 
with  a  time  constant  x  =  0.01.  From  equation  (8)  it  is  seen  that  the 
basic  integral  expression  evaluated  for  the  two  methods  is  identical 
when  T  is  the  same.   There  is  a  significant  difference  however  in  the 
way  the  integral  expression  is  evaluated  in  each  method.   ^'Jhen  a  rate 
evaluation  is  made  the  integral  in  equation  (8)  is  determined  by 
integration  from  time  zero  to  the  time  of  sampling.   TJhen  a  repetitive 
integral  evaluation  is  made  the  integral  is  determined  by  integration 
from  time  nAt  to  time  (n  +  1)  At  and  the  result  is  normalized  as 
discussed  above.   Figures  5a,  5b,  5c  and  5d  show  that  even  with  the 
same  value  of  t  the  repecitive  integral  method  reaches  the  step  value 
sooner  than  the  rate  method.   The  repetitive  integral  method  overesti- 
mates the  step  in  Figure  5b  and  rises  too  quickly  due  to  integration 
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over  an  interval  At.   The  rate  method  also  overestimates  the  step 
when  measurements  are  made  at  the  same  interval  At  apart  and  the  output 
is  represented  by  a  straight  line  between  points.   This  area  where 
overestimation  of  the  step  occurs  is  directly  dependent  upon  the  time 
interval  between  sampling  with  the  rate  measurement  and  the  time 
interval  over  which  integration  is  performed  in  the  case  of  the 
repetitive  integral  method.  The  repetitive  integral  method  is  clearly 
superior  to  the  rate  method  when  t  =  0.5,  although  again  the  step  is 
overestimated . 

Figure  5e  shows  the  response  of  the  repetitive  integral  method 
to  the  step  input  when  t  >  5.0.   It  was  found  that  Figure  5e  did  not 
change  for  values  of  t  >  5.0,  which  means  that  the  exponential  term 
in  equation  (16)  is  essentially  unity  and  thus  a  simple  integration  of 
the  input  function  is  performed.  The  response  when  t  >  5.0  is  then 
that  which  is  produced  by  the  repetitive  integral  method  as  it  was 
originally  defined,  that  is,   a  measure  of  the  average  of  the  input 
function  over  the  interval  of  measurement.   The  response  to  the  step 
input  of  the  repetitive  integral  method  when  t  <  5.0,  as  shown  in 
Figures  5b  and  5d,  is  not  a  measure  of  the  average  of  the  input  function 
but  an  average  of  the  product  of  the  exponential  function  for  the 
appropriate  x  and  the  input  function  over  the  interval  of  measurement. 
Only  with  x  >  5.0  does  the  repetitive  integral  method  produce  an  output 
which  is  proportional  to  the  average  input  current.   As  the  average 
input  current  is  the  only  reasonable  estimator  of  the  input  within  the 
Interval,  the  repetitive  integral  method  would  in  practice  always  use 
a  X  >  5.0. 

The  repetitive  integral  method,  as  shown  in  Figure  5e,  responds 
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most  quickly  to  the  step  when  the  step  occurs  at  the  beginning  of  the 
integrating  interval  in  which  case  the  step  value  is  reached  within 
one  measuring  interval.   The  worst  case  in  response  is  when  the  step 
occurs  in  the  middle  of  the  integrating  interval.   In  this  case,  two 
intervals  are  required  to  reach  the  step  value.   Results  which  have 
been  obtained  with  values  of  t  not  shown  in  Figure  5  indicate  that 
the  repetitive  integral  method  with  t  >  5.0  always  reaches  the  value 
of  the  step  before  the  rate  method  with  t  >  0.02.   The  repetitive 
integral  method  with  t  >  5.0  reaches  the  step  value  at  least  two 
intervals  earlier  than  the  rate  method  with  t  >  0.05. 

All  the  results  shown  in  Figure  5  and  the  above  discussion  of  the 
figure  pertain  to  measurements  made  with  a  perfect  point  detector, 
which  is  physically  impossible  to  achieve.   The  dimensions  of  the 
detector  become  important  when  evaluating  its  response  to  a  changing 
signal.   Ionization  chambers  are  most  often  cylindrical  because  of 
certain  electrical  characteristics  and  because  of  simplicity  of 
fabrication.   The  response  of  the  chamber  is  proportional  to  the  volume 
irradiated  or  if  a  two-dimensional  cross  section  is  considered,  the 
response  is  proportional  to  the  area  irradiated.   It  is  sufficient  to 
consider  the  two-dimensional  case  since  the  entire  length  of  the 
sensitive  volume  is  equally  irradiated  and  the  cross-sectional  area 
irradiated  is  then  proportional  to  the  volume  irradiated. 

Cylindrical  Detector  Response 

Figure  6  shows  a  cross-sectional  view  of  a  cylindrical  ionization 
chamber  of  radius  R.   The  shaded  section  represents  that  portion  of  the 
chamber  irradiated.   The  dividing  line  between  the  irradiated  and  the 
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unirradiated  portions  forms  a  chord  of  distance  x  from  the  center  of 
the  cylinder.   The  response  of  the  chamber  is  proportional  to  the  area 
irradiated  and  the  fraction  of  full  response  is  the  ratio  of  the  shaded 
area  to  the  total  area.   The  area  of  the  shaded  section  is  given  by  the 
expression^^ 


A(x)  =^-  [xVr'  -  x= 


+  r2  Arcsin  (x/R) ] 


for  X  in  the  interval  [-R,  R].   The  ratio  of  the  shaded  area  to  the 
total  area  is 

_,  .    A(x)         1         I 

^^''■^  -  ~1^   "  °-5  -  7  [(x/R)  yi  -  (x/R)  2  +  Arcsin  (x/R)] 

for  x  in  the  interval  [-R,  R].   The  response  of  the  cylindrical  detector 
is  described  by  F(x)  as  the  chamber  moves  across  an  irradiation  step. 
If  the  step  occurs  at  time  tg  and  the  chamber  moves  at  a  rate  of  k 
distance  units  per  time  unit 


t  <  tg  -  R/k 


t  -  R/k  <  t  <  t^  +  R/k 


t  >  tg  +  R/k 


is  the  expression  of  the  chamber  response  to  a  step  input.   The  function 
g(t)  is  used  in  equations  (15)  and  (17)  to  compare  the  results  of  the 
two  measurement  methods  for  the  cylindrical  detector. 

Figure  7  shows  the  response  of  the  cylindrical  detector  to  a 
step  input.   The  chamber  diameter  in  this  example  was  taken  to  be  the 
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Figure  7.   Cylindrical  Detector  Response  to  Step  Function  Input, 
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distance  transversed  during  one  measuring  interval  or,  in  time  units, 
it  corresponds  to  the  time  between  individual  measurements.   When 
Figures  5a  and  7a  are  compared,  it  is  seen  that  the  cylindrical 
detector  is  delayed  in  reaching  the  step  value  and  in  the  worst  case 
requires  three  intervals  to  reach  the  step  amplitude  whereas  the  point 
detector  requires  only  two  intervals.   The  repetitive  integral  method 
response  shown  in  Figure  7b  reaches  the  step  value  sooner  than  the  rate 
method,  but  overestimates  the  step  more  than  the  rate  method.   With 
the  large  value  of  t  =  0.5  in  Figures  7c  and  7d,  the  response  is  quite 
similar  to  that  with  the  point  detector,  but  both  methods  anticipate 
the  step  before  it  arrives  due  to  the  finite  size  of  the  detector.   The 
response  of  the  repetitive  integral  method  with  i  >  5.0  is  shown  in 
Figure  7e.   In  the  worst  case,  three  intervals  are  required  to  reach 
the  step  amplitude  and  in  the  best  case,  two  intervals  are  needed. 

As  the  cylindrical  chamber  is  made  smaller,  the  response  will 
become  closer  to  that  of  the  point  detector.   If  a  larger  chamber  is 
used,  the  number  of  measuring  intervals  required  to  reach  the  step 
value  increases.   The  dimensions  of  the  cylindrical  ionization  chamber 
used  in  depth  dose  scanning  are  approximately  the  same  as  those  used  in 
the  example  in  Figure  7  and  thus  its  response  is  well  described  by 
Figure  7. 

Figure  7e  shows  that  the  repetitive  integral  method,  with  t  >  5.0, 
is  superior  to  the  rate  method  in  approximating  the  step  input  except 
for  the  rate  method  with  x  =  0.01,  and  it  is  almost  as  good  in  that  case 
also. 

The  results  shown  in  Figures  5  and  7  and  some  similar  results, 
which  are  not  plotted,  with  other  values  of  x  indicate  that  the 
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repetitive  integral  method  with  t  >  5.0,  which  is  really  a  measure  of 
the  average  input,  is  far  superior  to  the  rate  method  with  x  >  0.02  in 
responding  to  a  rapidly  changing  signal. 

The  Size  of  the  Measuring  Interval 

The  closeness  of  approximation  of  the  step  by  both  methods  depends 
on  two  factors,  the  time  interval  between  measurements  and  the  diameter 
of  the  detector.   The  effect  of  reducing  the  measuring  interval  is 
easily  seen  with  the  point  detector.   Figure  8  shows  the  result  of 
reducing  the  measuring  interval  by  ten  times.   Figures  8a  and  8b  show 
the  response  of  both  methods  of  measurement  for  t  =  0.1  for  the  rate 
method  and  x  >  5.0  for  the  repetitive  integral  method.   Figure  8a 
shows  the  best  case  when  the  step  occurs  at  the  beginning  of  the 
interval.   Figure  8b  illustrates  the  worst  case  when  the  step  occurs 
in  the  middle  of  the  interval  of  measurement.   Figures  8c  and  8d  show 
the  same  results  for  x  =  0.01  for  the  rate  method.   Comparisons  of 
Figures  8  and  5  show  the  improvement  in  approximating  the  step.   It 
is  possible  to  approximate  the  step  as  closely  as  desired  with  the 
repetitive  integral  method  and  a  point  detector  by  choosing  a  suf- 
ficiently small  interval  of  measurement. 

In  the  case  of  the  cylindrical  detector,  the  closeness  of  approxi- 
mation of  the  step  is  limited  by  the  detector  diameter.   Reducing  the 
measurement  interval  improves  the  approximation  but  no  matter  how  small 
the  interval  becomes,  the  step  amplitude  is  first  reached  at  a  point 
inside  the  step  by  a  distance  equal  to  the  detector  radius. 
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Figure  8.   Point  Detector  Response  to  Step  Function  Input  with 
a  Tenfold  Reduction  in  the  Size  of  the  Measuring 
Interval. 
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Lateral  Dose  Function  Response 

The  lateral  dose  function  given  by  equation  (9)  represents  the 
input  current  expected  at  the  electrometer  during  a  scan  across  a 
radiation  field.   This  function  provides  a  means  of  evaluating  both 
measurement  methods  as  they  would  be  applied  in  depth  dose  scanning. 
Equations  (15)  and  (17)  were  evaluated,  by  numeric  integration,  using 
the  lateral  dose  function  as  the  input  function.   The  solid  curve  in 
Figure  9  shows  a  plot  of  the  lateral  dose  function.  The  dotted  and 
dashed  curves  are  the  result  of  evaluating  equation  (15)  and  correspond 
to  scans  across  the  field  in  opposite  directions.  The  amplifier  time 
constant  was  taken  as  0.1  seconds,  a  value  typical  in  depth  dose 
scanning,  and  the  scanning  speed  was  2.54  centimeters  per  second.   The 
curves  for  the  scans  in  opposite  directions  show  the  shifts  from  the 
true  value  which  are  characteristic  of  the  rate-measurement  method. 

The  repetitive  integral  measurement  method  was  applied  to  the 
lateral  dose  function  under  the  same  conditions  as  with  the  rate  method 
except  that  the  time  constant  for  the  repetitive  integral  method  was 
chosen  equal  to  50  seconds  so  that  the  electrometer  amplifier  output 
was  proportional  to  the  average  current.   Figure  10  shows  a  plot  of 
the  lateral  dose  function  and  the  results  of  the  evaluation  of  equation 
(17)  for  scans  in  opposite  directions.  The  three  curves  superimposed 
show  no  deviation  greater  than  the  line  width. 

Figure  10  shows  that  the  repetitive  integral  measurement  method 
produces  identical  results  for  scans  in  opposite  directions,  completely 
eliminating  the  shift  introduced  by  the  rate  method.   The  results  shown 
in  Figure  10  were  determined  with  a  scanning  speed  of  2.54  centimeters 
per  second  and  At  =  0.125  seconds. 
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The  accuracy  of  the  repetitive  integral  measurement  method  in 
reproducing  the  lateral  dose  function  was  evaluated  by  consideration 
of  two  types  of  errors,  the  function  error  and  the  displacement  error. 
Figure  11  shows  a  diagram  of  a  function  F,  over  interval  [nAt,  (n  +  1) 
At],  the  integral  approximation  of  the  midpoint  and  the  two  types  of 
errors.   The  repetitive  integral  method  approximates  the  function  value 
at  the  midpoint  of  the  interval  by  the  integral  of  the  function  over 
the  interval.   Only  rarely  does  the  function  value  at  the  midpoint 
equal  the  integral  value  and,  thus,  the  integral  value  lies  above  or 
below  the  function  value.  The  difference  between  the  function  value 
and  the  integral  value  is  denoted  by  AF  and  is  called  the  function 
error.   The  distance  from  the  midpoint  of  the  interval  to  the  point 
on  the  abscissa,  where  the  function  value  and  integral  value  are  equal, 
is  denoted  by  AX  and  is  called  the  displacement  error. 

Both  AF  and  AX  were  determined  for  both  measurement  methods,  as 
applied  to  the  lateral  dose  function.   These  errors  were  determined 
from  the  same  data  used  in  plotting  Figures  9  and  10.   The  time  constant 
was  taken  as  0.1  second  for  the  rate  method,  as  this  corresponds  to 
the  minimum  value  used  in  depth  dose  scanning.   These  errors  were 
determined  for  each  measuring  interval,  as  the  lateral  dose  function 
increased  from  0  to  100  percent.   In  the  case  of  the  repetitive 
integral  method,  the  value  of  the  integral  in  each  interval  [nAt,  (n  +  1) 
At]  was  normalized  to  a  100  percent  maximum  by  dividing  each  integral 
value  by  the  integral  value  obtained  when  the  input  function  remained 
at  its  maximum  during  the  entire  measuring  interval.   This  corresponds 
to  normalization  to  the  value  obtained  at  the  central  axis. 

Table  2  gives  the  values  of  AF  and  AX  for  both  the  repetitive 
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Integral  method  and  rate  method,  as  applied  to  the  lateral  dose 
function.  All  errors  are  the  difference  between  the  actual  value 
and  the  experimental  value.   With  the  rate  method,  the  function  error 
AF  exceeds  21  percent  in  the  region  where  the  function  is  increasing 
rapidly.   The  displacement  error  AX  continually  increases,  a  fact 
observed  in  Figure  9.   Both  the  function  errors  and  the  displacement 
errors  are  large  with  the  rate  method  showing  the  unsuitabllity  of 
the  method  in  this  application.   Both  Che  function  errors  and  dis- 
placement errors  with  the  repetitive  integral  method  are  quite  small 
showing  that  the  method  very  closely  approximates  the  lateral  dose 
function. 

The  Effect  of  Scanning  Speed 

It  has  been  noted,  from  experience  in  depth  dose  scanning,  that 
the  rate  measurement  results  are  directly  dependent  upon  the  speed 
of  scanning.   The  effect  of  scanning  speed  may  be  observed  by  consider- 
ation of  equation  (8),  which  may  be  rewritten  as 


/ 


V(t)  =^  (  exp[(x  -  t)/T]  I(x)  dx  (18) 


where  t  =  RC. 

The  above  equation  gives  the  results  of  measuring  the  current  I(x) 
under  the  influence  of  a  time  constant  t .   If  the  scanning  speed  is 
increased  m  times,  the  current  function  in  equation  (18)  becomes  I(mx) 
and  a  larger  range  of  function  values  is  included  in  the  integral.   If 
the  measurement  of  the  current  is  performed  at  time  t/m,  then 
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•^  0 


t/m 
V(t/m)  =^  I   exp[(x  -  t/m)/T]  I(mx)  dx  (19) 


In  order  to  compare  the  results  of  equations  (18)  and  (19) ,  it  is 
convenient  to  make  a  simple  change  of  variable  in  equation  (19)  by 
letting  y  =  mx,  then 


/: 


V(t/m)  =  -sx      j  exp[(y  -  t)/mT]  I(y)  dy  (20) 


raC 

■  0 


Equation  (20)  shows  that  the  expression  under  the  integral  sign  is 
identical  with  that  in  equation  (18),  except  that  the  time  constant  is 
effectively  increased  by  m  times,   l^^en  the  rate  measurement  method  is 
used,  the  effect  of  increasing  the  scanning  speed  m  times  is  an  increase 
in  the  time  constant  by  m  times,  relative  to  V(t),  as  given  by  equation 
(18). 

With  the  repetitive  integral  measurement  method,  the  value  pro- 
portional to  the  average  current  over  a  single  measuring  interval  is 
given  by 

(n  +  1)  At 
V[(n  +  1)  At]  =  ^      I   I(x)  dx  (21) 


1      r  ^"^ 

c    I    i(- 

J     nAt 


If  the  scanning  speed  is  increased  m  times,  the  current  function 
becomes  I(mx),  as  above.   If  the  measuring  interval  At  is  divided  by  m 
then 

r  (n  +  1)  At 

V[(n  +  1)  ^]  =  ^   I     i(mx)'"dx  (22) 

•'   nAt 

m 
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Again  changing  variables  by  letting  y  =  mx  and  substituting  into 
equation  (22) 


/(n  + 
nAt 


(n  +  1)  At 


V[(n  +  1)  Ak]  =  _^   I      i(y)  dy  (23) 


Comparing  equations  (23)  and  (21),  it  is  seen  that  the  only  effect  of 
increased  scanning  speed  on  the  repetitive  integral  method  is  a 
reduction  of  the  magnitude  of  the  output  voltage  by  1/m.  Theoretically, 
this  is  of  no  consequence,  since  all  the  results  of  the  repetitive 
integral  method  are  normalized  and,  thus,  this  factor  disappears.   It 
is  of  practical  importance,  however,  since  the  magnitude  must  be  large 
enough  to  be  distinguishable  from  noise.   It  may  be  said  that  the 
repetitive  integral  method  is  independent  of  scanning  speed  as  only 
the  signal  magnitude  changes. 

Equations  (18)  through  (23)  apply  equally  well  in  describing 
the  effect  of  reduced  scanning  speed  on  both  measurement  methods. 

Advantages  of  the  Repetitive  Integral  Measurement  Method 

The  previous  results  in  this  chapter  indicate  certain  advantages 
of  the  repetitive  integral  measurement  method,  when  compared  to  the 
rate  measurement  method.   These  advantages,  which  are  summarized  below, 
are  of  particular  importance  in  depth  dose  scanning  and  provide  the 
basis  for  the  development  of  a  totally  new  approach  in  scanning. 

Elimination  of  the  Effect  of  the  Time  Constant 

When  the  rate  method  is  applied  in  depth  dose  scanning,  the  lag  in 
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output  response,  a  result  of  the  nonzero  time  constant,  causes  the 
plot  of  output  voltage  versus  detector  position  to  be  shifted  away 
from  the  correct  response  toward  the  direction  of  the  scan.   The 
magnitude  of  the  shift  is  directly  dependent  upon  the  value  of  the 
time  constant.   Figure  9  shows  results  typical  of  those  obtained  in 
depth  dose  scanning  with  the  smallest  possible  time  constant  employed. 
The  results  of  scans  in  opposite  directions  show  the  lack  of  symmetry 
about  the  central  axis,  which  leads  to  a  zig-zag  effect  when  the 
recorded  data  are  plotted.  Figure  10  and  Table  2  show  that  the 
repetitive  integral  method  applied  to  the  same  input  data  produces  . 
much  smaller  displacement  and  function  errors  and  reproduces  the  plot 
of  the  input  function  so  well  that  no  deviations  are  evident.   Figures 
5  and  7  show  the  superior  response  of  the  repetitive  integral  method 
applied  to  a  step  input. 

Application  to  Pulsed  Radiation  Sources 

Pulsed  radiation  sources  such  as  betatrons  and  linear  accelerators 
pose  a  particular  problem  in  radiation  measurement  because  the  radiation 
output  is  not  constant,  as  with  cobalt,  but  appears  in  short  pulses  of 
microsecond  duration  and,  thus,  determination  of  the  dose  rate  is 
uncertain.   The  most  reasonable  measure  of  the  dose  rate  would  be  the 
average  dose  rate  over  some  interval.   The  repetitive  integral  method 
is  ideally  suited  for  this  type  measurement,  since  the  integral  result 
is  proportional  to  the  average  input  current. 

Elimination  of  Amplifier  Zero  Drift 

This  advantage  of  the  repetitive  integral  method  does  not  result 
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from  the  theoretical  development  of  the  method  but  becomes  evident 
when  the  electrical  analog  of  the  repetitive  integral  procedure  is 
developed.   All  DC  amplifiers  have  a  certain  amount  of  drift  from 
zero,  which  is  included  as  part  of  the  value  measured  by  the  rate 
method.   The  repetitive  integral  method  can  completely  eliminate  all 
long-term  zero  drift  (drift  over  a  period  much  longer  than  the  measuring 
interval).   Amplifier  zero  drift  may  be  eliminated  by  making  two 
measurements  during  integration,  the  first  soon  after  the  start  of 
integration  and  the  second  at  the  end  of  the  measuring  interval,  and 
using  the  difference  of  the  two  measurements  as  the  integral  value  over 
the  interval.   Any  zero  drift,  which  has  occurred  prior  to  the  interval 
of  measurement,  is  eliminated.   Drift  occurring  during  the  interval  of 
measurement  is  not  eliminated  by  this  method. 

Reduction  of  Signal  Noise 

It  is  well  known  that  the  RC  circuit  shown  in  Figure  1  has  a 
certain  filtering  effect  on  the  input  signal.  When  the  value  of  t  is 
large  the  circuit  functions  as  an  integrator  and  produces  a  filtering 
effect  which  severely  attenuates  high-frequency  signals  while  passing 
low-frequency  signals.   The  integrator  forms  what  is  called  a  low-pass 
filter.   The  filtering  effect  attenuates  high-frequency  noise  components 
superimposed  on  the  input  and  reduces  total  signal  noise.   A  further 
discussion  of  noise  is  included  in  the  next  chapter. 

Independence  of  Scanning  Speed 

The  repetitive  integral  method  is  theoretically  independent  of 
scanning  speed,  as  only  the  magnitude  of  the  output  is  reduced.   The 
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practical  realization  of  the  method  by  electronic  means,  however,  is 
limited  by  the  minimum  signal  which  may  be  distinguished  from  noise. 
The  method  is,  then,  practically  limited  by  the  magnitude  of  the  output 
which  is  dependent  upon  speed.   The  rate  method  is  quite  dependent 
upon  speed,  as  the  time  constant  is  effectively  increased.   The  effective 
increase  in  time  constant  is  proportional  to  the  increase  in  speed. 

Summary 

The  repetitive  integral  measurement  method  has  been  shown  to  be 
superior  to  the  rate  method  in  responding  to  a  step  input  and  to- the 
lateral  dose  function,  which  simulates  the  input  current  expected 
during  depth  dose  scanning. 

The  repetitive  integral  method,  when  applied  in  depth  dose  scanning, 
has  the  following  advantages: 

1.  The  effect  of  the  time  constant  is  eliminated. 

2.  It  is  better  suited  for  measurement  of  pulsed 
radiation  sources  because  the  result  produced  is 
proportional  to  the  average  current  over  the 
measuring  interval. 

3.  Amplifier  long-term  zero  drift  may  be  eliminated 

by  taking  the  difference  of  two  measurements  during 
integration. 

4.  The  integrating  circuit  acts  as  a  low-pass  filter, 
thus  attenuating  high-frequency  noise  components. 

5.  It  is  theoretically  independent  of  scanning  speed 
as  only  the  magnitude  of  the  signal  is  reduced. 

The  actual  implementation  of  the  repetitive  integral  method  by 

electronic  means  is  presented  in  the  following  chapter. 


CHAPTER  6 
THE  REPETITIVE  INTEGRAL  ELECTROiMETER  AMPLIFIER 

The  magnitudes  of  the  currents  produced  in  small  ionization 
chambers  such  as  those  used  in  depth  dose  scanning  are  of  the  order 
of  10    to  10  ^   amperes,  and  thus  sensitive  amplifying  devices  are 
required  for  their  measurement.   Amplifiers  used  for  current  measure- 
ment in  this  range  are  usually  called  electrometer  amplifiers  or  simply 
electrometers.   The  voltage  output  of  the  electronic  electrometer 
amplifier  is  usually  proportional  to  the  input  current. 

Frequently,  operational  amplifiers  are  used  in  the  construction 
of  electrometers.   An  operational  amplifier  is  simply  a  high-gain, 
direct-coupled  amplifier.   The  operational  amplifier  often  has  a 
differential  input  and  is  used  with  external  feedback  networks.   The 
ideal  operational  amplifier  is  represented  by  the  symbol 


•  e. 


where  e-j^  and  69  are  the  input  voltages,  eg  is  the  output  voltage  and  K 
is  the  gain  of  the  amplifier.   The  idealized  amplifier  properties 
usually  assumed  are 
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Gain  (K)  =  » 

eQ  =  0,  when  e^  =  ej 

Input  impedance  =  =° 

Output  impedance  =  0 

Bandwidth  =  <= 

No  temperature  drift. 

These  properties  can  never  be  realized  in  practice  but  the  assumption 
of  idealness  is  useful  in  the  analysis  of  feedback  circuits. 

The  Integration  Circuit 

The  electrical  analog  of  the  mathematical  operation  of  integration 
can  be  easily  obtained  with  the  use  of  the  inverting  input  of  an 
operational  amplifier  with  the  addition  of  a  feedback  capacitor.   The 
integration  circuit  is  diagrammed  below 


1  •-vVV^^^ 


Since  the  input  impedance  is  assumed  infinite,  no  current  flows  into 
the  amplifier.  Applying  Kirchhoff's  law  at  the  junction  at  the  amplifier 
input. 


ei  -  ej  ^  d(eo  -  ej)  ^^^ 


R 


dt 


87 


The  output- voltage  may  be  expressed  as  eg  =  -Ke]_.   Hence,  e^^ ~. 

As  K  in  the  ideal  case  is  infinite,  e^  =  0  and 


^  =  -11 
dt      RC 


1  f' 


^'^  "  "  RC   I  ^1  ^t  +  eo(t  =  0)  (24) 

0 

where  egCt  =  0)  is  the  voltage  across  the  capacitor  at  the  start  of 
integration.   The  output  voltage  is  seen  to  be  proportional  to  the 
integral  of  the  input.   In  this  discussion,  the  noninverting  amplifier 
input  is  connected  to  voltage  common  so  the  amplifier  acts  as  a  single 
input  device.  This  technique  of  Integration  has  been  used  for  many 
years  with  analog  computers  for  the  solution  of  differential  equations. 

The  repetitive  integral  measurement  method  described  in  Chapters 
4  and  5  approximates  the  input  function  by  a  series  of  integrations 
performed  over  short  time  intervals.   Each  of  the  individual  inte- 
grations is  independent  of  any  other  integration,  which  means  the 
feedback  capacitor  is  discharged  to  zero  voltage  before  integration 
begins  in  each  interval.   In  the  case  of  depth  dose  scanning,  the 
ionization  current  is  approximated  at  n  points  across  the  radiation 
field  by  repetitively  integrating  in  n  time  intervals  whose  length  is 
is  At  with  the  integrating  capacitor  being  reset  to  zero  before  inte- 
gration begins  in  each  interval.   An  operational  amplifier  with  feedback 
capacity  is  utilized  to  perform  the  integration  with  the  repetitive 
integral  electrometer. 
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The  Switching  Electrometer 

The  repetitive  integral  measurement  method  demands  that  the 
integrating  capacitor  be  reset  to  zero  and  held  at  zero  until  inte- 
gration of  the  input  is  to  begin,  at  which  time  the  capacitor  must 
be  allowed  to  charge.   Therefore,  the  electrometer  must  be  switched 
at  the  appropriate  time  in  order  that  the  amplifier  output  and  input 
are  shorted  and  the  integrating  capacitor  discharged.   This  procedure 
is  simple  to  describe  but  poses  quite  a  problem  to  accomplish  with  an 
electrometer.   Switching  at  the  input  of  such  a  sensitive  electrometer 
is  quite  difficult  since  transients  accompanying  both  making  and  - 
breaking  of  the  switch  may  send  the  operational  amplifier  into  immediate 
saturation.   Leakage  currents  through  the  switch  during  integration 
are  also  a  problem  with  a  sensitive  electrometer. 

The  methods^°~^^  of  digital  integration  discussed  in  Chapter  4 
employed  switching  at  the  input  of  the  electrometer,  but  the  currents 
involved  were  1,000  to  10,000  times  larger  than  those  encountered  in 
depth  dose  scanning.   Switching  at  the  electrometer  input  is  easily 
accomplished  with  these  large  currents  since  leakage  currents  through 
the  switching  apparatus,  even  though  present,  cause  no  problems  due  to 
the  low  sensitivity  of  the  amplifier. 

The  Switching  Method 

The  connection  of  an  electrometer  output  to  input  for  the  purpose 
of  resetting  may  be  accomplished  by  either  mechanical  or  electronic 
means.   Mechanical  resetting  using  relays  was  tried  and  immediately 
abandoned  because  the  switching  noise  drove  the  amplifier  into  saturation. 
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A  method  of  resetting  bv  electronic  switching  was  successful  and  has 
been  used  in  the  construction  of  the  repetitive  integral  electrometer. 
The  reset  network  is  illustrated  in  Figure  12. 

The  switching  of  the  integrator  to  the  reset  mode  is  accomplished 
with  a  pair  of  p-channel  enhancement-type  metal  oxide  silicon  field- 
effect  transistors  (MOSFET) .   These  transistors,  indicated  as  A  and  B 
in  Figure  12,  are  normally  nonconducting  with  zero  gate  voltage  and 
become  conducting  when  a  negative  voltage  is  applied  at  the  gate.   The 
integrating  capacitor  is  reset  during  the  conducting  phase  when  the 
electrometer  output  is  connected  to  the  input  through  the  transistors. 

The  first  requirement  of  such  a  switching  arrangement  is  that  the 
leakage  through  the  switch  be  so  small  as  to  produce  a  negligible 
effect  on  the  integration.   It  was  shown  in  a  previous  chapter  that 
an  ionization  current  of  10"' '^  amperes,  which  is  typical  in  depth  dose 
scanning,  will  produce  10  volts  at  the  output  of  the  electrometer 
when  a  feedback  resistor  of  10^^  ohms  is  used.   A  leakage  of  just  1 
picoampere  into  the  input  will  result  in  an  output  of  100  millivolts 
which  is  1  percent  of  the  value  measured. 

Two  transistors  are  used  together  as  shown  to  reduce  leakage 
problems.   The  drain  of  one  transistor  is  connected  to  the  source  of 
the  other  and  this  junction  goes  to  voltage  common  over  a  small 
resistance  keeping  the  junction  point  of  the  two  transistors  essentially 
at  zero  volts  during  integration.   The  source  of  transistor  A  in 
Figure  12  goes  to  the  input  of  electrometer  F.   The  input  voltage  to 
F  never  exceeds  a  fraction  of  a  millivolt  so  the  voltage  across 
transistor  A  is  very  small  during  integration.   Furthermore  during 
the  integration  cycle  the  gate  voltage  on  both  transistors  is  zero. 
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thus  the  potential  for  leakage  currents  is  greatly  reduced  due  to  the 
small  voltage  differences.   The  source  of  transistor  B  is,  like  the 
drain  of  A,  essentially  at  zero  volts  but  the  voltage  at  the  drain  of 
transistor  B  increases  to  a  maximum  of  ten  volts  during  integration. 
Any  leakage  current  from  the  drain  of  transistor  B  to  its  source  goes 
to  voltage  common  through  the  resistor  and  does  not  effect  the  input. 
Any  leakage  current  would  not  effect  the  amplifier  output.   The  "cut-off" 
resistance  between  drain  and  source  on  these  transistors  is  in  the 
order  of  10^  ohms  so  the  maximum  leakage  across  transistor  B  is  10  ^ 
amperes  which  has  an  insignificant  effect  on  the  output  of  electrome- 
ter F  but  which,  if  present  at  the  input,  would  constitute  a  signal 
100  times  larger  than  the  maximum  current  encountered  during  depth 
dose  scanning. 

The  above  discussion  shows  the  necessity  for  two  transistors 
connected  as  shown  in  Figure  12.   If  just  one  transistor  were  used 
the  full  10  volts  developed  at  the  output  would  be  the  voltage  across 
the  source  and  drain  and  the  potential  for  leakage  would  be  great. 
Another,  possibly  more  serious,  problem  of  using  just  one  transistor  is 
the  capacity  between  the  source  and  drain.  The  gate  voltage  is  zero 
during  integration  and  thus  this  capacity  is  also  charging  during  inte- 
gration.  This  capacity  is  in  parallel  with  the  integrating  capacity 
C  and  serves  as  a  second  integrating  capacity.   The  problem  is  the  fact 
that  this  capacity  between  source  and  drain  varies  with  the  voltage 
and  thus  the  charge  transferred  to  the  input  varies  during  the  inte- 
gration, producing  a  result  which  is  not  proportional  to  the  input. 
The  capacity  between  drain  and  source  for  this  type  of  transistor  is 
0.5  picofarad.   Since  the  effective  integrating  capacitor  C  is  about 
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this  value  the  varying  capacity  of  the  transistor  is  quite  important 
in  this  circuit  particularly  at  higher  gain. 

The  two  field-effect  transistors  resetting  network,  is  theoreti- 
cally all  that  is  needed  to  accomplish  the  reset  of  the  integrating 
capacity.   A  problem,  however,  does  exist  when  the  transistors  turn 
off.   The  transistors  cease  to  conduct  at  some  negative  voltage  below 
zero  volts  but  the  reset  pulse  at  the  gate  goes  entirely  to  zero  volts. 
The  capacity  between  the  source  and  gate  of  transistor  A  becomes 
important  after  the  transistor  ceases  to  conduct,  when  the  integrating 
capacity  begins  to  charge  and  integration  begins.   The  capacity  between 
the  source  and  gate  of  transistor  A  causes  charges  to  be  transferred 
to  the  input  during  the  time  after  the  transistor  ceases  to  conduct 
and  before  the  gate  voltage  reaches  zero.   This  transfer  of  charges 
causes  an  immediate  change  in  the  output  of  the  amplifier.   This 
problem  may  be  overcome  by  inserting  a  compensating  pulse  at  the  input 
of  the  electrometer.   The  same  reset  pulse  which  causes  the  transistors 
to  become  conducting  is  inverted  in  the  level  buffer  E  and  is  sent 
over  a  potentiometer  and  capacity  to  the  input.  The  output  of  the 
inverter  E  through  the  capacity  D  produces  the  opposite  level  compen- 
sating pulse  to  that  resulting  from  the  capacity  of  transistor  A. 
As  one  pulse  goes  from  low  to  high  the  other  goes  from  high  to  low. 
The  potentiometer  provides  adjustment  of  the  compensating  pulse  ampli- 
tude to  ensure  correct  compensation  at  different  gain  settings  of  the 
amplifier. 

Reset  Network  Analysis 

An  analysis  of  the  reset  network  is  made  in  order  to  obseirve  the 
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functioning  of  the  network  during  the  reset  and  the  integration  phases. 
Figure  13  shows  an  equivalent  circuit  of  the  resetting  network  with 
the  transistors  replaced  by  resistors.   The  constant  input  current  to 
the  electrometer  is  indicated  by  I  and  the  output  voltage  by  Vq.   The 
voltage  at  the  junction  between  the  two  resistors  is  indicated  by  V  . 
Using  Kirchhoff 'slaw  at  the  junction  of  the  reset  resistors  v'  is  found 


t  _  1 

°  1  +  R3/R1  +  R3/R2  (25) 


Applying  Kirchhoff's law  to  the  junction  at  the  electrometer  input. 


I  =  _  dVoc  _  v;. 

dt     R_L  (26) 


Substituting  equation  (25)  into  equation  (26)  and  rearranging. 


1   +  Vq i =  -  I  (27) 


dt      [Ri   +  R3  +  R1R3/R2IC     C 


Integrating  equation  (27), 


Vq  =  -  Il[  1  -  exp(-t/T)]  (28) 


where  t  =  [R^  +  R3  +  RiR3/R2]C 

During  the  time  of  reset,  transistors  A  and  B  in  Figure  12  become 
conducting  and  in  the  conducting  phase  each  has  a  resistance  of  about 
200  ohms.   Hence,  Rj^  =  R2  =  R3  =  200  ohms  during  reset.   Then, 
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T  =  1.8  X  10"  seconds  for  C  =  30  picofarads.   With  such  a  small 
value  of  T  the  output  Vq  drops  immediately  to  zero,  discharging  the 
capacitor  C.   The  capacitor  remains  discharged  as  long  as  the  tran- 
sistors remain  conducting,  that  is,  as  long  as  R^  =  R2  =  R3  =  200  ohms. 

During  the  integration  period  the  transistors  A  and  B  are  not 
conducting  and  each  has  a  resistance  of  about  10^  ohms.   Then  the 
value  of  R]^  =  R^  =  10*^  ohms  and  R2  =  200  ohms.   In  this  case,  t 
=  1.5  X  10^  seconds  for  C  =  30  picofarads  and  v'  =  Vq/(5  x  10^). 
During  integration  the  maximum  voltage  at  the  drain  of  transistor  A 
is  -2  X  10   volts  and  thus  the  voltage  difference  between  the  source 
and  drain  is  very  small.   If  the  exponential  term  in  equation  (28)  is 
expanded  as  a  Taylor  power  series. 


V  =  -  Il[i  -  1  +  1  -    tl_   +   t^  -  .  .  .)  (29) 

u      C        T   ^2   2!   ji    3! 


For  T  large,  equation  (29)  reduces  to 


V  =  -  It 


which  is  the  result  obtained  in  equation  (5)  by  integrating  the  constant 
input  current  I. 

It  is  clear  that  the  resetting  network  acts  to  immediately  reset 
and  hold  the  integrating  capacitor  in  reset  during  the  time  the  tran- 
sistors are  conducting.   The  resetting  network,  furthermore,  does  not 
disturb  the  integrating  characteristics  of  the  amplifier  during  its 
integration  phase. 

The  resetting  network  has  been  used  in  two  repetitive  integral 
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electrometers  and  has  provided  very  reliable  and  stable  resetting 
characteristics.  The  only  disadvantage  of  this  method  of  resetting 
is  the  necessity  for  the  adjustment  of  the  compensating  pulse.   Al- 
though this  adjustment  is  easily  made,  it  does  complicate  the  use 
of  the  electrometer  since  the  compensating  pulse  must  be  readjusted 
whenever  the  electrometer  gain  is  changed.   Once  the  adjustment  is  made, 
however,  it  need  not  be  repeated  as  long  as  the  same  gain  is  used. 

The  Repetitive  Integral  Amplifier  Circuit 

Figure  lA  shows  the  complete  repetitive  integral  amplifier  circuit 
which  is  implemented  in  the  two  amplifiers  which  have  been  fabricated. 
Several  aspects  of  this  circuit  bear  discussion. 

The  Operational  Amplifier 

The  operational  amplifier  is  the  ICH8500A*  and  was  chosen  because 
of  its  extremely  low  input  bias  current  of  less  than  0.01  picoamperes 
and  its  operating  temperature  range  of  -25  degrees  to  +85  degrees  Centi- 
grade.  The  amplifier  is  of  the  field-effect  transistor  input  type, 
guaranteeing  a  high  input  impedance.   The  output  voltage  range  is  ±12 
volts  and  the  large  signal  open  loop  voltage  gain  is  typically  10^.   The 
amplifier  shunt  capacity  is  0.1  picofarad.   The  slew  rate  is  0.5  volts 
per  microsecond.   The  amplifier  characteristics  make  it  particularly 
suitable  for  electrometer  design  in  picoampere  and  subpicoampere  current 
measuring  circuits.   The  amplifier  is  packaged  in  a  standard  TO-5  package, 


*Manufactured  by  Intersil,  Inc.,  Cupertino,  California. 
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The  Resetting  Field-Effect  Transistors 

The  function  and  design  of  the  resetting  network  has  been  previ- 
ously  discussed.   The  switching  components  used  in  the  resetting 
network  are  M108*  p-channel  enhancement-type  metal  oxide  silicon 
field-effect  transistors.   These  transistors  are  the  dual  type  and 
are  conveniently  packaged  together  in  a  TO-99  package.   The  pertinent 
characteristics  of  these  transistors  have  been  previously  discussed. 

The  Compensating  Pulse  Level  Buffer 

The  level  buffer,  or  inverter  as  it  is  sometimes  called,  for  the 
reset  network  compensating  pulse  is  a  CD4007AE**  buffer.  This  is  not 
an  inverter  from  the  standpoint  of  voltage  inversion  from  positive 
to  negative  or  vice  versa,  but  is  a  level  inverter.   The  buffer  output 
is  the  high  level  for  the  circuit  when  the  low  level  is  its  input  and 
vice  versa.   This  is  a  digital  device;  that  is,  it  operates  on  one  of 
two  levels,  low  or  high.   The  function  of  the  buffer  in  this  circuit 
is  to  invert  the  level  of  the  reset  pulse  before  it  reaches  the  input. 

The  Inverting  Operational  Amplifier 

A  second  operational  amplifier,  indicated  by  the  number  741  in 
Figure  14,  is  used  in  the  electrometer  circuit.   This  amplifier  is  a 
yA741,***  commonly  used  in  many  types  of  circuits.   Its  characteristics 
are  not  particularly  important  to  the  circuit  in  Figure  14  since  it 
is  used  only  as  a  voltage  inverter  and  does  not  handle  small  input 


♦Manufactured  by  Siliconix,  Inc.,  Santa  Clara,  California. 
♦♦Manufactured  by  RCA,  Inc.,  Somerville,  New  Jersey. 
♦♦♦Manufactured  by  Fairchild  Semiconductor,  Mountain  View,  California. 
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currents  as  does  the  electrometer  amplifier.   The  input  to  the 
electrometer  amplifier  (ICH8500A)  is  positive  and  thus  the  output  is 
negative.   The  analog-to-digital  converter,  which  is  used  at  the 
output  of  the  repetitive  integral  electrometer  in  Figure  14  accepts 
only  positive  input,  and  so  the  operational  amplifier  (741)  is  used 
to  convert  the  voltage  output  of  the  electrometer  amplifier  to  a 
positive  voltage.   The  amplifier  is  used  in  the  standard  summer 
configuration  with  a  feedback  resistor.   Since  the  gain  in  this 
application  is  unity  only  a  sign  change  is  evident  at  the  output. 

The  Input  Protection  Diodes 

The  two  reversed  diodes  in  Figure  14  at  the  input  of  the  electrome- 
ter are  for  protection  of  the  electrometer  from  excessive  input  voltage 
which  could  result  if  the  ionization  chamber  electrode  became  shorted 
to  the  high  voltage  at  the  chamber  cap.   The  diodes  are  the  low 
leakage  JT2000*  type  and  do  not  disturb  the  input.   If  high  voltage 
is  present  at  the  input,  the  diodes  begin  conducting  at  some  positive 
voltage  well  within  the  range  tolerated  by  the  electrometer  input  and 
thus  protect  the  input.   The  voltage  difference  at  the  electrometer 
inputs  may  be  ±100  volts  without  damage. 

The  Voltage-Divider  Gain  Adjustment 

The  gain  of  an  operational  amplifier  with  feedback  loops  is 
governed  by  the  characteristics  of  the  feedback  element.   This  is 
evident  in  equation  (24) ,  where  the  value  of  the  feedback  capacity 
determines  the  proportionality  constant.   If  a  fixed  value  of  the 


*Manufactured  by  Fairchild  Semiconductor, 
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feedback  capacity  is  used,  the  gain  is  fixed.   It  is  very  desirable  to 
have  a  means  of  gain  adjustment  so  that  the  maximum  possible  output 
signal  may  be  produced  for  a  variety  of  radiation  intensities.   Gain 
adjustment  may  be  easily  obtained  by  means  of  the  voltage-divider 
network  shown  in  Figure  14  at  the  electrometer  output.   In  order  that 
the  effective  feedback  capacity  may  be  varied,  the  feedback  capacity 
is  connected  from  the  electrometer  input  through  the  divider  network 
to  the  electrometer  output.   Figure  14,  for  illustrative  purposes, 
shows  only  four  switch  positions,  but  in  reality  ten  positions  are 
used  providing  ten  different  gain  settings. 

Figure  15  shows  a  diagram  of  the  electrometer  and  a  voltage- 
divider  network  connected  to  the  feedback  capacitor  C.   Since  the 
shunt  capacitance  Cg  inherent  in  the  amplifier  is  important  in  gain 
calculations,  it  is  indicated  in  parallel  with  the  feedback  capacitor. 
Considering  the  currents  flowing  into  the  input. 


I  =  -  ^C_  -  ^C  (30) 

dt      dt 


Now 


R2      V0R2 
^  =  ^O^^^TTJ  =  ^'  "h^^^  Rt  =  Rl  +  R2  (31) 


Substituting  equation  (31)  into  equation  (30)  and  integrating 


Vn  =  -  ^ .    I   r.    ,^         I  I  dx  (32) 

J    ( 


Cg     +     C     R2/RT       J  Q 
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The  gain  if  defined  as  unity  when  R2  =  Rj.,and  thus  the  ratio  of  the 
coefficient  of  the  integral  in  equation  (32)  to  the  coefficient  for 
unity  gain  is  the  amplifier  gain  for  different  positions  on  the 
divider  network.   If  the  gain  is  designated  by  A,  then 


Cg  +  C  R2/R^     Cs  +C 
A  =  ~ 


C3  +  C  R2/RJ, 


^s  ■•"  C  (33) 

_  Cs/C  +  1     ^   R  +  1 
Cg/C  +  R2/RT   R  +  R2/RX 


where  R  =  Cg/C 


Then  A  may  be  expressed  in  terms  of  R^,  the  resistance  between  the 
output  and  capacitor,  by  substituting  R2  =  R^  -  R^^  into  equation  (33) 


Then 


R  +  1 
A  =  


Rt       Rj  [i  +  rJ 


(34) 


Equation  (34)  may  be  solved  for  R-, 

R-L  =  Rt(R  +1)  (1  -  1/A)  (35) 

Equation  (35)  enables  R^  to  be  chosen  for  a  given  gain  A  with  a  total 
divider  network  resistance  R^. 
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Table  3  gives  the  gain  characteristics  for  ten  switch  positions, 
as  used  in  the  first  repetitive  integral  electrometer  amplifier 
constructed.   A  fixed  feedback  capacity  of  30  picofarads  was  used 
with  a  voltage  divider  network  total  resistance  of  21.76  kilohms. 
The  divider  network  provides  the  same  results  which  would  be  obtained 
by  changing  the  value  of  the  feedback  capacity.   An  effective  feedback 
capacity  Cg  may  be  determined  from 


Ce  =  Cs  +  C(l-^)  (36) 


Rt' 


The  effective  feedback  capacity  for  each  switch  position  is  shown 

in  Table  3.   The  value  of  the  shunt  capacitance  Cg  is  given  in  the 

operational  amplifier  (ICH8500A)  specifications  as  0.1  picofarad, 

but  after  the  complete  repetitive  integral  electrometer  amplifier 

was  wired,  the  value  of  Cg  determined  experimentally  was  0.26  picofarad. 

This  value  was  used  in  the  calculations  for  Table  3  since  the  shunt 

capacitance  is  important  in  calculating  larger  gain  values. 

After  the  first  repetitive  integral  electrometer  was  constructed 
and  used  in  depth  dose  scanning,  it  was  found  that  a  ten-position  gain 
switch  was  somewhat  inconvenient  because  no  setting  between  fixed 
switch  positions  was  possible.  If  the  gain  at  one  switch  position  was 
too  high,  the  next  lower  setting  was  used,  but  often  providing  a  less 
than  optimum  gain.   To  eliminate  this  problem  in  the  second  amplifier 
construction,  the  voltage-divider  network  was  replaced  by  a  ten-turn 
potentiometer.   This  provided  a  continuous  gain  selection  from  the 
minimum  to  the  maximum  and  was  much  more  satisfactory  than  fixed  gain 
switches.   The  analysis  of  the  divider  network  applies  equally  well  to 


104 


TABLE  3 
GAIN  CHARACTERISTICS  OF  THE  REPETITIVE  ELECTROMETER  A.MPLIFIER 


Effective 
Feedback 
Capacity 
(pF) 


Gain  Switch 

Ri 

Added 

Resistance 

Position 

(vn) 

Each 

Position 

1 

0 

0 

1.0 

30.26 

2 

8.2 

8.2 

1.6 

18.95 

3 

13.1 

4.9 

2.5 

12.20 

4 

16.0 

2.9 

3.7 

8.20 

5 

16.9 

0.9 

4.3 

6.96 

6 

17.8 

0.9 

5.3 

5.72 

7 

18.7 

0.9 

6.8 

4.48 

8 

19.6 

0.9 

9.3 

3.24 

9 

20.5 

0.9 

15.2 

2.00 

10 

21. A 

0.9 

0.36 

40.0 

.76 

21.76 


R^  =  21.76  kn;     C  =  30  pF;     C   =  0.26  pF 
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the  potentiometer.   An  additional  resistance  of  10^  ohms  was  added  in 
series  with  the  potentiometer  so  that  the  gain  values  shown  by  the 
available  gain  selection  in  Figure  16  are  obtained  when  the  potenti- 
ometer is  varied.   Figure  16  shows  the  graph  of  gain  versus  potenti- 
ometer setting  obtained  with  Rj.  =  2  x  10^  ohms,  C  =  10  picofarads  and 
Cg  =  0.26  picofarads.   The  total  potentiometer  resistance  was  10^  ohms. 

Amplifier  Construction 

Due  to  the  small  currents  originating  in  the  ionization  chamber, 
it  is  desirable  to  locate  the  electrometer  amplifier  as  close  as 
possible  to  the  ionization  chamber.   Some  separation  is,  however, 
unavoidable  because  the  ionization  chamber  must  be  submerged  in  water 
at  the  center  of  the  radiation  field.   The  amplifier  must  be  located 
far  enough  from  the  ionization  chamber  so  that  it  is  never  exposed 
to  the  direct  radiation  beam,  even  with  the  largest  field  size  to 
be  measured.   Figure  17  shows  the  ionization  chamber  and  amplifier 
housing.  The  amplifier  is  located  at  the  end  of  the  stylus  in  the 
housing,  which  remains  outside  the  water  and  outside  the  direct 
radiation  beam.   All  the  required  operating  voltages  and  the  reset 
pulses  are  carried  by  a  cable  to  the  plug  on  the  side  of  the  housing. 
The  amplifier  output  is  also  carried  by  the  cable  to  the  analog-to- 
digital  converter. 

The  repetitive  integral  electrometer  amplifier  is  constructed 
on  a  printed  circuit  board  and  is  mounted  in  the  housing  on  spacers 
so  that  no  component  touches  the  housing.   The  amplifier  housing  is 
not  located  in  the  primary  radiation  beam,  but  does  receive  considerable 
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scattered  radiation.   The  space  inside  the  housing  surrounding  the 
printed  circuit  board  must  be  filled  with  some  material  to  displace 
the  air,  since,  if  it  remains,  ions  produced  in  the  air  surrounding 
the  electrometer  will  be  collected.   The  resulting  current  would  be 
indistinguishable  from  that  produced  at  the  ionization  chamber. 
Ceresin  wax  is  used  to  fill  the  spaces  around  the  amplifier.   Ceresin 
wax,  refined  from  the  mineral  ozokerite,  has  a  melting  point  around 
75  degrees  Centigrade  and  a  volume  resistivity  greater  than  5  x  10^° 
ohm-centimeters,  and  hence  is  well  suited  for  this  application.   The 
liquid  wax  easily  displaces  all  air  inside  the  housing  and  quickly 
solidifies.   However,  the  wax  may  be  easily  removed  from  the  housing 
by  application  of  heat,  with  no  damage  to  the  components.   The  wax  also 
eliminates  most  vibrational  noise  called  microphonics. 

Although  the  major  source  of  extraneous  ionization  is  eliminated 
by  filling  the  housing  with  ceresin  wax,  two  other  small  sources  of 
ionization  exist.   Both  the  electrometer  amplifier  (ICH8500A)  and 
the  field-effect  transistors  are  packaged  in  standard  capped  packages. 
A  large  portion  of  the  space  inside  the  cap  on  each  package  is  filled 
with  some  gas  at  the  time  of  manufacture.   These  gas-filled  spaces  are 
possible  sources  of  ionization  which  would  be  registered  at  the  input 
of  the  electrometer  and  may  be  eliminated  by  decapping  the  packages 
before  installing  the  amplifier  in  the  housing.   Tools  for  decapping 
these  packages  are  commercially  available. 


Amplifier  Control  Logic 


The  repetitive  integral  measurement  method,  as  described  in 
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Chapter  4,  approximates  the  input  current  from  the  ionization  chamber 
at  discrete  points  by  repetitively  integrating  over  short  time  inter- 
vals,  resetting  the  integrator  and  beginning  a  new  integration  cycle. 
The  voltage  developed  during  each  integration  cycle  is  related  to  the 
center  of  the  interval  over  which  the  ionization  chamber  moves  during 
the  time  of  integration.   It  is  clear  that  the  cycles  of  integration 
and  movement  of  the  ionization  chamber  must  be  closely  correlated 
if  reliable  measurements  are  to  be  obtained.   The  correlation  of 
movement  and  measurement  cycles  must  be  reproducible  on  each  scan  line. 
Such  close  correlation  of  these  two  activities  is  easily  accomplished 
with  digital  logic. 

The  Digital  Clock 

A  highly  accurate  crystal-controlled  oscillator  is  used  as  a 
system  clock.   The  clock  frequency  is  16.384  kilohertz  and  is  accurate 
to  0.001  percent.   The  clock  runs  continuously  and  supplies  negative 
pulses  to  both  the  scanning  control  logic  (to  be  discussed  in  a  later 
chapter)  and  the  amplifier  control  logic.   The  clock  frequency  is 
reduced  to  256  hertz  before  entering  the  logic  sections.  The  accuracy 
of  the  clock  ensures  the  proper  correlation  between  movement  and 
measurement  cycles  from  one  scan  line  to  another.   The  clock  is 
sufficiently  accurate  so  that  the  time  required  for  a  particular  number 
of  pulses  to  occur  is  always  identical.   The  clock  pulses  are  used  to 
govern  the  reset  and  integration  cycles  of  the  amplifier  as  well  as 
the  movement  of  the  ionization  chamber. 
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The  Switching  Logic 

Figure  18  shows  a  diagram  of  the  amplifier  switching  logic.   The 
counters  F  and  G,  the  level  inverters  E,  H,  I  and  P,  and  gates  J,  K, 
and  L  are  all  of  the  RCA  COSMOS  4000  family*  of  digital  integrated 
circuits.   Operational  amplifiers  M  and  N  are  the  previously  discussed 
Fairchild  yA7Al  type. 

A  single  repetitive  integral  amplifier  measuring  cycle,  which 
includes  both  the  integration  and  resetting  phases  consists  of  32  clock 
pulses.   The  clock  pulses  at  the  amplifier  have  a  frequency  of  256 
hertz  and  thus  one  measuring  cycle  requires  0.125  seconds.   The  entire 
measuring  cycle  is  repeated  an  integral  number  of  times  while  the 
ionization  chamber  is  moved  across  the  radiation  field. 

Before  scanning  begins,  the  reset  inputs  R  of  the  7-stage  binary 
counters  F  and  G  (see  Figure  18)  are  in  the  high  (0  volts)  state  and 
the  counter  outputs  are  all  low  (-10  volts).   This  causes  the  output 
of  NOR  gate  L  to  be  high,  but  inverter  P  converts  this  to  the  low 
state,  and  the  repetitive  integral  electrometer  remains  reset  until 
the  A  output  of  counter  G  goes  high  at  the  second  clock  pulse.   The 
moment  the  reset  is  removed,  the  amplifier  begins  to  integrate  and 
continues  until  all  the  outputs  A,  B,  C  and  D  of  counter  G  go  low 
on  pulse  32.   The  outputs  of  counter  G  also  control  the  measurement 
timing  of  the  integrated  signal.   Counter  F  serves  only  as  a  pulse 
divider,  so  that  counter  G  counts  only  at  every  second  clock  pulse. 

Figure  19  shows  a  pulse  diagram  for  a  single  measuring  cycle,  which 
begins  with  the  amplifier  in  the  reset  condition.   The  reset  pulse  is 


*For  specifications  and  details  of  operation,  consult  RCA  Solid 
State  Data  Book  SSD-203A,  1973  edition. 
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removed  on  the  second  clock  pulse  and  integration  begins.   The 
amplifier  continues  integration  until  the  reset  pulse  goes  low  again 
at  pulse  32.   Gates  J  and  K  control  the  sampling  of  the  electrometer 
output.   On  clock  pulse  4  NOR  gate  K  goes  low  causing  outputs  S  and  U 
to  rise  from  0  to  5  volts.   This  level  change  initiates  the  convert 
command  to  the  analog-to-digital  (A/D)  converters  and  enables  the 
storage  registers  to  accept  the  converted  data.   NAND  gate  J  goes 
high  on  clock  pulse  28,  initiating  a  second  convert  command  and 
enabling  the  storage  registers  to  accept  a  second  set  of  converted  data. 

All  the  outputs  S,  T,  U  and  W  are  inputs  to  retriggerable 
monostable  multivibrators.*  These  devices  are  shown  on  a  later  figure. 
The  operational  amplifiers  M  and  N  serve  only  as  interface  devices 
to  convert  the  outputs  of  gates  J  and  K  into  voltage  levels  compatible 
with  the  1-shots.   Table  A  gives  a  truth  table  for  the  amplifier  logic 
shown  in  Figure  18. 

The  level  inverter  E,  in  Figure  18,  is  used  as  a  buffer,  delaying 
the  clock  pulses  by  one-half  cycle.   This  ensures  that  the  motor  moving 
the  ionization  chamber  completes  its  step  before  the  electrometer  output 
is  sampled.   The  action  in  counters  F  and  G  occurs  only  on  the  negative- 
going  pulse. 

Figure  20  shows  a  diagram  of  a  complete  measuring  cycle,  which  is 
divided  into  several  sections  to  illustrate  the  separate  events 
occurring  during  the  total  measuring  interval.   Although  each  measuring 
interval  extends  over  32  clock  pulses  and  a  time  interval  of  0.125 
seconds,  the  actual  interval  over  which  integration  is  performed  is 


*Type  TTuL9601,  manufactured  by  Fairchild  Semiconductor  and 
commonly  called  1-shots. 
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somewhat  smaller.   The  electrometer  output  Is  sampled  at  two  points 
C  and  D,  shown  in  Figure  20,  so  that  the  difference  of  the  two  results 
may  be  used  as  the  measure  of  the  average  current  in  the  interval. 
The  process  of  taking  the  difference  eliminates  any  problem  of  zero 
drift,  which  might  be  present  in  the  amplifier.   The  points  C  and  D 
are  symmetric  with  respect  to  the  extremes  of  the  measuring  interval 
so  that  scans  across  the  radiation  field  in  opposite  directions  will 
produce  data  corresponding  to  the  same  point  in  space  across  the  field. 
The  actual  interval  of  integration  is  then  that  between  points  C  and  D. 
Table  5  lists  the  elapsed  clock  pulses,  elapsed  time,  and  the  corre- 
sponding distance  traveled  at  each  of  the  points  labeled  in  Figure  20. 
The  results  in  Table  5  were  used  to  generate  Table  6,  which  lists  the 
corresponding  data  for  each  subinterval  between  the  points  in  Figure 
20.   Each  complete  measuring  interval  begins  with  the  integrator  being 
reset  and  remaining  reset  for  7.8125  milliseconds,  which  is  adequate 
time  for  the  capacitor  to  discharge  and  for  any  switching  transients 
to  disappear.   The  reset  pulse  is  removed  at  point  B,  integration 
begins  and  continues  without  interruption  from  point  B  to  point  E  with 
the  accumulated  value  being  sampled  at  two  points  during  the  interval 
BE.   The  first  measurement  is  made  at  point  C,  7.8125  milliseconds 
after  integration  begins.   The  second  measurement  is  made  at  point  D. 
The  total  time  of  integration  between  points  C  and  D  is  93.75  milli- 
seconds and  represents  the  measurement  time,  since  the  difference 
between  the  two  measurements  is  used  to  approximate  the  input  current 
midway  between  these  points.   The  ionization  chamber  moves  0.254 
centimeters  during  the  time  interval  between  points  C  and  D. 

Integration  continues  after  the  second  measurement  for  a  period 
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of  15.625  milliseconds  before  the  integrator  is  reset  at  point  E. 

The  interval  DE  is  of  the  proper  length  to  make  the  measuring  interval 

CD   symmetric  about  the  center  of  the  complete  measuring  cycle  AE,   The 

ionization  chamber  moves  0.3387  centimeters  during  a  complete  measuring 

cycle.   The  measuring  cycles  are  continually  repeated  an  integral 

number  of  times,  as  the  ionization  chamber  moves  across  the  radiation 

field. 

Amplifier  Testing 

After  construction  was  completed,  the  repetitive  integral  amplifier 
was  tested  for  linearity  and  stability.   During  all  testing  the  ampli- 
fier was  continually  cycling  under  the  control  of  the  system  clock 
and  the  difference  of  the  two  voltage  measurements  performed  in  each 
measuring  cycle  was  displayed.   The  A/D  converter  at  the  amplifier 
output  has  12-bit  resolution  and  thus  the  full  scale  amplifier  output 
voltage  (10  volts)  corresponds  to  a  displayed  value  of  A095.   A  single 
digit  of  the  display  represents  approximately  2.4  millivolts. 

Linearity  Testing 

For  this  test  procedure,  a  current  source  was  connected  at  the 
input  of  the  electrometer  amplifier.   This  required  that  the  test  be 
performed  before  the  amplifier  housing  was  filled  with  ceresin  wax 
and  before  the  ionization  chamber  was  connected  at  the  amplifier  input. 

The  current  source  was  obtained  by  connecting  one  lead  of  a 
high-ohmic  resistor  to  the  electrometer  input  and  the  other  lead  to  a 
highly  stabilized  adjustable  power  supply.   The  power  supply  was  varied 
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in  steps  and  the  response  of  the  repetitive  integral  amplifier  to 
the  corresponding  input  current  was  recorded. 

Both  repetitive  integral  amplifiers  were  tested  for  linearity  by 
this  method  with  amplifier  gain  ranging  from  unity  to  the  maximum. 
Figure  21  shows  that  both  amplifiers  have  a  linear  response  over  the 
full  range  of  amplifier  output  voltages  and  over  the  full  range  of 
gain. 

Stability  Testing 

The  ionization  chamber  was  connected  at  the  electrometer  input 
and  the  amplifier  housing  was  filled  with  ceresin  wax  before  stability 
tests  were  performed.   In  order  to  observe  long-term  resetting  stability 
the  amplifier  compensating  pulse  was  adjusted  so  that  the  result  of 
the  first  voltage  reading  of  the  measurement  cycle  was  about  100  on  the 
display  with  no  irradiation  of  the  chamber.   The  amplifier  was  allowed 
to  run  through  its  cycles  although  only  the  first  voltage  reading  was 
displayed  rather  than  the  difference  of  the  two  measurements.   The 
output  display  was  observed  four  times  each  hour  during  several 
eight-hour  days.   The  displayed  values  varied  no  more  than  ±1  digit 
in  the  least  significant  digit  during  the  test.   The  amplifier  was 
turned  off  each  night  and  restarted  the  next  day.   Within  1  to  2 
minutes  after  turn-on  the  amplifier  output  display  was  within  ±1  digit 
of  the  output  value  the  previous  day.   This  test  was  performed  with 
the  ionization  chamber  and  amplifier  housing  in  a  fixed  position  with 
no  movement  during  the  test.   The  resetting  characteristics  were 
similarly  evaluated  under  scanning  conditions  by  attaching  the  amplifier 
and  ionization  chamber  to  the  mechanical  scanning  unit  (to  be  discussed 
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in  the  following  chapter)  and  scanning.  The  deviations  again  did 
not  exceed  ±1  digit.  The  resetting  network  is  thus  stable  and  is 
not  affected  by  movement  of  the  amplifier. 

The  reproducibility  of  the  normal  output  data  (difference  of  2 
measurements  in  each  cycle)  was  evaluated,  again  with  no  irradiation 
of  the  chamber.   A  total  of  3,072  data  values  were  recorded  during 
scanning.   The  maximum  observed  deviation  from  the  mean  of  all  values 
was  1.5  digits.   The  3a  error  was  1.8  digits  and  corresponds  to  a 
voltage  deviation  of  4.4  x  10~^  volts  which  is  0.05  percent  of  the  full 
scale  value  corresponding  to  100  percent  or  maximum  radiation  exposure. 
Ninety-nine  percent  of  all  deviations  should  fall  within  the  3a  error 
and  thus  fluctuations,  which  are  due  to  noise,  are  small.   The  same 
test  was  performed  but  with  the  amplifier  fixed.   The  deviations 
observed  in  this  case  were  of  the  same  order  and  the  mean  of  all 
observations  was  not  significantly  different  at  the  1  percent  level 
from  the  mean  of  the  values  obtained  during  scanning.   Thus  movement 
and  vibration  of  the  amplifier  and  ionization  chamber  did  not  change 
the  amplifier  output. 

Electrical  Noise 

Any  electrical  signal  consists  of  two  components,  the  information 
being  transmitted  and  noise.   A  noise  component  is  always  present  and 
may  arise  from  many  sources  both  man-made  and  natural.   Electrical 
noise  imposes  the  ultimate  limit  on  the  performance  of  any  communications 
system.   Several  important  sources  of  noise  are  discussed  below. 
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Power  Hum 

Power  hum  consists  of  the  AC  power  frequency  (usually  60  cycles) 
with  harmonics  and  is  a  familiar  constituent  of  amplifier  outputs. 
Hum  results  when  extraneous  voltage  signals  are  coupled  into  the 
circuit  from  the  surroundings.   This  type  of  noise  may  be  greatly 
reduced  by  properly  shielding  the  circuit  with  a  grounded  conducting 
shield. 

The  housing  for  the  repetitive  integral  amplifier  is  made  from 
metal  and  is  connected  to  voltage  common  of  the  amplifier  and  thus 
serves  as  a  very  effective  shield.   The  amplifier  output  passes 
through  a  shielded  cable  to  the  A/D  converter  input.   Highly  stablized 
power  supplies  are  used  to  power  the  amplifier  and  thus  no  hum  origi- 
nates at  this  source.   With  proper  precautions  power  hum  may  be 
eliminated  to  the  extent  that  it  becomes  unimportant  as  a  noise  source. 

Microphonics 

Microphonic  noise  results  when  noise  voltages  are  generated  by 
movement  of  electrical  components  and  connecting  wiring.  This  type 
of  noise  is  particularly  important  if  the  circuit  will  be  subjected 
to  significant  vibration. 

The  mechanical  scanning  apparatus  is  driven  by  stepping  motors  at 
a  frequency  of  256  hertz.   This  means  the  amplifier  and  ionization 
chamber  are  accelerated  and  deaccelerated  256  times  each  second  during 
scanning,  a  fact  which  made  it  necessary  to  take  special  precautions 
to  avoid  microphonic  noise. 

The  amplifier  housing  was  filled  with  ceresin  wax  for  the  purpose 
of  displacing  all  air  surrounding  the  amplifier.  The  wax  serves  another 
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purpose,  that  of  firmly  holding  all  circuit"  components  and  wiring. 
The  remaining  source  of  microphonic  noise  is  the  ionization  chamber 
itself.   The  center  electrode  projecting  into  the  chamber  volume  is 
about  20.8  millimeters  in  length  and  is  free  at  one  end.   Any  movement 
of  the  center  electrode  causes  a  change  in  chamber  capacity  and  thus 
Induces  noise  at  the  amplifier  input.   In  order  to  reduce  electrode 
flexing  the  entire  center  electrode  was  fabricated  from  a  solid 
aluminum  rod  2.5  millimeters  in  diameter.   The  section  of  the  rod 
projecting  into  the  chamber  is  1.7  millimeters  in  diameter. 

The  previous  section  on  amplifier  testing  demonstrated  that  any 
microphonic  noise  present  is  not  of  sufficient  magnitude  to  cause  a 
change  in  the  amplifier  output.   The  deviations  observed  in  the 
amplifier  output  were  found  to  be  identical  when  measurements  were 
made  with  and  without  scanning  thus  implying  that  microphonic  noise  is 
unimportant. 

A  theoretical  treatment  of  microphonics  of  ionization  chambers 
has  not  been  included  in  this  discussion  but  such  a  treatment,  based 
on  the  mathematical  theory  of  elasticity,  has  been  considered  by 
Kilbourne.^S 

Amplifier  Drift 

Drift  is  classified  as  the  noise  of  very  low  frequency  in  the 
near-infinite  spectrum  of  frequencies.   Prolonged  tests  by  operational 
amplifier  manufacturers  have  shown  that  at  room  temperature  the 
long-term  drift  over  many  months  averages  near  zero  and  does  not  show 
any  steady  progression  or  systematic  pattern.   Drift  may  have  a  period 
of  days  or  weeks  or  even  longer. 
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The  repetitive  integral  amplifier  makes  2  measurements  during 
each  measuring  cycle  and  the  difference  of  the  2  values  is  recorded 
as  the  result  corresponding  to  the  value  of  the  input  current  at  the 
center  of  the  measuring  interval.   By  taking  the  difference,  drift 
with  a  period  much  greater  than  the  measuring  interval  is  eliminated. 
No  zero  adjustment  is  therefore  required  on  the  amplifier  and  drift 
need  not  be  considered. 

Input  Leakage 

Leakage  currents  reaching  the  electrometer  input  are  of  particular 
Importance  when  very  small  input  currents  are  measured.   Any  leakage 
current  which  varies  with  temperature  or  other  parameters  constitutes 
a  noise  source.   In  order  to  avoid  leakage  into  the  electrometer  input 
particular  care  is  taken  to  insulate  the  input  from  other  parts  of  the 
circuit.   The  electrometer  amplifier  (ICH8500A)  shown  in  Figure  14  is 
mounted  on  a  printed  circuit  board.   A  hole  is  placed  in  the  circuit 
board  at  the  electrometer  input  lead  so  that  the  lead  can  pass  through 
the  board  without  contacting  it.   Likewise  the  source  lead  of  the 
transistor  which  connects  with  the  input  passes  through  the  circuit 
board  without  contacting  it  and  is  joined  to  the  electrometer  input. 
The  feedback  capacity  and  the  protection  diodes  shown  in  Figure  14 
connect  at  the  electrometer  input  without  touching  the  printed  board. 

Ground  Loops 

Ground  loops  are  generated  when  a  chassis,  frame  or  bus  is  used 
as  a  common  ground  for  two  or  more  pieces  of  equipment.   No  two  ground 
points  have  the  same  potential  and  more  than  one  ground  on  a  signal 
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circuit  or  shield  creates  a  current  flow  between  grounds.   The  current 
flow  through  a  ground  loop  can  produce  electrical  noise  in  any  low- 
level  circuitry. 

In  an  effort  to  avoid  ground  loops  all  power  supply  cotnmon 
terminals  were  connected  at  one  point  and  optical  couplers  (discussed 
in  the  following  chapter)  were  used  to  isolate  the  circuits  of  the 
amplifier,  amplifier  logic  and  scanning  logic  from  circuits  located 
in  the  data  recording  unit.   Isolation  of  the  two  units  is  necessary 
because  they  may  draw  their  AC  power  from  different  electrical  circuits 
with  differing  ground  potentials.   The  two  units  are  located  in  separate 
rooms  during  scanning. 

Random  Noise 

Even  when  there  is  no  input  to  an  amplifying  system  random 
fluctuation  of  the  output  voltage  is  observed.   These  output  variations 
are  due  to  3  types  of  random  noise  mechanisms  usually  referred  to  as 
thermal  noise,  low-frequency  (1/f)  noise  and  shot  noise. 

Thermal  noise  is  inherent  in  all  systems  operating  at  a  temperature 
exceeding  absolute  zero  and  is  caused  by  the  thermal  excitation  of  the 
charge  carriers  in  a  conductor.   The  thermal  agitation  of  electrons 
in  a  conductor  produces  random  current  pulses  which  constitute  noise. 
The  noise  voltage  or  current  has  a  gaussian  probability  density  function 
and  its  spectrum  has  been  shown  to  have  a  uniform  distribution.^^ 
That  is,  the  power  density  versus  frequency  is  flat  and  thus  it  is 
said  to  be  white  noise. 

Noise  with  low-frequency  spectral  characteristics,  usually  called 
1/f  noise,  is  observed  in  tubes,  transistors,  diodes  and  resistors. 
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This  type  noise  is  traceable  to  properties  of  the  surface  of  the 
material.   It  is  thus  highly  dependent  on  manufacturing  techniques 
and  varies  from  one  device  to  another  even  with  devices  of  the  same 
type.   The  spectral  density  of  1/f  noise  increases  as  the  frequency 
decreases.   The  noise  power  is  equal  in  each  decade  of  frequency  and 
thus  1/f  noise  is  not  white. 

Shot  noise  is  observed  in  tubes,  transistors  and  diodes  and  is 
the  result  of  a  discontinuous  current  flowing  in  the  devices.   A  current 
is  the  result  of  movement  of  discrete  carriers  each  with  one  electronic 
charge.   The  random  arrival  of  these  discrete  carriers  at  a  point  gives 
rise  to  shot  noise,  which  has  a  discrete  distribution  (Poisson).- 

Even  if  it  were  possible  to  eliminate  all  other  sources  of  noise, 
thermal  noise  contributes  an  irreducible  minimum  amount  of  noise 
power  to  the  output  of  an  amplifier.   It  is  thermal  noise  which 
ultimately  limits  the  fidelity  of  amplification  when  the  input  signal 
is  small. 

The  true  noise  output  of  an  amplifier  and  the  thermal  noise 
output  are  related  together  by  what  is  called  the  noise  figure.   The 
noise  figure  F  of  a  two-port  linear  device  is  the  ratio  of  the  total 
available  output  noise  power  per  unit  bandwidth  to  the  total  available 
output  noise  power  due  to  thermal  noise. 


„    Nactual 
r  = 


^thermal 

N  ^  1  is  the  sum  of  the  noise  power  contributed  by  1/f  noise,  shot 
actual  '^ 

noise,  thermal  noise  and  any  other  type  noise  present. 
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As  F  depends  upon  the  thermal  noise,  it  therefore  varies  with 
temperature.   With  the  ideal  device,  the  only  noise  is  thermal  due  to 
the  device  being  in  an  environment  in  which  the  temperature  is  above 
absolute  zero.   Thus,  the  noise  figure  of  an  ideal  amplifier  is  unity.. 
Any  additional  noise  sources  within  the  device  cause  F  to  exceed  unity 
and  thus  the  excess  of  F  over  unity  represents  a  degradation  from  the 
ideal  case.   The  noise  powers  N  are  often  measured  over  a  band  of 
frequencies  centered  on  some  operating  frequency  of  interest  and  the 
resulting  noise  figure  F  is  called  a  spot  noise  figure  as  opposed  to 
an  average  noise  figure  taken  over  a  wider  range  of  frequencies. 

The  noise  figure  is  principally  used  as  a  f igure-of-merit  for 
comparison  of  amplifiers  and  is  not  necessarily  the  appropriate 
indicator  for  optimizing  overall  performance  since  the  amplifier 
gain  may  not  be  at  a  maximum  at  the  point  where  the  noise  figure  is 
at  a  minimum. 

Amplifier  Noise 

Any  noise  response  analysis  of  the  repetitive  integral  electrome- 
ter may  be  limited  to  an  analysis  of  the  noise  response  of  the 
electrometer  amplifier  as  this  is  where  the  input  signal  and  input 
noise  is  amplified.   The  electrometer  amplifier  performs  the  function 
of  integration. 

Figure  19  gives  a  diagram  of  a  single  amplifier  measuring  cycle. 
The  voltage  generated  at  the  amplifier  output  is  sampled  two  times 
during  each  cycle,  the  first  time  shortly  after  integration  begins 
and  the  second  0.09375  seconds  later.   The  amplifier  measuring  cycles 
are  continuously  repeated  during  scanning,  each  cycle  requiring 
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0.125  seconds.   The  difference  of  the  two  voltage  samples  is  recorded 
and  is  used  to  estimate  the  average  current  during  the  measuring 
interval. 

For  the  purposes  of  noise  analysis,  it  is  sufficient  to  consider 
a  single  measuring  cycle  and  use  an  ensemble  average  to  obtain  the 
RMS  error  in  the  data  points  due  to  the  noise. 

Considering  a  single  measuring  cycle,  assume  that  integration 
begins  at  time  tn   and  that  the  two  voltage  samples  are  made  at  times 
t^  and  t2.   Also  assume  that  the  amplifier  input  current  during  the 
measuring  cycle  is  given  by  I(t)  for  tQ  <  t  s  tn  +  At  where  At  is 
the  time  between  voltage  samples.   If  X  represents  the  difference  of 
the  two  sample  points  then 


X  =  W   I(t)  dt  -  i  J   I(t)  dt 


■0 

where  C  is  the  feedback  capacity.   The  current  I(t)  =  Is(t)  +  Ifi^'^^ 
where  IgCt)  is  the  signal  current  from  the  ionization  chamber  and 
In(t)  is  the  current  produced  by  random  noise.   Since  integration  is 
a  linear  operation  it  is  possible  to  separate  X  into  a  signal  term 
and  a  noise  term.   If  x  is  that  portion  of  the  output  voltage  due  to 
noise  then 

t2  t^ 

X  =^  I   In(t)  dt  -  ^  I   In(t)  dt 

The  voltage  has  only  an  AC  component  since  any  DC  noise  is  cancelled 
out  in  the  subtraction.  Thus,  the  variance  is  given  by 
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o2  =  E(x2) 


where  E(x^)  denotes  the  ensemble  average  of  x^ 


E(x2)  =  E{[^  I  I„(t)  dt  -  ^   I   In(t)  at]2} 

•'to  •'to 


=  -^E{[  I   I„(t)  dt]2}  _  -^  E{   j   In(t)  dt   I   In(t)  dt} 
J     ^0  ^     J  tg        J  to 

+^  E{[   I  In(t)  dt]2} 

/t2   #»t2 
I  E[In(t)  In(t')]  dt  dt' 
^0  ^    ^0 

2     f'2  f  '1 
-  77   I     I   E[In(t)  In(t')]  dt  dt'  (37) 

^    J  to  J  to 

1     f   f'l 

+  ^    I     I  E[In(t)  In(t')]  dt  dt' 

J     toJ  to 


When  the  measuring  process  is  ergodic  the  autocorrelation  function 
of  two  jointly  distributed  gaussian  random  variables  N(t)  and  N(t')  is 
equal  to  the  ensemble  average.   In  this  case,  if  y  =  N(t)  and  z  =  N(t') 


/7 


E(yz)  =  R(t  -  t')  =    I     I  yz  p(y,z)  dy  dz 

where  p(y,z)  is  the  joint  probability  density  function  and  R(t  -  t')  is 
the  autocorrelation  function. 

Then  the  ensemble  averages  in  equation  (37)  may  be  replaced  by  the 
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autocorrelation  function. 
Thus 


1  r^  f '2  2  r^  f  1 

E(x2)    =  ^       I  I  R(t   -    t')    dt   dt'    -  ^      I  I        R(t    -    t')    dt   dt' 

Jto     Jto  ^       Jto      Jto 


+  ^   I    I  R(t  -  t')  dt  dt' 


The  autocorrelation  function  and  the  power  spectral  density  are  Fourier 
transform  pairs  and  thus  the  variance  for  any  random  noise  may  be 
calculated  from  equation  (38)  if  either  the  autocorrelation  function 
or  the  power  spectral  density  is  known  for  the  noise  process. 

The  power  spectral  density  for  white  noise  is  known  to  be  constant 
over  frequencies  of  interest  in  noise  analysis.   If  the  power  spectral 
density  is  denoted  by  n,  a  constant,  then  the  autocorrelation  function 
of  white  noise  is  given  by 

R(t  -  t')  =  n  6  Ct  -  t')  (39) 

where  5(t  -  t')  is  the  delta  function  defined  as 


H 


,,   /  0  t  *  t' 

t  -  f)  =  t  „  t  =  t' 


The  area  under  the  delta  function  is  unity. 

Substituting  equation  (39)  into  equation  (38)  and  integrating,  the  noise 

power  at  the  amplifier  output  is 


E(x2)  =  -gz  (to  -  t^) 
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and  the  RMS  noise  voltage  due  to  thermal  noise  is 

a   =  yE(x^)   =  ^yn(t2  -  t^)  (40) 

It  is  well  documented  experimentally  as  well  as  from  quantum 
mechanical  considerations  that  the  available  thermal  noise  power  is 
given  by 

P  =  k  T  B  (41) 

where  k  =  1.38  x  10~^^  watt-seconds  is  Boltzman's  constant,  T  is  the 
temperature  in  Kelvin  and  B  is  the  noise  bandwidth  of  the  measuring 
system. ^^ 

The  ionization  chamber  supplying  current  to  the  amplifier  input 
may  be  represented  by  the  equivalent  circuit  below 


where  A  is  an  ideal  current  source  supplying  a  current  I  and  R  is  the 
internal  resistance  of  the  chamber.   The  available  thermal  noise  power 
given  by  equation  (41)  is  dissipated  in  a  matched  load  condition 


P  =  k  T  B  =  14   R 


■I 


The  value  I  may  be  considered  as  normalized  power  P 
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,   4  k  T  B 


Thus,  in  the  case  of  thermal  noise,  a  white  noise,  the  power  spectral 
density  n  is 


2  k  T 
n  =  — 5 


The  minimum  effective  feedback  capacity  given  in  Table  3  is  0.76 
X  10"^^  farads  and  the  time  difference  between  the  two  measurements 
in  each  cycle  is  0.09375  seconds.   The  internal  resistance  of  the 
ionization  chamber  was  determined  to  be  10^ ^  ohms.  Thus,  from  equation 
(40),  the  RMS  noise  voltage  due  to  thermal  noise  is 


. I V 

0.76  X  10"12  y 


2(1.38  X  10  -')    (290)  (0.09375) 
IQlS 


=  I.IA  X  10"^  volts 

Thus,  thermal  noise  is  negligible  since  the  values  recorded  (corre- 
sponding to  the  difference  of  the  2  measurements)  exceed  800  millivolts 
in  this  application. 

Summary 

The  repetitive  integral  measurement  method  approximates  the  input 
function  by  a  series  of  integrations  performed  over  short  time  intervals. 
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Each  period  of  integration  is  independent  of ■  any  other,  which  means 
that  the  integrating  capacitor  must  be  discharged  to  zero  voltage  before 
integration  begins  in  each  interval. 

Problems  encountered  when  using  traditional  methods  to  reset  the 
integrator  were  overcome  through  the  use  of  dual  field-effect  transistors. 

The  cycles  of  integration  and  reset  are  controlled  and  correlated 
with  ionization  chamber  movement  by  a  system  clock  and  digital  logic. 
A  single  measuring  cycle  extends  over  a  time  of  0.125  seconds  during 
which  time  two  measurements  are  made,  the  first  shortly  after  integration 
begins  and  the  second  0.09375  seconds  later.   The  difference  of  the  two 
measurements  is  recorded  and  is  used  as  a  measure  of  the  average  -current. 

Actual  run  testing  has  found  the  amplifier  to  be  linear,  stable  and 
free  of  microphonics.   A  noise  analysis  has  shown  that  thermal  noise 
does  not  present  a  problem  with  the  repetitive  integral  amplifier. 

The  following  chapter  discusses  the  total  depth  dose  scanning 
system. 


CHAPTER  7 
THE  DEPTH  DOSE  SCANNING  SYSTEM 

The  development  of  the  repetitive  integral  measurement  method  and 
its  implementation  in  the  repetitive  integral  electrometer  amplifier 
has  made  it  possible  to  eliminate  the  most  serious  problem  encountered 
in  depth  dose  scanning,  that  of  the  time  constant  error.   As  measuring 
accuracy  is  increased  with  the  new  method,  precise  positioning  of  the 
detector  and  reproducibility  of  movement  on  each  scan  become   important. 
Previously  used  methods  of  scanner  control  by  limit  switches  are  not 
sufficiently  accurate  to  ensure  precise  detector  positioning.   Scanner 
drive  mechanisms,  which  allow  variation  in  scanning  speed,  may  not  be 
used  with  a  measuring  system  whose  cycles  are  based  on  time,  as  no 
correlation  exists  between  measurement  cycles  and  actual  detector 
position.   Although  the  application  of  the  repetitive  integral  method 
eliminates  the  most  serious  problem  encountered  in  depth  dose  scanning, 
it  introduces  other  problems  which  necessitate  redesign  of  mechanical 
scanning  components.   Furthermore,  the  new  measurement  method  requires 
a  more  sophisticated  system  for  data  handling,  since  arithmetic  processing 
is  performed  on  the  two  measurements  made  in  each  measurement  cycle. 

The  complete  depth  dose  scanning  system  may  be  considered  as  a 

combination  of  three  subsystems:  the  mechanical  scanning  system,  the 

data  processing  and  recording  system,  and  the  detector  and  amplifier 

system. 
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The  Mechanical  Scanning  System 

Measuring  cycles  of  the  repetitive  integral  amplifier  are  governed 
by  pulses  which  originate  from  an  accurate  crystal  clock.   The 
measuring  cycles  are  thus  strictly  time  dependent.   If  good  correlation 
of  the  measuring  cycles  and  detector  position  is  to  be  maintained,  the 
mechanical  scanning  movement  must  also  be  time  dependent.   This  means 
that  the  length  of  each  individual  lateral  scan  must  be  identical  and 
the  speed  of  scanning  must  always  be  the  same  on  each  lateral  scan. 
Perfect  correlation  of  the  movement  and  measuring  cycles  may  be  main- 
tained if  the  movement  is  also  controlled  by  the  clock  pulses.   Motors 
which  respond  to  pulses  rather  than  running  continously  are  called 
stepping  motors. 

Stepping  Motors 

Stepping  motors  are  specially  designed  motors  which  rotate  through 
a  particular  angle  upon  application  of  a  pulse,  and  then  stop  and  hold 
their  position  until  another  pulse  is  received.   The  stepping  motor  may 
be  likened  to  a  series  of  electromagnets  or  solenoids  arranged  in  a 
circle.   When  sequentially  energized,  the  electromagnets  or  solenoids 
react  with  a  rotor,  causing  it  to  turn  in  either  direction,  depending 
on  the  input  command.   The  stepping  angle  is  determined  by  the  design 
of  the  motor  and  may  be  as  large  as  90  degrees.   The  speed  at  which  a 
stepping  motor  can  run  without  missing  a  step  is  determined  by  the  load 
on  the  motor.   When  the  motor  is  accelerated  beyond  the  maximum  step 
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response  the  motor  enters  the  slewing  mode  and  is  no  longer  capable  of 
instantaneous  stopping  or  change  of  direction  without  loss  of  step 
integrity.   When  the  motors  are  energized  but  are  not  receiving  pulses, 
they  strongly  resist  movement  and  do  not  drift  so  that  position  is 
maintained.   As  long  as  pulse-to-step  integrity  is  maintained,  the 
stepping  error  is  noncummulative.   After  a  stream  of  pulses  is  fed  to 
a  stepping  motor,  its  final  shaft  position  is  known  to  within  a  small 
percentage  of  a  single  step. 

The  close  correlation  between  detector  position  and  measuring 
cycles  may  easily  be  obtained  through  the  use  of  stepping  motors  to  drive 
the  detector  assembly.   Any  stepping  error  is  small  and  the  scan  speed 
Is  always  the  same  since  the  clock  frequency  is  constant. 

Two  identical  stepping  motors*  were  chosen  for  the  x-drive  and  the 
y-drive,  the  former  moving  the  detector  laterally  across  the  radiation 
field  and  the  latter  in  increasing  depths.  These  particular  motors 
are  approximately  1  inch  in  diameter  and  2  inches  in  length  and  were 
chosen  because  of  their  small  size  and  stepping  characteristics.   The 
motors  rotate  15  degrees  with  each  input  pulse,  have  a  maximum  stepping 
rate  of  900  steps  per  second  under  no  load  and  will  rotate  in  either 
direction  as  desired. 

Drive  Assembly 

Since  the  chosen  stepping  motors  are  small,  they  have  a  limited 
load  capacity.   It  was  essential  to  keep  the  drive  assembly  as  friction 
free  as  possible  and  therefore  ball  screws  were  chosen  for  both  the 


♦Manufactured  by  IMC  Magnetics  Corporation,  Maywood,  California. 
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X  and  y-drive  to  reduce  friction  since  they  have  large  threads  and  the 
corresponding  ball  nut  moves  up  and  dovm  the  threads  on  a  cyclic  stream 
of  ball  bearings.   This  type  of  screw  has  so  little  friction  that  the 
ball  nut  moves  freely  down  the  threads  solely  under  the  influence  of 
gravity.   Their  design  is  such  that  movement  in  either  direction  requires 
the  same  torque.   The  ball  screws  chosen  are  about  0.3  inches  in  diameter 
and  are  about  22  inches  long,  the  large  diameter  being  essential  for 
stability.   The  screws  were  connected  directly  to  the  stepping  motors, 
no  gearing  being  required.   Twenty-four  pulses  to  the  motors  produce  one 
complete  rotation  of  the  ball  screw,  during  which  the  ball  nut  moves 
along  the  screw  a  certain  distance  called  the  screw  "lead."  The  only 
appropriate  ball  screws  which  could  be  obtained  had  a  lead  of  0.1  inch. 
The  mechanical  scanning  apparatus  consists  of  the  ball  screws  with 
stepping  motors.   The  ball  nuts  for  the  x  and  y-drive  are  connected  to 
the  detector  mounting  bracket  and  to  the  assembly  on  which  the  x-drive 
is  mounted,  respectively.   Spline  rods  are  used  with  each  drive  for 
stability.   Figure  22  shows  the  mechanical  scanning  apparatus  with  the 
drive  screws  and  motors. 

The  Scanning  Logic 

The  purpose  of  scanning  logic  is  to  provide  controlled  pulses  to 
the  stepping  motors  so  that  a  particular  pattern  of  scanning  is  negoti- 
ated by  the  detector.   The  scanning  logic  maintains  the  required 
correlation  between  scanner  movement  and  measuring  cycles  and  acts  as 
a  special  purpose  digital  computer  in  its  controlling  function. 

This  special  purpose  computer  was  programmed  to  automatically 
direct  the  movement  of  the  detector  in  lateral  scans  across  the  radiation 
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field  with  each  lateral  scan  being  separated  by  a  fixed  distance.   After 
reaching  a  fixed  depth  in  the  radiation  field,  the  detector  is  moved 
back  to  its  original  starting  position.   Pulses  going  to  the  drive  motors 
are  counted  so  that  each  lateral  scan  is  identical  in  length  and  the 
distance  between  scans  is  kept  constant.   By  this  digital  control  method, 
absolute  correlation  between  amplifier  measuring  cycles  and  detector 
movement  is  possible.   A  selection  of  several  scan  widths  and  scanning 
speeds  is  available.   The  ability  to  select  a  scanning  width  is  par- 
ticularly important  since  scanning  time  can  be  reduced  by  selecting  the 
minimum  scan  width  which  adequately  covers  the  radiation  field. 

The  action  of  the  scanning  logic  may  be  summarized  as  follows. 
Automatic  scanning  is  manually  initiated  with  a  start  switch  at  which 
time  pulses  are  immediately  directed  to  the  x-motor  controller.   The 
motor  controller  generates  the  high  current  pulses  required  by  the  motor 
and  also  governs  the  direction  of  rotation  depending  on  which  input  to 
the  controller  is  pulsed.   When  the  x-motor  receives  pulses,  the  detector 
moves  in  steps  until  the  number  of  pulses  corresponding  to  the  chosen 
scan  width  is  counted  at  which  time  pulses  to  the  x-motor  are  discon- 
tinued and  the  y-motor  is  pulsed  to  move  the  detector  to  a  greater  depth. 
The  x-motor  holds  its  last  position  while  the  y-motor  is  pulsed  until 
the  number  of  pulses  corresponding  to  the  depth  increment  is  counted. 
The  y-motor  pulses  are  then  stopped  and  the  x-motor  is  again  pulsed 
but  the  direction  of  travel  is  reversed.   The  x-motor  receives  the 
same  number  of  pulses  as  before  and  again,  the  y-motor  is  pulsed  to 
move  to  greater  depth.   This  process  is  repeated  an  even  number  of  times 
after  which  the  detector  is  returned  to  its  original  starting  position. 

The  amplifier  measuring  cycles  are  continually  repeated  at  a 
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frequency  of  8  cycles  per  second  during  the  time  the  x-motor  is  stepping. 
The  amplifier  does  not  cycle  during  the  time  the  detector  is  moving 
to  a  greater  depth  nor  during  the  time  the  detector  is  returning  to  its 
original  starting  position  after  the  scan  is  completed. 

Specific  operation  of  the  scanner  logic  may  best  be  seen  by 
consulting  Figure  23  which  shows  the  scanner  logic  diagram.   Throughout 
this  diagram,  the  high  logic  level  is  zero  volts  and  the  low  logic  level 
is  -10  volts. 

The  heart  of  the  scanning  logic  as  well  as  the  amplifier  logic 
(discussed  in  Chapter  6)  is  the  crystal-controlled  system  clock  (U  in 
Figure  23)  which  runs  continuously  whenever  power  is  supplied  to  the 
system.   The  clock  frequency  of  16.384  kilohertz  is  far  too  high  to  drive 
the  motors  directly  and  is  reduced  by  a  14-stage  binary  counter  (V). 
Four  selectable  outputs  are  available  on  this  counter  so  that  pulse 
rates  of  512,  256,  128  and  64  pulses  per  second  may  be  used  to  drive  the 
stepping  motors.   Although  the  maximum  no-load  stepping  rate  of  the 
motors  is  900  steps  per  second,  this  is  reduced  when  a  load  is  present. 
With  the  frictional  and  inertial  load  introduced  by  the  drive  screw  and 
detector  assembly,  it  was  found  that  the  motors  did  not  step  reliably 
at  a  pulse  rate  of  512  pulses  per  second  but  performed  with  no  loss  of 
steps  at  256  pulses  per  second. 

When  the  scanning  logic  is  energized  by  connection  to  the  AC  line 
voltage,  the  silicon-controlled  rectifier  (SCR  in  Figure  23)  fires  and 
becomes  conducting.   This  causes  the  output  of  gate  W  to  go  high 
providing  a  high  logic  level  to  all  counter  reset  inputs,  labeled  R,  in 
both  the  scanner  logic  and  the  amplifier  logic  (Figure  18) .   This  causes 
all  counters  to  reset  to  their  zero  state.   Input  5  to  NAND  gate  S  is 
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low  and,  thus,  this  gate  is  disabled  and  does  not  pass  clock  pulses 
to  the  counters  or  motors.   The  system  remains  in  a  reset  condition 
until  start  switch  S6  is  connected  to  the  negative  side  of  the  relay 
power  supply  causing  relay  R6  to  open.  When  R6  opens,  the  current 
through  the  SCR  drops  and  it  ceases  to  conduct.   Gate  W  immediately  goes 
low  and  the  counters  are  enabled  for  counting.   Input  5  to  gate  S 
simultaneously  goes  high  and  the  gate  passes  clock  pulses  to  the  logic. 
Inverter  R  acts  as  a  buffer,  delaying  the  pulses  by  one-half  pulse  to 
ensure  that  the  counters  have  time  to  recover  from  their  reset  condition 
before  being  pulsed.   The  clock  pulses  pass  through  gate  A  to  counter  B 
which  counts  the  number  of  pulses  supplied  to  the  x-motor.  The  clock 
pulses  pass  through  gate  J  and  are  sent  to  the  x-motor  causing  the 
detector  to  move.   The  detector  continues  to  step  with  the  clock  pulses 
until  the  number  of  counts  corresponding  to  the  selected  scan  width  is 
reached  in  counter  B.   When  the  output  of  counter  B  goes  high,  the 
EXCLUSIVE-OR  gate  E  goes  high  and  disables  gate  A  so  that  no  further 
pulses  go  to  the  x-counter.   Gate  G  is  enabled  and  passes  clock  pulses 
to  counter  H  and  to  the  y-motor.   The  detector  steps  to  increasing  depth 
at  each  clock  pulse  until  the  Ay-counter  reaches  the  number  of  counts 
corresponding  to  the  depth  increment.   When  this  count  is  reached,  the 
output  of  the  Ay-counter  goes  high  and,  now  that  both  inputs  of  the 
EXCLUSIVE-OR  gate  are  high,  its  output  goes  low,  allowing  pulses  to  pass 
again  to  the  x-counter  and  disabling  gate  G  so  that  no  further  pulses 
reach  the  Ay-counter.   Since  the  Ay-counter  output  stays  high,  gate  J 
does  not  pass  pulses  to  the  x-motor  but  pulses  pass  gate  C  and  cause 
the  x-motor  to  reverse  direction.   When  the  appropriate  number  of  pulses 
are  counted  by  the  x-counter,  the  counter  output  goes  low  and  the 
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EXCLUSIVE-OR  gate  E  disables  gate  A  and  enables  gate  G  so  that  pulses 
again  go  to  the  Ay-counter  and  to  the  y-motor.   When  the  appropriate 
number  of  pulses  is  counted  into  the  Ay-counter  its  output  goes  low 
and  the  EXCLUSIVE-OR  gate  disables  gate  G  and  again  enables  gate  A  so 
that  pulses  go  to  the  x-counter  and  x-motor,  causing  movement  in  the 
direction  opposite  to  the  previous  scan.   This  sequence  of  events  is 
repeated  16  times  during  a  complete  radiation  field  scan.   Only  pulses 
which  pass  through  gate  A  are  sent  to  the  amplifier  logic  so  that  the 
amplifier  runs  through  its  measurement  cycles  only  during  lateral 
scanning  and  not  during  the  time  that  the  detector  moves  to  greater 
depths. 

All  pulses  which  go  to  the  Ay-counter  are  also  counted  by  the 
y-counter  (N)  which  registers  the  cummulative  total  of  all  pulses  sent 
to  the  y-motor  as  the  detector  moves  to  greater  depths.   The  output  of 
this  counter  goes  high  after  the  detector  has  been  incremented  in  depth 
33  times,  at  which  time  gate  0  is  enabled  allowing  pulses  to  pass  to 
the  y-motor  so  that  its  direction  is  reversed.   Pulses  are  also  counted 
by  the  y-return  counter  (P)  until  the  same  number  of  pulses  as  recorded 
in  the  y-counter  are  counted,  causing  the  y-motor  to  return  the  detector 
to  its  original  starting  position.   When  the  y-return  counter  reaches 
the  appropriate  value,  its  output  goes  high  causing  the  SCR  to  fire  and 
reset  the  system.   When  the  y-counter  goes  high,  gate  A  and  all  the  gates 
which  pass  pulses  to  the  x-motor  are  disabled  and  thus  no  pulses  are  sent 
to  the  amplifier  logic  on  the  x-motor  during  the  return  scan.   After 
the  last  lateral  scan,  the  detector  is  in  the  same  x-position  as  it  was 
at  the  original  starting  position  since  an  even  number  of  lateral  scans 
are  made.   Figure  23  shows  all  pulses  to  the  motors  passing  through 
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operational  amplifiers  of  the  741  type  discussed  in  Chapter  6.   The 
operational  amplifiers  are  used  to  change  the  negative  pulses  to  posi- 
tive pulses,  as  required  by  the  motor  controllers. 

The  automatic  scanning  sequence  may  be  stopped  at  any  time  during 
the  scan  by  closing  switch  S5  which  causes  relay  R5  to  fire  the  SCR, 
resetting  the  system.   Once  stopped,  the  detector  must  be  moved  manually 
back  to  the  starting  point  before  automatic  scanning  can  resume. 

Manual  movement  of  the  detector  is  possible  by  actuating  one  of 
switches  SI,  S2,  S3  or  54,  which  causes  corresponding  relays  to  connect 
the  clock  pulses  to  the  motor  required  for  the  desired  movement.   When 
any  one  of  the  switches  is  activated,  the  SCR  is  immediately  fired  and 
the  system  logic  goes  to  the  reset  condition.  The  system  stays  reset 
during  any  manual  movement.   Switch  57  provides  a  selection  of  two 
speeds,  normal  or  slow,  for  manual  operation.   When  the  slow  speed  is 
used,  the  normal  clock  frequency  of  256  pulses  per  second  is  further 
reduced  to  facilitate  precise  placement  of  the  detector. 

Switch  S8,  when  set  in  position  2,  allows  clock  pulses  to  reach 
the  amplifier  logic  while  the  scanner  is  not  in  operation  (reset  con- 
dition).  This  option  is  very  useful  for  amplifier  testing  or  when 
measurements  are  to  be  made  at  a  fixed  point. 

For  safety  purposes,  limit  switches  are  placed  at  each  end  of  the 
X  and  y-drive.   All  these  switches  are  connected  to  relay  R7  which, 
when  energized,  disconnects  the  pulse  source  to  all  motors.   In  normal 
automatic  operation,  the  limit  switches  are  never  engaged  but  in  manual 
operation,  the  extremes  of  each  drive  screw  may  be  reached.   When  the 
detector  assembly  hits  a  limit  switch,  the  motor  pulses  are  immediately 
stopped,  thus  preventing  damage  to  the  scanning  mechanism  and  the  motors. 
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The  x-counter  (B  in  Figure  23)  is  composed  of  several  logic 
components  and  a  switch  which  allows  selection  of  4  scan  widths,  chosen 
so  that  an  integral  number  of  amplifier  cycles  occurs  during  the  lateral 
scan.   The  amplifier  cycles  consist  of  32  pulses  and  so  the  total  of 
pulses  for  each  scan  width  is  an  integral  multiple  of  32.   The  largest 
scan  width  was  chosen  to  be  the  maximum  permitted  by  the  mechanical 
system.   Each  of  the  other  scan  widths  is  decreased  by  about  5.5  centi- 
meters so  that  the  difference  between  the  minimum  and  maximum  scan  width 
is  over  16  centimeters.   Table  7  shows  the  number  of  pulses  corresponding 
to  each  width,  the  time  required  for  the  lateral  scan,  the  time  required 
for  the  total  scan,  the  number  of  amplifier  cycles  (points  measured)  for 
that  width  and  the  field  sizes  which  can  be  covered  with  each  width. 

Figure  24  shows  the  complete  x-counter  circuit.   The  heart  of  the 
circuit  is  the  4018  divide-by-n  counter  which  allows  pulse  division  by 
any  digit  2  through  9  by  connecting  the  appropriate  output  to  the  data 
input  of  the  4018.   For  example,  division  by  4  is  accomplished  by 
connecting  output  Q2  to  the  data  input.   Then  the  Q2  output  goes  high 
after  every  fourth  pulse.   Four  switch  positions  are  available  to  divide 
by  4,  5,  6  or  7.  The  switch  output  is  connected  to  a  14-stage  binary 
counter  through  a  pulse  inverter.   The  pulse  inverter  is  used  in  this 
case  so  that  the  14-stage  counter  will  count  simultaneously  with  the 
divide-by-n  counter,  whose  output  goes  from  low  to  high  after  the  number 
of  pulses  corresponding  to  the  divisor  is  reached.   The  14-stage  counter 
counts  only  on  the  negative  going  pulse  so  that  the  inverter  reverses 
the  level  and  causes  the  counter  to  advance  as  the  divide-by-n  counter 
goes  from  low  to  high.   The  14-stage  counter,  which  goes  high  on  the 
512th  pulse,  is  the  output  of  the  x-counter  shown  in  Figure  23.   The 
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Figure  Ik.      Diagram  of  X-Counter. 
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output  stays  high  until  the  1,024th  pulse,  at  which  time  it  again  goes 
low.   The  divide-by-n  counter  has  the  facility  for  a  preset  count  but 
since  this  was  not  used  in  this  application  all  the  preset  inputs  and 
the  preset  input-enable  were  disabled  by  connection  to  logic  level  low. 

The  Ay-counter,  the  y-counter  and  the  y-return  counter  shown  in 
Figure  23  function  in  a  manner  similar  to  the  x-counter.   Figure  25 
shows  the  diagram  of  these  counters,  which  use  the  divide-by-n  counter 
to  divide  the  clock  pulses  by  6.   The  Ay-counter  uses  a  7-stage  binary 
counter  rather  than  the  14-stage  counter  used  in  the  other  circuit.   The 
output  of  this  7-stage  counter  goes  high  after  96  clock  pulses  and 
remains  high  until  192  pulses  are  counted  at  which  time  it  goes  low. 
The  depth  increment,  then,  consists  of  96  pulses  and  corresponds  to  a 
distance  of  1.016  centimeters.   Each  lateral  scan  is  separated  by  this 
depth.   Both  the  y-counter  and  the  y-return  counter  are  configured  as 
shown  in  Figure  25.   These  counters  go  high  after  3,072  pulses.   This  is 
the  total  number  of  pulses  accumulated  during  the  total  of  the  32  depth 
increments . 

All  gates  and  counters  used  in  the  scanning  logic  are  of  the  RCA 
COSMOS  4000  family  of  digital  integrated  circuits.   The  SCR  is  the  type 
C5A*. 

The  Water  Phantom 

The  water  phantom  is  a  tank  filled  with  water  which  simulates  human 
tissue.   The  scanning  detector  is  moved  through  the  water  at  various 
depths  to  determine  the  relative  percentage  depth  dose  at  those  depths. 


*Manufactured  by  General  Electric  Company,  Rectifier  Components 
Department,  Auburn,  New  York. 


150 


Y-COUNTERS 


OR 


O  -10  V 


1 


4020 


(512) 


14 


6  Out 


A  Y-COUNTER 


X 


Clock  14 


Data 


2  S7  9:2 


PE 


4018       *^l^3^^4J5 

Qi    Qz   Q3  Q4  Q5 


-O  -10  V 


^~m 


/6 
400 


12  I 


I 


4004 


(I6i 


6  Out 


Figure  25.   Diagrams  of  Y-Counter  and  AY-Count 
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The  water  phantom  serves  as  the  base  for  mounting  the  scanning  mechanism 
consisting  of  the  drive  motors  and  ball  screws.   Figure  26  shows  the 
water  phantom  and  the  mechanical  scanning  apparatus.   The  scanning 
ionization  chamber  is  shown  mounted  in  its  scanning  position  and  the 
reference  ionization  chamber  is  shown  as  it  would  be  mounted  during 
scanning. 

The  water  phantom  has  dimensions  of  approximately  56  x  56  x  56 
centimeters^  and  is  large  enough  for  measurement  of  a  35  x  35  centimeter 
field.   The  deepest  lateral  scan  is  performed  at  approximately  32 
centimeters  in  depth,  which  is  sufficient  to  reduce  the  depth  dose  to 
less  than  20  percent  with  cobalt  radiation.  The  radiation  beam  enters 
the  water  phantom  from  the  front  and  so  "depth"  refers  to  the  distance 
of  the  scanning  chamber  from  the  front. 

The  front  and  rear  panels  of  the  water  phantom  through  which  the 
radiation  enters  and  exits  are  made  of  a  transparent  plastic  material* 
to  facilitate  positioning  in  the  radiation  beam.   These  panels  are 
one-half  inch  thick  except  for  one  window  section,  through  which  the 
radiation  enters,  which  was  milled  out  to  a  thickness  of  about  3  milli- 
meters to  minimize  beam  attenuation.  The  3  millimeters  wall  thickness 
and  the  separation  of  the  ionization  chamber  from  the  window  is  equiva- 
lent to  a  total  thickness  of  about  5  millimeters  of  water.   This  window 
may  be  seen  in  Figure  26. 

The  moveable  base  on  which  the  water  phantom  rests  houses  the 
scanning  logic  and  the  amplifier  logic.   Together  the  base,  water  phantom 
and  mechanical  scanning  apparatus  make  up  the  scanning  unit  which,  when 


*An  acrylic  resin  generally  known  as  Lucite  or  Plexiglas.   Lucite 
is  a  trade  name  of  the  E.I.  DuPont  de  Nemours  and  Company,  Wilmington, 
Delaware.   Plexiglas  is  a  trade  name  of  the  Rohm  and  Haas  Company, 
Philadelphia,  Pennsylvania. 
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Figure  26.  Water  Phantom  and  Mechanical  S 


canning  Apparatus. 
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in  use,  is  located  in  the  radiation  therapy  room  and  is  connected  by 
way  of  a  cable  to  a  second  unit,  the  control  and  recording  unit,  which 
is  located  outside  the  therapy  room.' 

System  Testing 

The  performance  of  the  mechanical  scanning  system  was  evaluated 
by  two  methods.   For  the  first  method,  a  pulse  counter  was  attached  at 
each  motor  input  and  automatic  scanning  was  initiated.   The  pulse  count 
for  each  lateral  scan  in  both  directions  was  compared  for  all  scan  widths. 
This  test  was  performed  at  two  different  times,  the  first  time  when  the 
mechanical  scanning  system  was  completed  and  the  second  time  several 
weeks  later  after  the  system  had  accumulated  several  hours  of  run  time. 
At  both  times  the  number  of  pulses  in  each  lateral  scan  in  both 
directions  was  identical  for  all  scan  widths.   Likewise,  the  number  of 
pulses  for  each  depth  increment  was  identical  for  each  of  the  32  depth 
increments. 

For  the  second  test,  a  sharp  pointer  and  a  measuring  scale  were 
fixed  at  the  scanning  detector  mount  and  the  starting  position  of  the 
pointer  accurately  determined.   The  mechanical  scanning  system  was 
then  automatically  run  through  its  complete  cycle,  including  its  return 
to  the  original  starting  position,  for  about  four  hours,  during  which 
time  32  complete  scans  at  the  largest  width  were  performed.   The  position 
of  the  pointer  after  these  scans  was  identical  with  the  original  starting 
position  showing  that  no  pulses  were  gained  or  lost.   This  test  was 
later  repeated  at  several  different  times  and  the  results  were  identical. 
These  testing  procedures  show  the  absolute  reproducibility  of  the 
mechanical  scanning  system  and  ensure  the  proper  correlation  between 
detector  movement  and  the  measuring  cycles. 
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Adaptive  Scanning  Rate 

The  amplifier  logic  and  scanning  logic  modules  previously  discussed 
were  designed  so  that  the  detector  is  moved  the  same  distance  during 
each  amplifier  cycle.   Thus,  the  data  recorded  in  a  single  lateral 
scan  correspond  to  points  equally  spaced  across  the  radiation  field. 
Figure  4,  which  shows  the  response  of  a  previously  used  scanning 
system,  illustrates  the  results  obtained  when  measurements  at  equally 
spaced  points  are  used  to  approximate  the  depth  dose.   It  is  noted 
that  few  points  are  recorded  in  areas  of  rapid  signal  change,  but 
in  "flat"  areas  many  points  are  measured. 

An  adaptive  scanning  mode,  which  would  regulate  scanning  speed 
in  accordance  with  the  input  signal  changes,  would  make  it  possible 
to  measure  more  points  in  areas  of  rapid  signal  change  and  reduce  the 
total  scan  time  by  making  fewer  measurements  in  "flat"  areas. 

A  mechanical  scanning  system  which  is  driven  by  stepping  motors 
is  ideally  suited  to  an  adaptive  scanning  procedure  since  the  speed 
of  movement  depends  on  the  pulse  rate  to  the  motors.   An  adaptive 
scanning  mode  would  require  modification  of  the  scanning  logic  so  that 
the  pulse  rate  to  the  x-motor  would  vary  inversely  with  the  signal 
rate  of  change.   In  "flat"  areas  the  detector  would  move  at  maximum 
speed  recording  data  at  widely  spaced  points  and  in  areas  of  rapid 
signal  change  the  detector  would  slow  so  that  more  points  are  recorded 
at  shorter  intervals  apart.   With  an  adaptive  scanning  rate,  it  would 
no  longer  be  possible  to  correlate  detector  position  with  the  timing 
of  amplifier  cycles  as  is  the  case  when  equally  spaced  points  are 
measured.   The  recording  system  would  thus  require  modification  so  that 
the  cummulative  total  of  pulses  issued  to  the  x-motor  is  recorded  with 
the  data  for  each  amplifier  cycle. 
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The  Data  Processing  and  Recording  System 

The  repetitive  integral,  amplifier  logic  was  designed  so  that  two 
measurements  of  the  amplifier  output  voltage  are  taken  during  each 
measuring  cycle.   The  difference  of  these  two  measurements  is  then 
recorded  as  the  data  corresponding  to  the  point  at  the  center  of  the 
measuring  interval.   The  data  processing  and  recording  system  is  con- 
cerned with  the  processing  of  these  two  measured  values  beginning  at 
the  amplifier  output  and  until  the  value  of  the  difference  is  recorded. 
The  operation  of  the  data  processing  and  recording  system  may  be 
summarized  as  follows. 

The  analog  outputs  of  both  the  scanning  and  reference  detectors 
are  sent  to  separate  A/D  converters.   The  convert  command  for  the  first 
measurement  is  generated  at  the  appropriate  time  in  the  amplifier 
logic  and  is  sent  simultaneously  to  each  A/D.   (The  amplifier  logic 
is  discussed  in  Chapter  6  and  is  diagrammed  in  Figure  18.)   The  digital 
output  of  the  A/Ds  is  in  binary  form  and  is  loaded  in  parallel  mode 
into  shift  registers.   The  I's  complement*  of  these  data  is  then 
shifted,  bit-by-bit,  out  of  the  registers  in  syncronization  with  the 
clock  pulses.   The  bits  are  shifted  through  a  cable  from  the  scanning 
unit  to  the  recording  unit  into  an  identical  set  of  shift  registers, 
where  the  data  are  stored  until  the  second  measurement  of  the  cycle 


*The  I's  complement  of  a  binary  number  is  the  number  resulting 
after  all  1  bits  of  the  original  number  are  changed  to  0  and  all  0 
bits  are  changed  to  1.   The  2's  complement  is  the  I's  complement  plus  1. 
Given  two  binary  numbers  A  and  B,  the  difference,  B  less  A,  is  most 
easily  obtained  by  adding  B  to  the  2's  complement  of  A.   The  result  of 
this  addition  is  the  algebraic  difference  of  the  two  numbers. 
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is  made.   When  the  amplifier  logic  issues  the  convert  command  for 
the  second  measurement  the  above  process  is  repeated  except,  in  this' 
case,  the  converted  data  (not  the  complement)  are  shifted  into  a 
second  set  of  registers  at  the  recording  unit.   The  parallel  outputs 
of  the  two  sets  of  shift  registers  are  connected  to  arithmetic  units 
whose  output  provides  the  algebraic  difference  of  the  two  measurements. 
This  difference  is  recorded  on  IBM  computer-compatible  magnetic  tape. 

Figure  27  shows  a  complete  diagram  of  the  data  processing  and 
recording  system.   The  same  system  clock  which  drives  the  amplifier 
logic  and  scanning  logic  is  used  also  in  the  data  processing  and 
recording  system  to  time  and  coordinate  events.   The  specific  operation 
of  the  circuit  is  best  described  by  consideration  of  the  operation  of 
individual  components. 

The  Analog-To-Digital  Converters 

Identical  A/D  converters  are  used  at  the  outputs  of  the  scanning 
and  reference  amplifiers.   Both  units  are  the  Model  ADC-QM*  12-bit 
binary  converters  which  were  chosen  because  of  their  small  physical 
size  (2  X  4  X  0.4  inches^),  the  12-bit  resolution,  with  an  accuracy  of 
one-half  the  least  significant  bit,  and  the  short  conversion  time  of 
25  microseconds.   Analog  inputs  between  0  and  10  volts  are  accepted. 
The  12-bit  outputs  are  presented  in  parallel  mode. 

The  Shift  Registers 

The  shift  registers  are  of  the  RCA  COSMOS  4000  family  of  digital 


*Manufactured  by  Analog  Devices,  Inc.,  Norwood,  Massachusetts, 
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integrated  circuits.   Each  separate  unit  consists  of  8  stages  which 
are  capable  of  parallel  or  serial  input  and  output.   Separate  units  • 
may  be  cascaded  for  register  expansion  to  a  larger  number  of  bits. 
Storage  of  the  12  bits  from  each  A/D  requires  cascading  3  registers. 
The  most  useful  feature  of  the  shift  registers  is  the  facility  for 
parallel  input  and  serial  output.   The  serial  output  feature  allows 
the  use  of  just  one  data  line  to  the  recording  unit  rather  than  24 
lines  if  the  data  were  transmitted  in  the  parallel  mode. 

The  Optical  Couplers 

Optical  couplers  are  used  to  electrically  isolate  the  scanning 
and  recording  units  and  thus  eliminate  the  possibility  of  ground  loops. 
The  portion  of  the  circuit  shown  in  Figure  24,  which  is  to  the  left 
of  and  including  the  optical  couplers,  is  located  in  the  scanning  unit. 
The  remainder  of  the  circuit  is  located  in  the  recording  unit. 

Each  optical  coupler  (Model  MOC  1003*)  is  made  up  of  a  photo- 
transistor  and  an  infrared  emitting  diode. 

A  signal  applied  to  the  diode  in  the  coupler  causes  it  to  emit 
and  optically  transfer  a  pulse  to  the  phototransistor  which  sends  the 
pulse  on  through  the  connecting  cable.   The  voltages  to  operate  the 
diode  and  the  phototransistor  come  from  the  scanning  unit  power  supplies 
and  the  recording  unit  power  supplies,  respectively.   All  voltages 
from  the  recording  unit  are  identified  by  an  asterisk  in  Figure  27. 
All  data  and  logic  signals  between  the  two  units  go  over  an  optical 
coupler  except  the  start  and  stop  switches  which  are  just  interruptions 


*Manufactured  by  Motorola  Semiconductor  Products,  Inc.,  Phoenix, 
Arizona. 
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in  the  same  circuit  and  are  not  electrically  connected  to  the  re- 
cording unit. 

The  Arithmetic  Processing  Units 

The  two  shift  registers  which  store  the  two  measurements  are 
connected  to  6  cascaded  4-bit  binary  adders.   Each  adder  has  a  forward 
carry  feature  so  a  set  of  3  cascaded  adders  provides  a  true  binary  sum 
of  12  corresponding  bits  from  each  of  the  shift  registers. 

The  adders  provide  parallel  output  with  the  sums  settling  within 
530  nanoseconds  after  the  inputs  are  present.   These  devices  are  of 
the  RCA  COSMOS  4000  family  of  digital  integrated  circuits. 

Binary-To- Binary-Coded-Decimal  Converters 

The  output  of  the  binary  adders  is  in  binary,  but  the  digital 
module  used  to  display  the  difference  of  the  two  readings  requires 
binary-coded-decimal  (BCD)  input.   The  required  conversion  to  BCD  was 
done  using  Model  SN74185A*  binary-to-BCD  converters.   The  units  may 
be  cascaded  to  provide  conversion  for  a  larger  number  of  bits.   Eight 
of  the  devices  are  required  for  12-bit  conversion. 

The  Data  Processing  Circuit 

The  data  processing  circuit  shown,  in  Figure  27,  functions  in 
conjunction  with  the  amplifier  logic  circuit  shown  in  Figure  18,  both 
being  driven  by  the  system  clock.   The  clock  pulses  originate  at  the 
scanning  logic  (Figure  23)  and  come  directly  from  the  clock  at  full 


♦Manufactured  by  Texas  Instruments,  Inc.,  Dallas,  Texas. 
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frequency.   The  negative  clock  pulses  are  made  positive  by  an  oper- 
ational amplifier  so  that  logic  levels  of  0  and  10  volts  are  available. 
Because  of  the  several  logic  level  shifts  required  in  this  circuit, 
these  logic  levels  are  more  convenient.   The  clock  pulses  go  through 
a  NAM)  gate  and  into  a  7-stage  binary  counter,  which  reduces  the 
frequency  and  also  allows  clock  pulses  to  reach  the  circuit  only  when 
shifting  of  data  bits  is  desired. 

The  analog  outputs  of  both  the  scanning  (signal)  and  the  reference 
amplifiers  are  connected  to  the  input  of  the  A/Ds,  but  conversion  occurs 
only  on  command  from  the  amplifier  logic.   At  the  appropriate  time 
for  the  first  measurement,  a  convert  command  is  issued  by  the  amplifier 
logic  and  is  sent  to  the  S  input  of  the  convert  1-shot  shown  in 
Figure  27.   On  the  second  convert  command,  the  T  input  of  the  convert 
1-shot  is  energized.   The  output  of  the  1-shot  goes  high  for  20  micro- 
seconds triggering  both  A/D  units.   After  conversion,  the  digital  result 
is  available  in  parallel  binary  form  at  the  A/D  outputs.  The  logic 
level  of  the  A/D  outputs  is  converted  to  the  logic  level  of  the  shift 
registers  by  the  use  of  a  transistor  at  each  A/D  output.   These 
transistors  are  noted  as  level  shifters  in  Figure  27  and  convert  the 
data  bits  to  the  proper  logic  level,  but,  in  the  process,  invert  the 
bits.   This  means  the  I's  complement  of  the  data  is  loaded  into  the 
shift  registers. 

When  the  convert  command  is  issued  to  the  A/Ds,  the  convert- 
complete  output  of  the  signal  A/D  goes  high,  enabling  the  shift 
registers  to  accept  parallel  data  input  and  at  the  same  time  resetting 
the  clock  counter.   After  conversion,  this  output  goes  low  fixing  the 
A/D  outputs  in  the  registers  and  enabling  the  registers  to  operate 
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in  the  serial  data  mode.   The  clock  counter  innnediately  begins  to 
count  and  sends  pulses  to  the  clock  input  of  the  shift  registers 
and  over  an  optical  coupler,  to  the  shift  registers  in  the  recording 
unit.   When  in  the  serial  mode,  the  shift  registers  transfer  the  data, 
bit-by-bit,  synchronously  with  the  clock  pulses  from  the  scanning  unit 
to  the  recording  unit. 

Simultaneously  with  the  A/D  convert  pulse,  a  pulse  is  sent  to 
either  the  first  reading-enable  1-shot  or  the  second  reading-enable 
1-shot  (indicated  by  U  and  W  in  Figure  27),  according  to  which 
measurement  is  being  made.  These  pulses  cause  the  appropriate  1-shot 
to  send  a  pulse  to  a  NAND  gate  which  allows  clock  pulses  to  pass  to 
the  correct  shift  registers  at  the  recording  unit. 

The  clock  counter  disables  the  clock  pulses  after  24  pulses  are 
transmitted  to  the  shift  registers.   During  these  24  pulses,  all  24 
bits  from  the  A/D  conversion  are  transmitted  to  the  appropriate  shift 
registers  at  the  recording  unit.   The  reading-enable  pulse  from  the 
appropriate  1-shot  remains  high  for  about  4  milliseconds  which  is 
adequate  time  for  all  24  bits  to  be  shifted  into  the  registers. 

The  difference  of  the  two  measurements  is  determined  by  adding  the 
2's  complement  of  the  first  reading  to  the  second  reading.   The  first- 
reading  registers  contain  the  I's  complement,  and  the  2's  complement 
may  be  obtained  by  setting  the  carry  bit  high  on  the  least  significant 
binary  adder  for  each  set  of  12  bits.   The  binary  adders  then  sum  the 
I's  complement  of  the  first  reading,  plus  1  (carry  bit),  plus  the 
second  reading. 

The  parallel  outputs  of  the  adders  pass  through  level  buffers  to 
the  binary-to-BCD  converters  and  then  to  the  display  which  accepts  input 
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only  when  the  display  latch  is  opened.   When  the  second-reading  enable 
(W  in  Figure  27)  is  issued  by  the  amplifier  logic,  the  tape  trigger 
1-shot  goes  high  for  4  milliseconds.   When  this  pulse  goes  low,  a 
second  1-shot  is  triggered  actuating  the  tape  recording  unit  and 
opening  the  display  latch  for  70  microseconds.   This  1-shot  is  triggered 
only  after  both  readings  are  shifted  and  the  difference  is  available 
at  the  arithmetic  units.   The  display  latch  closes  again  after  70 
microseconds,  and  the  display  shows  the  latched  result  until  the  latch 
is  opened  again.   The  display  is  changed  during  each  measuring  cycle 
and  so  changes  8  times  each  second  which  is  too  rapid  for  manual 
recording.   Thus,  the  display  is  used  only  to  observe  the  operation 
of  the  system. 

The  Magnetic  Tape  Recorder 

The  magnetic  tape  recorder  (Model  1337*)  used  to  record  the 
difference  of  the  two  measurements  in  a  measurement  cycle  is  an 
incremental  digital  recorder  which  advances  one  step,  recording  8 
binary  digits  and  a  check  bit  at  each  trigger  pulse.   The  unit  is  a 
9-track  recorder  which  is  so  designed  that  the  spacing  between  each 
set  of  8  recorded  bits  is  IBM  compatible.   The  tapes  written  on  the 
recorder  may  be  read  with  any  IBM  computer  in  the  360/370  series  which 
is  equipped  with  9-track  tape  units.   The  recording  density  is  800 
bytes**  per  inch  of  tape  and  the  unit  may  record  at  up  to  300  steps 
per  second . 


*Manuf actured  by  Digi-Data  Corporation,  Bladensburg,  Maryland. 
**The  byte  is  the  basic  unit  of  data  and  memory  with  the  IBM 
360/370  series  computers.   The  byte  is  equivalent  to  8  binary  digits. 
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During  each  measurement  cycle,  the  difference  between  the  second 
reading  and  the  first  reading  for  both  the  scanning  and  reference 
amplifier  is  recorded.   Thus,  a  total  of  24  bits  is  recorded  during 
each  cycle,  12  for  each  amplifier.   Only  8  bits  are  recorded  at  each 
step  of  the  recorder  and  so  3  steps  are  required  to  record  the  results 
of  each  measurement  cycle. 

System  Testing 

The  only  testing  required  of  the  data  processing  and  recording 
system  was  to  verify  that  the  data  present  at  the  A/D  outputs  were 
transferred  to  the  storage  registers  and  that  the  difference  of  the 
two  readings  was  displayed  and  recorded.   A  constant  voltage  was 
applied  to  the  A/D  inputs  and  the  outputs  were  compared,  bit-by-bit, 
with  the  shift  register  outputs.   After  it  was  verified  that  the  data 
were  being  transferred  correctly,  a  ramp  voltage,  which  was  synchronized 
with  the  measuring  cycles,  was  connected  to  the  A/D  inputs.   The  constant 
difference  of  the  two  measurements  was  then  displayed  and  the  proper 
operation  of  the  arithmetic  units  could  be  verified  visually. 

The  Detector  and  Amplifier  System 

The  detector  and  amplifier  system  consists  of  two  ionization 
chambers  and  two  electrometer  amplifiers  which  operate  in  the  repeti- 
tive integral  mode  under  the  direction  of  the  amplifier  logic.   The 
amplifier  and  amplifier  logic  have  been  discussed  in  detail  in  Chapter 
6  and  are  mentioned  here  only  for  the  sake  of  completeness.   One 
ionization  chamber  and  amplifier  is  mounted  on  the  mechanical  scanning 
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apparatus  and  is  moved  through  the  water  phantom  during  scanning. 
The  other  ionization  chamber  and  amplifier  is  mounted  outside  the 
phantom  so  that  the  ionization  chamber,  which  serves  as  the  reference 
detector,  is  at  the  center  of  the  radiation  field.   Both  the  reference 
and  scanning  ionization  chambers  are  identical  in  construction  and  may 
be  interchanged.   Each  is  located  at  the  end  of  a  stylus  connected 
to  a  housing  containing  the  amplifier  (See  Figure  17).   This  separation 
of  ionization  chamber  and  amplifier  is  essential  so  that  the  amplifier 
is  outside  the  radiation  beam. 

The  reference  detector  is  used  when  measurements  are  made  with 
electron  accelerator  types  of  radiation  therapy  equipment.   With  this 
mode  of  x-ray  production,  there  is  continual  variation  of  radiation 
output  and  thus  measured  data  do  not  reflect  average  absorbed  dose, 
the  entity  which  should  be  represented  in  depth  dose  measurements. 
When  the  reference  chamber  is  used,  output  variations  cause  the  same 
percentage  change  in  the  reference  detector  signal  as  with  the  scanning 
detector.   If  the  ratio  of  the  two  signals  is  taken,  the  effect  of 
output  variation  is  eliminated.   The  reference  detector  is  not  needed 
when  measurements  are  made  with  isotope  therapy  units,  since  their 
outputs  do  not  vary. 

Detector  Construction 

The  construction  of  the  detector  assembly  is  important  since 
the  performance  of  the  entire  depth  dose  scanning  system  depends  upon 
the  basic  ionization  measurement.   Figure  28  shows  a  breakdown  of  the 
ionization  chamber  from  the  sensitive  volume  to  the  point  of  connection 
to  the  electrometer.   The  center  electrode,  which  collects  the  ions. 
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was  constructed  from  a  solid  aluminum  rod  and  was  highly  polished  to 
eliminate  any  air  cavities.   The  section  of  the  rod  which  projects 
inside  the  sensitive  volume  is  about  1.7  millimeters  in  diameter.   All 
but  this  section  of  the  rod  was  covered  with  insulating  tubing  which 
is  heat  shrinkable.   Heat  was  applied  to  the  tubing  so  that  it  shrunk 
tightly  to  the  electrode  eliminating  all  air  spaces  between  the  electrode 
and  the  insulating  tubing,  which  was  sprayed  with  conductive  paint.   A 
second  insulating  tubing  was  shrunk  into  tight  contact  with  the  con- 
ductive paint.   The  layer  of  conductive  paint  serves  as  a  guard 
electrode  and  was  connected  to  voltage  common.   The  second  tubing  was 
also  sprayed  with  conductive  paint  and  covered  with  a  third  tubing 
which  was  shrunk  tightly  to  the  paint.   This  second  layer  of  con- 
ductive paint  serves  as  the  high-voltage  electrode.   The  electrode 
assembly  was  then  placed  inside  a  stainless  steel  stylus  which  was 
connected  to  the  amplifier  housing. 

The  ionization  chamber  cap  was  made  from  a  phenolic  material 
which  is  nonconducting.   The  inside  walls  of  the  cap  were  painted  with 
a  carbon-base  conductive  paint  which  makes  contact  with  the  high-voltage 
electrode  when  the  cap  is  attached  to  the  electrode  assembly. 

The  ionization  chamber  center  electrode  was  connected  to  the 
electrometer  input  and  the  high-voltage  electrode  was  connected  to  a 
stabilized  160  volt  power  supply.   The  chamber  operates  in  the 
ionization  chamber  voltage  region  under  a  voltage,  which  is  sufficiently 
high  to  ensure  saturation  under  any  anticipated  irradiation  conditions 
but  is  not  so  high  as  to  initiate  ion  multiplication. 
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Design  Considerations 

The  ionization  current  produced  in  an  ionization  chamber  is 
directly  proportional  to  the  volume  of  the  chamber.   A  large  ionization 
current  is  very  desirable  in  depth  dose  scanning,  but  the  chamber 
volume  must  be  small  to  achieve  good  resolution.   As  a  compromise, 
the  ionization  chamber  was  made  cylindrical  in  shape  so  that  the 
diameter  of  the  chamber  moving  across  the  field  gradient  is  small. 
The  inside  diameter  of  the  chamber  is  3.8  millimeters,  and  the  chamber 
length  is  26.5  millimeters.   The  sensitive  volume  is  0.21  cubic 
centimeters.   The  length  of  the  center  electrode  inside  the  sensitive 
volume  is  20.8  millimeters.   The  solid  electrode  construction,  with 
its  large  diameter,  minimizes  microphonics  introduced  by  repeated 
acceleration  and  deacceleration  of  the  chamber  by  the  stepping  motors. 

System  Testing 

The  testing  of  the  repetitive  integral  electrometer  amplifier 
and  the  amplifier  logic  have  been  discussed  in  Chapter  6  and  will  not 
be  mentioned  further  here. 

Possibly  the  most  important  consideration  in  chamber  construction 
is  the  elimination  of  air  spaces  along  the  length  of  the  central 
electrode.   Any  air  spaces  present  provide  a  source  of  ionization  which 
contribute  to  the  ionization  current  from  the  chamber.   This  leads  to 
significant  stem  effect  and  makes  the  chamber  unsuitable  for  depth  dose 
scanning. 

The  stem  effect  of  a  radiation  detectot  should  be  less  than  0.5 
percent  when  used  in  depth  dose  scanning,  since  the  central-axis  depth 
dose  is  identical  with  rectangular  fields.   The  response  of  the 
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detector  must  be  identical,  to  within  0.5  percent,  no  matter  how  much 
of  the  detector  stylus  is  irradiated. 

The  scanning  detector  was  evaluated  for  stem  effect  by  measuring 
its  response  on  exposure  to  cobalt  radiation.   The  ionization  chamber 
was  fixed  in  position  and  was  not  moved  during  the  measurements.   The 
sensitive  volume  of  the  detector  was  placed  at  the  center  of  an 
8x8  centimeter  radiation  field,  and  the  difference  of  the  two 
measurements  made  in  each  measuring  interval  was  recorded.   The  field 
size  was  increased  to  8  x  12  centimeters,  with  the  12  centimeter  length 
along  the  stylus  connecting  the  sensitive  volume  to  the  amplifier 
housing,  and  the  difference  was  again  recorded.   The  collimator  Of 
the  cobalt  therapy  unit  was  then  rotated  90  degrees  so  that  the 
8  centimeter  side  of  the  field  was  along  the  stylus  and  the  difference 
recorded.   This  procedure  was  repeated  in  several  steps  up  to  a  field 
size  of  8  X  54  centimeters.   At  this  large  field,  the  entire  stylus 
was  irradiated  when  the  54  centimeter  side  was  along  the  stylus. 
Figure  29  shows  a  plot  of  the  results  obtained.   These  curves  show  that 
the  stem  effect  is  small  and  does  not  exceed  0.28  percent  and  thus  the 
chamber  is  satisfactory  for  scanning.   The  stem  effect  of  the  reference 
detector  was  later  evaluated  in  a  similar  manner  with  comparable  results, 

The  Total  System 

The  three  subsystems  composed  of  the  mechanical  scanning  system, 
the  data  processing  and  recording  system  and  the  detector  and  amplifier 
system  are  integrated  into  the  total  depth  dose  scanning  system  which  is 
housed  in  the  scanning  unit  and  the  recording  unit.   Figure  30  shows 
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a  diagram  of  the  total  system.   The  interconnection  of  the  three 
subsystems,  the  data  flow  and  the  logic  signal  flow  is  shovm.   The 
unit  location  of  each  component  part  is  also  included.   The  analog 
signals  from  each  detector  and  amplifier  to  the  A/D  imits  are  the 
only  analog  signals  utilized  in  the  system.   The  remainder  of  the 
depth  dose  scanning  system  operates  in  a  digital  mode. 

Figure  31  shows  the  scanning  unit  which  consists  of  the  mechanical 
scanning  apparatus,  the  water  phantom,  an  adjustable  table  on  which  the 
phantom  rests,  and  the  base  which  houses  the  amplifier  and  scanning 
logic  systems.   The  unit  has  wheels  for  easy  movement.   The  adjustable 
table  facilitates  fine  adjustments  during  scanner  setup.   The  table 
rotates,  moves  in  both  lateral  directions  and  may  be  raised  or  lowered. 

Figure  32  shows  the  control  and  recording  unit  which  houses  a 
portion  of  the  data  processing  and  recording  system  and  the  IBM- 
compatible  tape  recorder.   The  system  start  and  stop  switches  and  the 
display  units  are  on  the  front  panel  of  the  tape  controller  which  is 
shown  in  close-up  in  Figure  33.   Twelve  10-position  thumb  switches  on 
the  tape  controller  are  for  tape  identification  purposes  and  have  been 
coded  to  ensure  proper  identification  of  each  scan  recorded  on  the  tape. 
The  number  set  on  these  switches  may  be  manually  recorded  on  the  tape. 
Table  8  shows  the  coding  adopted. 

Computer  Processing 

A  computer  program  was  written  in  FORTRAN  IV  for  processing  the 
scan  data.   A  subroutine  was  written  in  IBM  370  assembly  language  to 
read  the  data  from  the  magnetic  tape  and  to  format  them  for  easy  access 
by  the  FORTRAN  language.   When  the  reference  chamber  is  not  used,  the 
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Figure  31,   Scanning  Unit. 
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Figure  32.   Control  and  Recording  Unit. 
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TABLE  8 
SCAN  IDENTIFICATION  SWITCH  CODING 

Switch  Position  Description 

1-2  Field  width  (cm) 

3-4  Field  length  (cm) 

5-7  SSD  (cm) 

8  Speed 

1  5.A18667  cm/sec 

2  2.709333  cm/sec 

3  1.354667  cm/sec 

4  0.677333  cm/sec 

5  0.338667  cm/sec 

9  Scan  width 

1  21.674667  cm 

2  27.093333  cm 

3  32.512  cm 

4  37.930667  cm 

10  Not  used 

11  Therapy  unit  scanned 

12  Scan  conditions 

0  with  trimmer 

1  no  trimmer 

2  wedged  field 
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data  recorded  from  the  scanning  chamber  are' normalized  to  100  percent 
at  the  depth  of  maximum  buildup.   Isodose  curves  are  then  plotted  on  a 
digital  plotter  under  computer  control  directly  from  the  measured  data 
with  no  smoothing  or  other  processing  techniques  being  used.   When  the 
reference  chamber  is  used,  the  ratio  of  the  scanning  chamber  output  to 
the  reference  chamber  output  is  determined  at  each  measurement  point. 
These  ratios  are  then  normalized  to  100  percent  and  isodose  curves  are 
plotted  as  above.   In  either  case,  the  data  used  to  plot  isodose  curves 
are  also  stored  for  later  use  in  treatment  planning. 

Summary 

The  depth  dose  scanning  system  may  be  considered  as  being  composed 
of  three  subsystems,  the  mechanical  scanning  system,  the  data  processing 
and  recording  system,  and  the  detector  and  amplifier  system,  which 
function  together  under  the  control  of  the  system  clock.   The  mechanical 
scanning  system  includes  the  water  phantom,  the  stepping  motors  and 
ball  screws  driving  the  detector  and  the  digital  logic  required  to 
direct  scanning  in  a  prescribed  pattern.   The  data  processing  and 
recording  system  is  concerned  with  processing  of  measured  data  from  the 
point  of  the  amplifier  output  to  the  point  of  recording  the  processed 
data.   The  difference  of  the  two  measurements  made  in  each  measurement 
cycle  is  determined  and  this  difference  is  recorded  on  IBM-compatible 
magnetic  tape.   The  detector  and  amplifier  system  consists  of  the 
scanning  and  reference  ionization  chambers,  their  associated  repetitive 
integral  amplifiers  and  the  logic  controlling  the  measurement  cycles. 
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A  large  digital  computer  (IBM  370)  is  used  to  process  the 
recorded  data  and  produce  isodose  curves  on  a  computer- controlled 
digital  plotter.   The  digital  data  are  preserved  for  later  use  in 
computer  treatment  planning. 

The  results  obtained  with  the  depth  dose  scanning  system  and 
conclusions  are  presented  in  the  next  chapter. 


CHAPTER  8 

RESULTS  AND  CONCLUSIONS 

Depth  Dose  Scanning  System  Results 

The  performance  of  the  depth  dose  scanning  system  was  evaluated 
by  run  testing  under  actual  scan  conditions  with  a  cobalt  teletherapy 
unit.   The  presence  or  absence  of  signal  lag  was  observed  by  comparison 
of  the  responses  obtained  when  scanning  across  the  radiation  field  in 
opposite  directions  at  the  same  depth.   To  facilitate  this  comparison, 
the  y-drive  motor  was  disconnected  so  that  no  depth  increments  could  be 
made.   A  10  x  10  centimeter  field  was  set  on  the  therapy  unit  and  the 
scanner  was  set  up  for  automatic  scanning.   The  detector  moved  across 
the  field  32  times  at  the  same  depth,  16  times  in  each  direction,  at 
the  maximum  scanning  speed  of  2.71  centimeters  per  second.  The  measured 
data  were  plotted  at  5  to  90  percent  isodose  levels  as  if  each  lateral 
scan  was  separated  by  a  depth  increment  of  1.016  centimeters.   The 
results,  given  in  Figure  34,  consist  of  a  series  of  parallel  lines 
and  show  that  the  responses  obtained  with  the  16  sets  of  scans  in 
opposite  directions  do  not  depend  on  the  direction  of  travel.   Hence, 
there  is  no  effect  of  signal  lag  in  the  measurements.   This  result  is 
of  interest  when  compared  with  Figure  4  in  which  the  rate  measurement 
results  for  scans  in  opposite  directions  show  a  deviation  of  about 
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Figure  34.   Plot  of  32  Lateral  Scans  at  Same  Depth, 
16  in  each  Direction. 
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4  millimeters  even  at  a  much  slower  scanning  speed  (about  1.0  centimeter 
per  second) . 

The  reproducibility  of  the  depth  dose  scanning  system  was 
evaluated  under  actual  scan  conditions.   A  10  x  10  centimeter  field 
was  scanned  under  automatic  control,  this  time  with  all  depth  incre- 
ments, at  maximum  scanning  speed.   Neither  the  scanning  unit  nor  the 
therapy  unit  was  moved  and  the  same  radiation  field  was  completely 
scanned  automatically  4  more  times  for  a  total  of  5  scans  for  the  same 
field.   The  first  lateral  scan  of  each  was  plotted  as  shown  in  Figure 

35.  The  maximum  deviation  of  any  set  of  5  data  values  corresponding 
to  the  same  point  on  the  scan  was  0.8  percent.   The  data  corresponding 
to  the  central  axis  of  each  of  the  5  scanned  fields  were  plotted  to 
observe  the  reproducibility  of  the  depth  dose  and  are  given  in  Figure 

36.  The  points  are  almost  coincident  and  the  maximum  deviation  of  any 
set  of  5  data  points  corresponding  to  the  same  depth  was  0.6  percent. 
A  set  of  11  isodose  curves  from  5  to  100  percent  was  plotted  from  the 
data  for  each  for  the  5  scanned  fields  and  the  curves  were  superimposed 
upon  one  another  as  shown  in  Figure  37.   All  5  sets  of  isodose  curves 
show  excellent  reproducibility. 

The  accuracy  of  a  set  of  isodose  curves  is  difficult  to  determine 
because  of  a  lack  of  data  with  which  to  compare.   Standard  central-axis 
depth  dose  tables,  such  as  those  found  in  the  British  Journal  of 
Radiology,  Supplement  11,'^'^  have  been  derived  by  combining  the 
measured  results  of  many  researchers  for  many  different  types  of  cobalt 
therapy  units  and  are  generally  accepted  as  representative  of  most 
cobalt  therapy  units.   The  central-axis  depth  dose  data  plotted  in 
Figure  36  were  compared  with  the  standard  tables  and  the  maximum 
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Figure  37.   Five  Complete  Isodose  Distributions  for  a  10  x  10 
Centimeter  Field  Superimposed. 


deviation  at  any  depth  was  found  to  be  1.1  percent.   No  standard 
tables  are  available  for  comparison  in  the  penumbra  region  of  the  field, 
where  the  depth  dose  changes  very  rapidly,  because  the  depth  dose  in 
this  region  is  dependent  upon  the  geometry  of  the  beam  collimating 
device  used  with  the  therapy  unit.   It  can,  however,  be  concluded  that 
the  results  obtained  in  this  region  are  also  accurate  to  better  than 
1  percent  because  (1)  the  theoretical  response,  shown  in  Figure  10 
and  Table  2,  of  the  repetitive  integral  measurement  method  to  the 
lateral  dose  function  was  almost  identical  in  the  region  of  rapid 
change  and  (2)  the  measured  data  plotted  in  Figure  3A  show  no  dependence 
on  the  direction  of  travel. 

The  depth  dose  scanning  system  has  been  used  thus  far  to  measure 
approximately  100  different  radiation  fields  with  2  cobalt  units.   The 
system  has  performed  in  a  flawless  manner  and  the  data  from  only  one 
radiation  field  were  unusable  because  of  a  bad  spot  on  the  magnetic 
tape  used  to  record  the  data.   Figures  38  -  40  show  examples  of  the 
plots  resulting  from  data  measured  with  the  system.   Figure  38  is  the 
plot  of  a  10  X  10  centimeter  field  with  a  different  therapy  unit  from 
the  one  used  to  record  the  data  plotted  in  Figure  37.   Figure  40  is 
the  plot  of  a  wedged  field  where  a  light  gap  of  1  centimeter  was  used 
with  the  wedge.   In  this  case,  the  wedged  field  was  normalized  to 
100  percent  for  the  same  size  open  field. 

Conclusions 

The   goals   set   forth   at    the  beginning  of   this   project  were   the 
development  of   a   fast,    accurate,    self-contained  and  easily   transportable 
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Figure   38.      Isodose  Curves   for  10  x  10  Centimeter  Cobalt 
Radiation  Field. 
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Figure   39.      Isodose   Curves    for  5  x  10  Centimeter  Cobalt 
Radiation  Field. 
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Figure  40.   Isodose  Curves  for  4x4  Centimeter  Cobalt 
Radiation  '..'edged  Field. 
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depth  dose  scanning  system.   In  order  to  meet  these  goals,  it  was 
necessary  to  develop  a  new  measurement  technique  and  a  new  detector 
transport  system. 

A  new  measurement  technique  called  the  repetitive  integral  method 
has  been  shown  theoretically  to  be  superior  to  the  rate  method  for  the 
real-time  dynamic  current  measurements  required  in  depth  dose  scanning. 
Electronic  instrumentation  which  functions  in  the  repetitive  integral 
n\ode  has  been  developed  and  experimental  results  have  been  shown  to 
follow  the  theory.   The  repetitive  integral  measurement  method  applied 
to  depth  dose  scanning  has  the  following  advantages: 

1.  The  effect  of  the  amplifier  time  constant  is  eliminated.' 

2.  The  method  is  better  suited  for  measurement  of  pulsed 
radiation  sources. 

3.  Long-term  amplifier  drift  may  be  eliminated. 

4.  The  integrating  circuit  acts  to  attenuate  high-frequency 
noise. 

5.  The  error  in  the  data  points  caused  by  thermal  noise  was 
found  to  be  negligible. 

6.  The  method  is  independent  of  scanning  speed. 
Application  of  the  repetitive  integral  measurement  method  in 

conjunction  with  the  new  detector  transport  system  has  made  possible 
accurate  and  reproducible  measurements  at  scanning  speeds  greater  than 
previously  possible.   Speed  of  the  depth  dose  scanning  system  is 
currently  limited  by  the  mechanical  scanning  apparatus  rather  than  the 
amplifier  response  and  may  be  further  increased  by  using  (1)  larger 
scanning  motors  to  drive  the  scanning  assemblies,  or  (2)  additional 
ionization  chambers  so  that  the  total  scan  width  is  reduced,  or  (3)  an 
adaptive  scanning  rate.   Some  modification  of  the  scanning  logic  and 
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the  recording  system  would  be  required  to  implement  the  last  two 
above  possibilities. 

The  digital  techniques  employed  in  the  mechanical  scanning  system, 
the  scanning  logic  and  the  amplifier  logic  provide  the  means  of 
guaranteeing  the  accuracy  and  reproducibility  of  the  depth  dose 
scanning  system.   Such  techniques  permit  precise  positioning  of  the 
radiation  detector  to  within  a  small  fraction  of  a  millimeter  and 
e^nsure  correlation  of  detector  position  with  measurement  cycles. 

The  depth  dose  scanning  system  is  easy  to  set  up  for  scanning 
since  fine  adjustment  of  the  height  and  phantom  position  may  be  made 
after  the  phantom  is  filled  with  water.   The  entire  system  is  completely 
self-contained  and  is  housed  in  two  units,  each  of  which  has  wheels 
for  easy  movement  within  a  working  area.   The  system  is  small  enough 
to  be  easily  transportable  using  a  small  van. 

All  the  goals  set  forth  for  the  development  of  a  depth  dose 
scanning  system  have  been  realized.  The  system  developed  provides  a 
means  of  rapidly  determining  complete  depth  dose  distribution  for  a 
variety  of  radiation  therapy  equipment. 
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